Matlab设计报告Word文档下载推荐.docx
《Matlab设计报告Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《Matlab设计报告Word文档下载推荐.docx(16页珍藏版)》请在冰豆网上搜索。
1-82-1;
0-593;
-70-45];
>
b=[0;
2;
-1;
6];
x=inv(A)*b
结果为:
x=
-0.6386
-0.4210
-0.3529
0.0237
2.分别用solve和fsolve函去数求方程组的解:
solve解法的程序为:
f='
x^2*y^2-2*x-1=0'
;
g='
x^2-y^2-1=0'
[x1,y1]=solve(f,g,'
x'
'
y'
)
结果为:
x1=
1/2+1/2*5^(1/2)
1/2-1/2*5^(1/2)
-1/2+1/2*i*3^(1/2)
-1/2-1/2*i*3^(1/2)
y1=
1/2*(2+2*5^(1/2))^(1/2)
-1/2*(2+2*5^(1/2))^(1/2)
1/2*(2-2*5^(1/2))^(1/2)
-1/2*(2-2*5^(1/2))^(1/2)
1/2*(-6-2*i*3^(1/2))^(1/2)
-1/2*(-6-2*i*3^(1/2))^(1/2)
1/2*(-6+2*i*3^(1/2))^(1/2)
-1/2*(-6+2*i*3^(1/2))^(1/2)
fsolve的解法的程序为:
%ex1.m
functionz=ex1(p)
z
(1)=p
(1)^2*p
(2)^2-2*p
(1)-1;
z
(2)=p
(1)^2-p
(2)^2-1;
x=fsolve('
ex1'
[11])
Optimizationterminated:
first-orderoptimalityislessthanoptions.TolFun.
1.61801.2720
3.用MATLAB计算
。
formatlong;
y=quad('
sin(x)'
0,pi/2)
y=
0.99999999787312
4.绘制[0,4π]区间上的x1=10sint曲线,并要求:
(1)线形为点划线、颜色为红色、数据点标记为加号;
(2)坐标轴控制:
显示范围、刻度线、比例、网络线;
(3)标注控制:
坐标轴名称、标题、相应曲线文本;
>
t0=0:
pi/10:
4*pi;
x1=10*sin(t0);
t1=[0pi/2pi3*pi/22*pi5*pi/23*pi7*pi/24*pi];
x=interp1(t0,x1,t1,'
cubic'
);
plot(t0,x1,'
r-.'
t1,x,'
+'
grid
legend('
10*sin(t0)'
数据点'
axis([0,4*pi,-11,11])
xlabel('
t0'
ylabel('
title('
x1=10*sin(t0)'
x1'
图形为:
5.用for循环和while循环,统计某学生10门功课的平均值。
for循环的程序如下:
a=[30405060706580907585];
s=0;
k=0;
forn=a;
s=s+n;
end
average=s/10
average=
64.50000000000000
While循环的程序如下:
num=0;
i=1;
while(i<
=10)
s=s+a(i);
i=i+1;
num=num+1;
6.已知某MIMO系统模型如下:
试求其零极点增益模型。
num={[11],[10.5];
[32],[12]};
den={[121],[123];
[1021],[152]};
sys=tf(num,den)
Transferfunctionfrominput1tooutput...
s+1
#1:
-------------
s^2+2s+1
3s+2
#2:
s^3+2s+1
Transferfunctionfrominput2tooutput...
s+0.5
s^2+2s+3
s+2
s^2+5s+2
sys1=zpk(sys)
Zero/pole/gainfrominput1tooutput...
(s+1)
-------
(s+1)^2
3(s+0.6667)
-----------------------------------
(s+0.4534)(s^2-0.4534s+2.206)
Zero/pole/gainfrominput2tooutput...
(s+0.5)
---------------
(s^2+2s+3)
(s+2)
--------------------
(s+4.562)(s+0.4384)
7.已知某闭环系统结构框图如下:
图中,
,
,试求其闭环系统仿真模型的传递函数。
num1=[51];
den1=[10];
num2=[25];
den2=[1225];
num3=[0.1];
den3=[1];
[na,da]=series(num1,den1,num2,den2);
[nb,db]=feedback(na,da,num3,den3,-1);
[num,den]=cloop(nb,db);
printsys(num,den)
num/den=
125s+25
----------------------------
s^3+2s^2+162.5s+27.5>
num3=[0.1];
num4=[0];
den4=[0];
[na,da]=series(num1,den1,num2,den2);
[nc,dc]=series(num4,den4,nb,db);
s^3+2s^2+162.5s+27.5
tf(nb,db)
Transferfunction:
--------------------------
s^3+2s^2+37.5s+2.5
8.已知某系统状态空间模型为:
建立simulink仿真模型,并求其单位阶跃响应。
矩阵:
a=[0100;
0010;
-20-13-22-5],
b=[0;
0;
10],c=[1000],d=0。
仿真结果如下图:
9.已知典型二阶系统传递函数如下:
试求
=6,
为0.4、0.7、1.0、1.3时系统的bode图。
1、当
=0.4时,在命令窗口输入:
num1=36;
den1=[14.836];
bode(num1,den1)
得到系统的bode图如下:
2、当
=0.7时,在命令窗口输入:
num2=36;
den2=[18.436];
bode(num2,den2)
3、当
=1.0时,在命令窗口输入:
num3=36;
den3=[11236];
bode(num3,den3)
4、当
=1.3时,在命令窗口输入:
num4=36;
den4=[12*1.3*636];
bode(num4,den4)
10.已知某闭环系统传递函数如下:
(1)
(2)
绘制系统零极点分布图,并判断系统稳定性。
1)z=[j-1-j-1];
k=[8];
p=[00-5-7j-j];
ii=find(real(p)>
0);
n1=length(ii);
if(n1>
0)
disp('
TheUnstablePolesare:
'
dis(p(ii));
elsedisp('
SystemisStable'
pzmap(z,p);
title('
Zero-PoleMap'
SystemisStable
num=100*[11];
den=conv([11],conv([0.21],conv([0.51],[0.81])));
[z,p]=tf2zp(num,den);
ii=find(real(p)>
if(ii)
theunstablepolesare:
disp(p(ii));
systemisstable'
结果:
systemisstable
再在命令窗口输入:
pzmap(num,den)
四、实验总结
本次实验收获颇多,可能底子不好,或者说根本就没底子,做这个报告的题目花的时间比较长,虽然时间花了很多,但我学到了很多这方面的知识。
本次实验内容包括了本学期所学的本课程的主要的内容,做这次实验的题目就需要把整本书的内容都浏览了一遍,所以通过做这些题目对书本的内容又加深了一些了解。
前面的内容都做过上机,有些印象,做起来还比较轻松,但从后面的simulink开始,上课没怎么听懂,接着后面的更听不懂,这样恶性循环,后面的内容几乎什么都不知道,课后虽然看了一会书,但没经过实践操作,所以对后面的内容还是很生疏,还是似懂非懂。
做后面的题目吃力程度是可想而知,但通过自己的研究和与专业里面其他人讨论,最终还是成功的完成了所有的题目,学会了后面没懂的知识,虽然还不是很熟悉,但我还是会通过多编程序来让自己熟悉起来。
所以我觉得任何一门课程都需要自己在底下下功夫才能学好,特别是像这种需要上机操作的课程更要通过实践才能学好,以后我要加强我自己的动手能力,上课专心听讲,课后勤动手,多思考,多问。
五、意见或建议
本课程采用做课程报告来作为考试,这样对于我这种虽然经历无数次考试但还是讨厌考试的人来说很欢喜。
对本课程有如下建议:
第一:
老师很有亲和力,我们大家都很喜欢,但上课的时候只是一味的讲课,跟我们互动比较少,课堂有点枯燥,所以希望能有多点的互动,多调节一下课堂气氛;
第二:
老师其他方面都很好,希望老师继续保持;
最后,本课程是由上机加上课,也就是理论加实践组成,这样很好,让我们学的更透彻,这样的教学方式很好。