《MATLAB及应用》参考答案

《MATLAB及应用》上机作业学院名称：机械工程学院专业班级：测控1201学生姓名：学生学号：201年4月

《MATLAB及应用》上机作业要求及规范一、作业提交方式：word文档打印后提交。二、作业要求：1．封面：按要求填写学院、班级、姓名、学号，不要改变封面原有字体及大小。2．内容：只需解答过程（结果为图形输出的可加上图形输出结果），不需原题目；为便于批阅，题与题之间应空出一行；每题答案只需直接将调试正确后的M文件内容复制到word中（不要更改字体及大小），如下所示：%作业1_1clcA=[1234;2357;1357;3239;1894];B=[1+4*i43677;233554+2*i;26+7*i5342;189543];C=A*BD=C(4:5,4:6)三、大作业评分标准：1．提交的打印文档是否符合要求；2．作业题的解答过程是否完整和正确；3．答辩过程中阐述是否清楚，问题是否回答正确；4．作业应独立完成，严禁直接拷贝别人的电子文档，发现雷同者都以无成绩论处。

作业11、用MATLAB可以识别的格式输入下面两个矩阵，再求出它们的乘积矩阵，并将矩阵的右下角子矩阵赋给矩阵。赋值完成后，调用相应的命令查看MATLAB工作空间的占有情况。解：A=[1234;2357;1357;3239;1894;]B=[1+4i43677;233554+2i;26+7i5342;189543;]B=[1+4i43677;233554+2i;26+7i5342;189543;]C=A*BD=C(4:5,4:6);2、设矩阵，求，，，，，并求矩阵的特征值和特征向量。解：A=[162313;511108;97612;414152;]det(A)inv(A)A.^32*A+inv(A)3*A-A'[V,D]=eig(A)abs(A)3、解下列矩阵方程:解：A=[010;100;001;];B=[100;001;010;];C=[1-43;20-1;1-20;];X=inv(A)*C*inv(B)4、求多项式当时的值，并求的导数。解：a1=[1

5

35

13]

a2=[0

1

45

3]

a3=[0

0

31

3]

p1=conv(a1,a2)

p1=conv(p1,a3)

polyval(p1,3)

polyder(p1)5、求多项式的根和导数。解:

p1=[1

-5

-14

-10

-3]

roots(p1)

polyder(p1)6、对于有理多项式(1)计算该多项式相除的结果;(2)将该多项式展开为部分分式的形式;(3)计算。解(1)A=[104567];B=A.*10;C=[11];D=[12];E=[13];F=conv(conv(C,D),E);G=deconv(B,F);(2)[r,p,k]=residue(B,F)(3)[p,q]=ployder(B,F)7、在某次传感器输入输出特性实验中测得输入输出的一组数据如下表所示:（输入）12345（输出）1。31。82。22。93。5已知输入和输出可以近似成线性关系，即，求系数和，并求当输入时输出的值。解x=[12345];y=[1.31.82.22.93.5];p=polyfit(x,y,1)a=polyval(p,8)8、根据人口理论的马尔萨斯模型可知，人口总数可以采用指数函数对人口数据进行拟合。据统计，六十年代世界人口数据如下(单位：亿)t196019611962196319641965196619671968y3。39183。42133。45033。46983。47633。49203。51333。53223。5505试求马尔萨斯模型中的a，b值，并画出拟合曲线图，同时预测一下2010年的人口数值。解year=[1960:1:1968];n=[3.39183.42133.45033.46983.47633.49203.51333.53223.5505];y=log(n);p=polyfit(year,y,1);plot(year,y,'g*');a=p(2)b=p(1)y2010=exp(polyval(p,2010))9、某实验测得强度随时间变化的一组数据:00。511。522。5300。47940。84150。99150。90930。59850。14111)利用二次曲线拟合求出秒处强度指标。2)利用样条曲线插值求出秒处强度指标。解：

