2、用MATLAB求解下面的的方程组。
21?2??x1??4??7?9153?2??x??7?????2??????2?2115??x3???1???????13213??x4??0? ? >> B=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13]; >> C=[4;7;-1;0]; >> r=B\\C
?x?y?z?1?x?2y?z?w?8???2x?y?3w?3??3x?3y?5z?6w?5
>> D=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6]; 2、用MATLAB求解下面的的方程组。
21?2??x1??4??7?9153?2??x??7????2???? (1)???2?2115??x3???1???????13213???x4??0?
>> B=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13]; >> C=[4;7;-1;0]; >> r=B\\C
?x?y?z?1?x?2y?z?w?8?(2)? ?2x?y?3w?3??3x?3y?5z?6w?5
>> D=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6]; >> E=[1;8;3;5]; >> S=D\\E
3、已知多项式a(x)=x2+2x+3,b(x)=4x2+5x+6,求a,b的积。
>> A=[1 2 3]; >> B=[4 5 6];
>> a=poly2str(A,'x'); >> b=poly2str(B,'x'); >> C=conv(A,B); >> c=poly2str(C,'x'); >> c=poly2str(C,'x')
d2ydydy 4、求解微分方程2?2?2y?0,y(0)?1,(0)?0,并画出相平面图。
dtdtdt.m文件
function y=fn(t,y);
y=[y(2);-2*y(1)-2*y(2)];
>> tspan=[0,30]; >> y0=[1;0];
>> [t,y]=ode45(@fn,tspan,y0); >> plot( y(:,1),y(:,2)) >> plot( y(:,1),y(:,2)); >> plot( y(:,1),y(:,2)) >> tspan=[0,30]; y0=[1;0];
[t,y]=ode45(@fn,tspan,y0); plot( y(:,1),y(:,2))
5、已知系统冲激响应为h(n)=[0.05,0.24,0.40,0.24,0.15,-0.1,0.1] ,系统输入u(n)由指令randn('state',1);u=2*(randn(1,100)>0.5)-1产生,该输入信号
的起始作用时刻为0。试用直杆图(提示:用stem指令)画出分别显示该系统输入、输出信号的两张子图。
>> h=[0.05,0.24,0.40,0.24,0.15,-0.1,0.1];
randn('state',1); u=2*(randn(1,100)>0.5)-1; c=conv(h,u); k=[0:105]; stem(k,c)
实验七 采用SIMULINK的系统仿真
一、实验目的:
1、熟悉SIMULINK 工作环境及特点
2、掌握线性系统仿真常用基本模块的用法 3、掌握SIMULINK 的建模与仿真方法 二、实验内容 1、练习
SIMULINK仿真实际应用.
d2x2dx??(1?x)?x?02dt利用SIMULINK求解微分方程 dt,方程的初始条件为x(0)=1,
dx(0)?0)dt。在增益模块‘Gain’取值分别为2和100的情况下(即??2,??100运行,
给出运行结果。
2、已经系统框图,求传递函数,并求阶跃响应。
UA1+-1s?1+-A2G21s?1A3+-G11s?1G31s?1G4Y1s?1G6-+A51s?1G5-+A4