上机作业

系统仿真上机作业姓名：屈佳龙学号：04103176班级：041032系统仿真上机作业一、 计算机辅助系统分析：r ：单位阶跃N G ：非线性器件)1025.0)(1625.0)(110()1()(0++++=s s s s s K s G1．1=N G ，40=K 时，用MATLAB 画出开环Bode 图，求出c ω、βθ。

由其估计出r t 、s t 、%δ答：>> den=conv(conv(conv([1 0],[10 1]),[0.625 1]),[0.025 1]); >> num=[40 40]; >> sys=tf(num,den); >> sys=feedback(sys,1); >> [y,t]=step(sys); >> ytr=find(y>=1); >> [ym,tp]=max(y); >> rise_time=t(ytr(1)) rise_time = 0.7772（上升时间） >>peak_time=t(tp) peak_time = 1.2953 (峰值时间) max_overshoot=ym-1 max_overshoot = 0.8168（超调量） s=length(t);while y(s)>0.98&y(s)<1.02 s=s-1; endsettling_time=t(s+1) settling_time = 17.2280（调整时间） >>bode(num,den)[mag,phase,w]=bode(num,den);[gm,pm,wcg,wcp]=margin(mag,phase,w) gm = 4.3168 （增益裕量） pm = 10.0158 （相位裕量） wcg = 5.1598 （交界频率）yu+ - + )(s G N)(0s Gr ewcp = 2.3975 （增益=0dB 时） Bode 图如下：可得：c ω=5.1598 βθ=10.0158r t =0.7772s s t =17.0985s %δ=81.68%2．1=N G ，40=K 时，用MATLAB 画出根轨迹图，并求出40=K 时的闭环极点；由其估计出r t 、s t 、%δ 答：绘制系统的根轨迹：num=[0 0 0 40 40];den1=conv([1 0],[10 1]);den2=conv([0.625 1],[0.025 1]); den=conv(den1,den2); sys=tf(num,den); [p,z]=pzmap(sys) rlocus(sys)根轨迹图形：系统闭环零极点:程序： sys=feedback(sys,1);[p,z]=pzmap(sys)极点：p = -40.1616-0.2274 + 2.4146i -0.2274 - 2.4146i -1.0837零点：z = -1求r t ，s t ，%δ：根据零极点的特点，取主导极点p=-0.2274+2.4146i 和p=-0.2274-2.4146i ，可以将此系统化为二阶系统个G （s ）=sys =有以下程序可求出r t 、s t 、%δ 程序：num=[0 0 5.882]; den=[1 0.4548 0]; sys=tf(num,den)sys=feedback(sys,1); [y,t]=step(sys); ytr=find(y>=1);rise_time=t(ytr(1)) [ymax,tp]=max(y); peak_time=t(tp);max_overshoot=ymax-1 s=length(t);while y(s)>0.98&y(s)<1.02 s=s-1; endsettling_time=t(s)结果：rise_time =0.7772max_overshoot = 0.7438settling_time = 16.9691可得：r t =0.7772s s t =16.9691s %δ=74.38%3．1=N G ，40=K ：仿真之，并由仿真结果求出r t 、s t 、%δ 答：建立系统模型用自适应变步长方法（ode45）仿真可得：观察上图可知：超调量约为80%、上升时间Tr 约为0.5S ；调整时间Ts 约为18S 。

Word 上机题1、 将下列文字复制3次，形成4个内容相同的段落，并为所有的段落同时加上“●”形的项目符号，不得作其他修改。

计算机是新技术革命的一支主力，也是推动社会向现代化迈进的活跃因素。

计算机科学与技术是第二次世界大战以来发展最快、影响最为深远的新兴学科之一。

2、 使用“段落”对话框，将下列段落的格式设定为：首行缩进1厘米、行距为2倍行距，并拷贝2次，成为3个段落。

要求最后一个回车符在在最后一个段落的尾部。

现代计算机是一种按程度自动进行信息处理的通用工具。

它处理的对象是信息，处理结果也是信息。

在这一方面上，计算机和人脑有某些相似之处。

因为人的大脑和五官也是采集、识别、转换、储存和处理信息的器官，所以，有许多人把计算机称为电脑。

3、 将下列文字复制5次，形成6个内容相同的段落，然后将前3个段落连成一个段落，其余的段落连成一个段落，不得作其他修改。

最后一个回车符在在最后一个段落的尾部。

计算机是新技术革命的一支主力，也是推动社会向现代化迈进的活跃因素。

计算机科学与技术是第二次世界大战以来发展最快、影响最为深远的新兴学科之一。

4、 将下列文字复制3次，并连成一个段落（文档中仅有一个回车符），然后按“居中”格式排版，不得作其他任何修改。

现代计算机是一种按程度自动进行信息处理的通用工具。

它处理的对象是信息，处理结果也是信息。

在这一方面上，计算机和人脑有某些相似之处。

因为人的大脑和五官也是采集、识别、转换、储存和处理信息的器官，所以，有许多人把计算机称为电脑。

5、 将下列文字复制3次，每次复制的内容各成一个段落，使文档共有4个内容相同的段落，将各段落的字体依次设置为宋体、仿宋体、楷体和黑体，不得作其他修改。

计算机是新技术革命的一支主力，也是推动社会向现代化迈进的活跃因素。

计算机科学与技术是第二次世界大战以来发展最快、影响最为深远的新兴学科之一。

6、7、 下列2行4列表格的行高是15磅，列宽是3厘米，表格修改成同样格式的4行4列的8、 使用“拆分单元格”对话框将表格二修改成表格一的式样。

计算方法上机作业上机实习题目1.某通信公司在一次施工中，需要在水面宽度为20米的河沟底部沿直线走向铺设一条沟底光缆。

在铺设光缆之前需要对沟底的地形进行初步探测，从而估计所需光缆的长度，为工程预算提供依据。

已探测到一组等分点位置的深度数据(单位:米)如下表所示：分点0 1 2 3 4 5 6深度9.01 8.96 7.96 7.97 8.02 9.05 10.1 3分点7 8 9 10 11 12 13深度11.1812.2613.2813.3212.6111.2910.22分点14 15 16 17 18 19 20深度9.15 7.90 7.95 8.86 9.81 10.810.93(1)请用合适的曲线拟合所测数据点；(2)估算所需光缆长度的近似值，并作出铺设河底光缆的曲线图；（1）算法思想分段多项式是由一些在相互连接的区间上的不同多项式连接而成的一条连续曲线，其中三次样条插值方法是一种具有较好“光滑性”的分段插值方法。

在本题中，假设所铺设的光缆足够柔软，在铺设过程中光缆触地走势光滑，紧贴地面，并且忽略水流对光缆的冲击。

计算光缆长度时，用如下公式：20()L f x ds=⎰200(f x =⎰191(k kk f x +==∑⎰= 本题采取三次样条插值的方法，因为三次样条插值方法是一种具有较好“光滑性”的分段插值方法。

根据提供的数据，只用x,y 值，不包含导数值，因此采用第三类三次插值多项式进行插值编程。

设计算法如下： 1.For n i ,,2,1,0⋅⋅⋅=1.1 i i M y ⇒2.For 2,1=k2.1 For k n n i ,,1, -=2.1.1 i k i i i i M x x M M ⇒----)/()(13. 101h x x ⇒-4.For 1-,,2,1n i =4.1 11++⇒-i i i h x x4.2ba c c h h h i i i i i i ⇒⇒-⇒+++2;1;)/(114.3 i i d M ⇒+165.000;;c M d M d n n ⇒⇒⇒λnn n b a b ⇒⇒⇒2;;20μ6. 1111,γμ⇒⇒d b7.获取M 的矩阵元素个数，存入m8.For m k ,,3,2 =8.1kk k l a ⇒-1/μ8.2 k k k k c l b μ⇒⋅-1- 8.3 kk k k l d γγ⇒⋅-1-9.m m m M ⇒μγ/10. For 1,,2,1 --=m m k 10.1 k k k k k M M c ⇒⋅-+μγ/)(1 11. 获取x 的元素个数存入s 12. k ⇒113. For 1,,2,1-=s i 13.1 ifix x ≤~ then k i ⇒；breakelse k i ⇒+114.xx x x x x h x x k k k k ˆ~;~;11⇒-⇒-⇒---yh x h M y x h M y x M x M k k k k k k ~/]ˆ)6()6(6ˆ6[2211331⇒-+-++---（3）源程序clear; clc;x=0:1:20; %产生从0到20含21个等分点的数组 X=0:0.2:20;y=[9.01,8.96,7.96,7.97,8.02,9.05,10.13,11.18,12.26,13.28,13.32,12.61,11.29,10.22,9.15,7.90,7.95,8.86,9.81,10.80,10.93]; %等分点位置的深度数据n=length(x); %等分点的数目 N=length(X);%% 求三次样条插值函数s(x)M=y;for k=2:3; %计算二阶差商并存放在M中 for i=n:-1:k;M(i)=(M(i)-M(i-1))/(x(i)-x(i-k+1));endendh(1)=x(2)-x(1); %计算三对角阵系数a,b,c及右端向量d for i=2:n-1;h(i)=x(i+1)-x(i);c(i)=h(i)/(h(i)+h(i-1));a(i)=1-c(i);b(i)=2;d(i)=6*M(i+1);endM(1)=0; %选择自然边界条件M(n)=0;b(1)=2;b(n)=2;c(1)=0;a(n)=0;d(1)=0;d(n)=0;u(1)=b(1); %对三对角阵进行LU分解y1(1)=d(1);for k=2:n;l(k)=a(k)/u(k-1);u(k)=b(k)-l(k)*c(k-1);y1(k)=d(k)-l(k)*y1(k-1);endM(n)=y1(n)/u(n); %追赶法求解样条参数M(i)for k=n-1:-1:1;M(k)=(y1(k)-c(k)*M(k+1))/u(k);ends=zeros(1,N);for m=1:N;k=1;for i=2:n-1if X(m)<=x(i);k=i-1;break;elsek=i;endendH=x(k+1)-x(k); %在各区间用三次样条插值函数计算X点处的值x1=x(k+1)-X(m);x2=X(m)-x(k); s(m)=(M(k)*(x1^3)/6+M(k+1)*(x2^3)/6+(y(k)-(M(k)*(H^2)/6))*x1+(y(k+1)-(M(k+1)*(H^2)/6))*x2)/H;end%% 计算所需光缆长度L=0; %计算所需光缆长度for i=2:NL=L+sqrt((X(i)-X(i-1))^2+(s(i)-s(i-1))^2);enddisp('所需光缆长度为 L=');disp(L);figureplot(x,y,'*',X,s,'-') %绘制铺设河底光缆的曲线图xlabel('位置X/测量点','fontsize',16); %标注坐标轴含义ylabel('深度Y/m','fontsize',16);title('铺设河底光缆的曲线图','fontsize',16);grid;（4）结果与分析铺设海底光缆的曲线图如下图所示：拟合结果表明，运用分段三次样条插值所得的拟合曲线能较准确地反映铺设光缆的走势图。

常 出现 的错 误 。 空 经 常 翻 看 逐 步 积 累经 验 ． 免 自 己 抽 避
论 学得好 ． 并不 代表 上机实践 就好 ， 事实 上 ， 上机 实 在 训 教学 中有很 多学 生都完成不 了上机作业 ．在这 里主
要从教学 方法的角度上来分析原 因
程序代 码违背 了 Ｃ语言 的语法规 定而 出现错误 ，
编译 程序一般会用英文给 出“ 出错 信 息 ”１ 示 。 告 ［ ２ 的提 并 诉 错 误 在 哪 一行 由于 学 生 英 语 水 平 普 遍 不 高 ． 难 明 很
：
《 数据结构》 上机作业难完成 的成 因分析
杨智 明 ． 漆 丽 娟
６７０ ） ５０ ０
／ ／ ／
（ ． 南保 山学 院 计 算 机 科 学 系 ， 山 ６ ８ ０ ； ． １云 保 ７０ ０ ２ 云南 昭通 师 专 计算 机 科 学 系 ， 通 昭
摘
要 ：通 过 对保 山 学 院 计科 系 ２ ０ ￣ ０ ０届 的 学 生难 以 完成 《 据 结 构 》 机 作 业 的情 况分 析 ， ０５２１ 数 上 从 教 学 方 法 的 角度 寻找 原 因 ， 据 实 际情 况提 出 上 机 实训 教 学 改 革 的 具 体 措 施 ， 根 以促 使 学 生 顺 利 完成 作 业 ． 高 用计 算 机 解 决 问题 的 能 力 。 提
１１ 算 法 转换 成 程 序 困难 ．
在 我 们 选 用 的严 蔚 敏 编 著 Ｃ版 本 的《 据 结 构 》 数 教
白编译程序提示语句错误 的原因 ．或者编译程序 提示 的语 句错误不属 于本行 ． 学生查错 、 使 改错存在 一定的
困难 。
材 中 ． 多经典 算法都 采用类 Ｃ语 言描述［ 要想 把算 许 １ ］ 。