t=[0

0.5

1

1.5

2

2.5

3]

y=[0

0.4794

0.8415

0.9915

0.9093

0.5985

0.1411]

p1=polyfit(t,y,2)

y1=polyval(p1,2.25)

y2=interp1(t,y,2.25,'spline')试用MATLAB求解下面的问题:(1)求极限；(2)不定积分；(3)对（2）的不定积分结果进行微分，看是否能还原原函数；(4)对函数做20项的Taylor幂级数展开；(5)求级数的和。解：（1）f=((sin(x))^2-(sin(a))^2)/(x-a);limit(f,x,a)（2）f=[x^2*e^(2*x)+sinx*e^x];int(f)(3)x=sym('x');

f=1/3*x^3+2*x+log(x^2+1)+3*atan(x);

diff(f)(4)f=x^2*sin(cos(x^2))*cos(x);taylor(f,x,21)(5)S=symsum(n^2,n,1,inf)S=symsum(n^2,n,1,100)11、求解下列方程组的解(1)(2)解（1）[xy]=solve('sin(x+y)-exp(x)*y=0','x^2-cos(y)=0','x,y')（2）[xy]=solve('x^2+x*y+y=0','x^2-4*x+3=0','x,y')12、求微分方程(组)的解(1)，，(2)，(3)解：(1)y=dsolve('D3y+Dy-x','y(1)=8','Dy(1)=7','Dy(2)=4','x')(2)x=dsolve('(Dx)^2+x^2-1','x(0)=0')(3)[x,y]=dsolve('Dy-2*x+y','Dx-4*x+2*y','t')作业21、一球从100米高度自由落下，每次落地后反跳回原高度的一半，再落下。求它在第10次落地时，共经过多少米?第10次反弹有多高?解：y=0.5*x;sum=0;fori=2:10,sum=sum+x+y;x=0.5*x;y=0.5*x;endsumy2、用MATLAB的M函数文件定义如下分段函数：解：functiony=f(x)ifx>10y=5;elseifx>=-10&x<=10y=1/2*x;elsey=-5;end3、分别用for和while循环编写程序，求出并考虑一种避免循环的简洁方法来进行求和，并比较各种算法的运行时间。解：x=2;K=0;fori=1:63,K=K+x^i;endK4、应用MATLAB语言及二分法编写求解一元方程在区间[3，6]的实数解的算法，要求绝对误差不超过0.001。解：A=[1-1459-70];a=3;b=6;c=0.001;whilec<0.5*(b-a),x1=(a+b)/2;f1=polyval(A,x1);fa=polyval(A,a);fb=polyval(A,b);iff1==0;x=x1;elseiff1*fa>0;a=x1;elsef1*fb>0;b=x1;endendendx=x15、二阶系统的单位阶跃响应为，在同一平面绘制分别为0，0。3，0。5，0。707的单位阶跃响应曲线。要求：四条曲线的颜色分别为蓝、绿、红、黄，线型分别为“——”、“……”、“oooooo”、“++++++”；(2)添加横坐标轴和纵坐标轴名分别为“时间t”和“响应y”，并在平面图上添加标题“二阶系统曲线”和网格；(3)在右上角添加图例（即用对应的字符串区分图形上的线），并分别在对应的响应曲线的第一个峰值处标示“zeta＝0”、“zeta＝0。3”、“zeta＝0。5”、“zeta＝0。707”。解：t=[0:0.1:10];kos=0;y=1-1/sqrt(1-kos^2)*exp(-kos*t).*sin(sqrt(1-kos^2)*t+2*cos(kos));plot(t,y,'b-')holdon;kos=0.3;y=1-1/sqrt(1-kos^2)*exp(-kos*t).*sin(sqrt(1-kos^2)*t+2*cos(kos));plot(t,y,'g.')holdon;kos=0.5;y=1-1/sqrt(1-kos^2)*exp(-kos*t).*sin(sqrt(1-kos^2)*t+2*cos(kos));plot(t,y,'ro')holdon;kos=0.707;y=1-1/sqrt(1-kos^2)*exp(-kos*t).*sin(sqrt(1-kos^2)*t+2*cos(kos));plot(t,y,'y+')holdon;xlabel('时间t');ylabel('响应y');gridon;title('二阶系统曲线');legend('kos=0','kos=0.3','kos=0.5','kos=0.707')gtext('kos=0')gtext('kos=0.3')gtext('kos=0.5')gtext('kos=0.707')6、绘制如下图所示的图形，把图形窗口分割为2行2列，窗口左上角画一正弦曲线；窗口右上角画3条单边指数曲线；窗口左下角画一矩形脉冲信号，脉冲宽度为1，高为2，开始时间为1；窗口右下角画一单位圆。解：t=0:0.01:2;subplot(2,2,1);plot(t,sin(2*pi*t));title('plot(x,2*PI*t)');gridon;subplot(2,2,2);plot(t,exp(-t));holdon;plot(t,exp(-2*t),'g');holdon;plot(t,exp(-3*t),'r');holdon;title('exp(-t),exp(-2*t),exp(-3*t)');gridon;x=[0112234];y=[0022000];subplot(2,2,3);plot(x,y);axis([04-0.53]);title('pulsesignal');symsxysubplot(2,2,4);ezplot(x^2+y^2-1);axis([-1.51.5-1.21.2]);title('circle');xlabel('')ylabel('')7、已知函数，试分别应用三维曲线图绘制命令plot3、三维网线图绘制命令mesh、三维曲面图绘制命令surf在同一窗口中绘制出3个子图。解：x=-2:0.1:2;[X,Y]=meshgrid(x);Z=1./sqrt((1-X).^2+Y.^2)+1./sqrt((1+X).^2+Y.^2);subplot(311);plot3(X,Y,Z);gridon;axis([-11-11-19]);subplot(312);mesh(X,Y,Z);axis([-11-11-19]);subplot(313);surf(X,Y,Z);axis([-11-11-19]);8、求解下面两个方程构成的联立方程组在区间内的解，并用绘图的方法绘出两曲线在同一坐标上的图，以验证求得的解的正确性。，解：clc[x,y]=solve('(x-2)^2+3*(y-1)^2-3=0','y-6*(x-2)^2=0')ezplot('(x-2)^2+3*(y-1)^2-3');holdonezplot('y-6*(x-2)^2=0');gridonaxis([04-14]);作业31、已知系统的传递函数为，试在MATLAB中建立其传递函数模型，并将其转化为零极增益模型和部分分式模型。解%[n,d]=numden(G);num=sym2poly(n);den=sym2poly(d);symss;G=4*(s+2)*(s^2+6*s+6)^2/(s*(s+1)^3*(s^3+3*s^2+2*s+5));[n,d]=numden(G);num=sym2poly(n);den=sym2poly(d);[z,p,k]=tf2zp(num,den)[r1,p1,k1]=residue(num,den)2、一反馈控制系统的前向通道传递函数和反馈通道传递函数分别为和，试在MATLAB中实现两个子系统的负反馈联结，给出该控制系统的传递函数和以零极点形式表示的数学模型，并绘出单位阶跃响应曲线图。解：G1=tf(3*[16],conv([1,1],[1,3,5]));H=tf([51],[151])Gtf=feedback(G1,H)Gzpk=zpk(Gtf)Gzpk=zpk(Gtf)step(Gtf);holdon;step(Gzpk,'r');3、二阶系统的传递函数为，设其固有频率，在阻尼时，分别画出其单位阶跃响应和单位脉冲响应曲线。解：wn=10;symsksforks=[0.10.30.71]num=[1];den=[120*ks100];figure(1);step(num,den);holdon;figure(2);impulse(num,den);holdon;endfigure(1);gridon;figure(2);gridon;4、已知某控制系统的开还传递函数为，，试绘制系统的开环频率特性曲线，即系统的Bode图。解：K=1.5;num=[K];den=conv([10],conv([11],[12]));bode(num,den)5、二阶系统的传递函数为，设其固有频率，在阻尼时，分别画出其Bode图，研究阻尼系数对二阶系统频率响应的影响。解wn=10;forks=[0.10.30.71]num=[1];den=[120*ks100];bode(num,den);holdon;endgridon;6、系统模型如下所示：，判断系统的稳定性，以及系统是否为最小相位系统。解：num=[3164128];den=[11411052814942117112];[z,p,k]=tf2zp(num,den)%求系统的零极点，并显示ii=find(real(p))>0n1=length(ii);jj=find(real(z)>0)n2=length(jj);if(n1>0)disp('thesystemisUnstable')elsedisp('thesystemisstable')end%判断系统是否为最小相位系统if(n2>0)disp('thesystemisanonminimalphaseone')elsedisp('thesyetemisaminimalphaseone')end%绘制零极点图pzmap(p,z)作业41、由周期信号的分解与合成可知方波分解为多次正弦波之和。如图所示的周期方波，其傅立叶级数为用MATLAB演示谐波合成情况。解：t=-2*pi:0.1:2*pi;

