matlab复习题及答案1

本文格式为Word版，下载可任意编辑——matlab复习题及答案1=0];

7、已知输入的信号是有限序列x?n??u?n??u?n?10?，而脉冲响应式无限序列：h?n???0.9?u?n?，求y?n??x?n??h?n?。请用matlab编写程序

n将输入序列，输出序列画出。

clearall;closeall;clc;b=[1];a=[1,-0.9];

n=-5:50;x=stepseq(0,-5,50)-stepseq(10,-5,50);x=u(n)-u(n-10);h=0.9.^n.*u(n);

subplot(311);stem(n,x);axis([-10,50,-0.5,2]);title('输入序列');xlabel('n');ylabel('x(n)');subplot(312);stem(n,h);axis([-10,50,-0.5,2]);title('脉冲响应无限序列');xlabel('n');ylabel('h(n)');h=filter(b,a,x);

subplot(313);stem(n,h);axis([-10,50,-0.5,8]);title('输出序列');xlabel('n');ylabel('y(n)');functionf=u(n);f=(n>=0);

function[x,n]=impseq(n0,n1,n2)n=[n1:n2];x=[(n-n0)==0];function[x,n]=stepseq(n0,n1,n2)n=[n1:n2];x=[(n-n0)>=0];

n????xn?0.5u?n?的离散傅里叶变??0,?8、将分为501个等间隔的点，计算

jwXe换

??并画出其模、相角、实部、虚部的曲线。

clearall;closeall;clc;

w=[0:1:500]*pi/500;%将[0,pi]区域分为501个点X=exp(j*w)./(exp(j*w)-0.5*ones(1,501));magX=abs(X);angX=angle(X);realX=real(X);imagX=imag(X);subplot(221);plot(w/pi,magX);grid

title('幅度部分');xlabel('以pi为单位的频率');ylabel('幅度');subplot(223);plot(w/pi,angX);grid

title('相角部分');xlabel('以pi为单位的频率');ylabel('相角');subplot(222);plot(w/pi,realX);grid

title('实部');xlabel('以pi为单位的频率');ylabel('实部');subplot(224);plot(w/pi,imagX);grid

title('虚部');xlabel('以pi为单位的频率');ylabel('虚部');

9、已知其传递函数为H?s??（1）用

4s。

s3?3s2?6s?4S函数编程对一个连续线性时不变系统进行仿真；

（2）请自己建立子系统来实现该系统仿真；

（3）请用Simulink模块库中提供的传递函数模块对该系统进行仿真。

function[sys,x0,str,ts]=ans9a(t,x,u,flag)%á?D??μí3×′ì?·?3ì;%x'=Ax+Bu%y=Cx+Du%?¨ò?A,B,C,D???óA=[010;001;-4-6-3];B=[0;0;1];C=[040];D=0;switchflag,

case0%flag=03?ê??ˉ

[sys,x0,str,ts]=mdlInitializeSizes(A,B,C,D);%?é??A,B,C,D???ó?íè?3?ê??ˉoˉêy

case1%flag=1????á?D??μí3×′ì?·?3ì(μ?êy)sys=mdlDerivatives(t,x,u,A,B,C,D);case3%flag=3????ê?3?sys=mdlOutputs(t,x,u,A,B,C,D);case{2,4,9}%????×÷2?′|àíμ?flagsys=[];%?Tó?μ?flagê±·μ??sys?a?????óotherwise%òì3￡′|àí

error(['Unhandledflag=',num2str(flag)]);end

%?÷oˉêy?áê?

%×óoˉêyêμ??￡¨1￡?3?ê??ˉoˉêyfunction[sys,x0,str,ts]=mdlInitializeSizes(A,B,C,D)%sizes=simsizes;%??è?

sizes.NumContStates=3;%á?D??μí3μ?×′ì?êy?a3

sizes.NumDiscStates=0;%à?é￠?μí3μ?×′ì?êy￡???óú±??μí3′????é2?ó?sizes.NumOutputs=1;%ê?3?D?o?êy??ê?1sizes.NumInputs=1;%ê?è?D?o?êy??ê?1sizes.DirFeedthrough=0;%òò?a???μí32?ê??±í¨μ?sizes.NumSampleTimes=1;%?aà?±?D??a1sys=simsizes(sizes);

str=[];%í¨3￡?a?????ó

x0=[0;0;0];%3?ê?×′ì????óx0￡¨á?×′ì??é??￡?ts=[00];%±íê?á?D?è??ùê±??μ?·???%3?ê??ˉoˉêy?áê?

%×óoˉêyêμ??￡¨2￡??μí3×′ì?·?3ìoˉêyfunctionsys=mdlDerivatives(t,x,u,A,B,C,D)%?μí3×′ì?·?3ìoˉêysys=A*x+B*u;%?aà?D′è??μí3μ?×′ì?·?3ì???óD?ê??′?é%?μí3×′ì?oˉêy?áê?

%×óoˉêyêμ??￡¨3￡??μí3ê?3?·?3ìoˉêy

functionsys=mdlOutputs(t,x,u,A,B,C,D)sys=C*x;%?aà?D′è??μí3μ?ê?3?

10、用S函数实现一个绝对值模块，即输出信号是输入信号求绝对值的

结果，试仿真得出输出结果。

function[sys,x0,str,ts]=ans10(t,x,u,flag)switchflag,case0

sizes=simsizes;%??è?

sizes.NumContStates=0;%á?D??μí3μ?×′ì?êy?a3sizes.NumDiscStates=0;%à?é￠?μí3μ?×′ì?êy￡???óú±??μí3′????é2?ó?

sizes.NumOutputs=1;%ê?3?D?o?êy??ê?1sizes.NumInputs=-1;%ê?è?D?o?êy??ê?1sizes.DirFeedthrough=1;%òò?a???μí32?ê??±í¨μ?sizes.NumSampleTimes=1;%?aà?±?D??a1sys=simsizes(sizes);str=[];%í¨3￡?a?????ó

x0=[];%3?ê?×′ì????óx0￡¨á?×′ì??é??￡?ts=[-10];%±íê?á?D?è??ùê±??μ?·???%3?ê??ˉoˉêy?áê?case3sys=abs(u);case{2,4,9}

sys=[]