sum=0;

for

a=1:2:300

a1=1./a;

sum=sum+4*a1*sin(a*t)/pi;

plot(t,sum)

grid

on

pause(0.25)

end00pi2pi2、考虑简单的线性微分方程，且方程的初值为，试求该方程的解析解和在区间上的数值解，并比较二者得出的曲线。解：syms

x1

y1

x

y;

y1=dsolve('Dy1+2*y1/x-4*x','y1(1)=2','x')

f=inline('4*x-2*y/x','x','y')

[x,y]=ode45(f,[1,2],2)

subplot(121);

ezplot(y1);

axis([1

2

0

4])

subplot(122);

plot(x,y,'r')

axis([1

2

0

4])3、设二阶连续系统，其特性可用常微分方程表示为，求其冲击响应。若输入为，求其零状态响应。解：num=[1];den=[1

2

8];impulse(num,den)symssg1XGtG=1/(s^2+2*s+8);X=laplace(3*t+cos(0.1*t))

;g1=G*X;[n,d]=numden(g1);num=sym2poly(n);den=sym2poly(d);figure(2);impulse(num,den)axis([0

5

0

2]);4、一阶低通电路的频率响应下图所示是一阶低通电路，若以为响应，求频率响应函数画出其幅频响应(幅频特性)和相频的响应(相频特性)。解：G(jw)=1/1+jwrcden=[2

1];num=[1];bode(num,den);gridon;5、级放大器，每级的传递函数均为，求阶跃响应，画出不同时的波形和频率特性。解：w0=input('ÇëÊäÈëw0µÄÖµ:')num=[w0];den=[1

w0];G0=tf(num,den);forn=1:6G=G0^n;

figure(1);

step(G);

gridon;

holdon;

figure(2);

bode(G);

gridon;

holdon;end6、下图试典型的二阶动态电路，其零输入响应有过阻尼，临界阻尼和欠阻尼三种情况。如已知，，，初始值，，求时的和的零输入响应，并画出波形。解:clear,formatcompactL=0.5;

C=0.02;R=12.5;

uc0=1;

iL0=0;dt=0.01;t=0:dt:3;num=[uc0,R/L*uc0+iL0/C];den=[1,R/L,1/L/C];[r,p,k]=residue(num,den);ucn=r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t);iLn=C*diff(ucn)/dt;figure(1),plot(t,ucn),holdon;gridon;figure(2),plot(t(2:end),iLn),holdon;gridon;

作业51、二阶系统的传递函数为，设其固有频率，在阻尼，输入信号为单位阶跃信号时，分别按以下要求利用Simulink绘制系统的响应曲线：