matlab考试题及答案.doc

《matlab考试题及答案.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、%1、编制一个解数论问题的函数文件：取任意整数，若是偶数，则用2除，否则乘3加1，重复此过程，直到整数变为1。functionf=NO_1(X);Y(1)=X;k=1;while(X~=1)k=k+1;if(mod(X,2)==0)X=X/2;elseX=3*X+1;endY(k)=X;endplot(Y,'b.')end%2、编制程序产生一个数组，满足：a1=1，a2=1，从第三个元素开始，每个元素等于前两个元素的和，直到数组的前后两个元素的比值比0.618小于1e-4，并且以红色点线的形式画出这个数组。clear;A(1)=1;A(2)=1;i=

2、3;Z=1;while(abs(Z)>=1e-4)A(i)=A(i-1)+A(i-2);Z=A(i-1)/A(i)-0.618;i=i+1;endplot(A,'r.')%3、编写一个函数，能够产生分段函数。functiony=test_3_1(X)if(X<=2)y=0.5*X;elseif(X>6)y=0.5;elsey=1.5-0.25*X;endend%调用分段函数，绘制曲线。clear;i=1;forj=0:0.01:2;x(i)=j;y(i)=test_3_1(j)*test_3_1(j+2);i=i+1;endplot(x,y)%4、在

3、2pi周期内画正弦函数曲线，并加注坐标轴标识和标题，然后在3pi/4，pi，5pi/4处分别加入带箭头的说明性文本，最后加注图例。clear;t=0:pi/50:2*pi;n=length(t);y=sin(t);plot(t,y,'-bo','linewidth',1)xlabel('X');ylabel('Y');title('正弦函数曲线');text(2.356,0.707,'leftarrow3pi/4','FontSize',18);text(3.142,0.000,'leftarrow3pi/4','FontSize',18);te

4、xt(3.927,-0.707,'leftarrow3pi/4','FontSize',18);hleg1=legend('sin(x)');%5、A为任意一个n*m矩阵，写程序来计算A中有多少个零元素，并输出个数。A=input('输入一个矩阵A=')n=length(find(A==0))%6、A为任意一个向量，写程序找出A中的最小元素，并且输出这个最小元素。A=input('输入一个向量A=');x=length(A);i=1;y=A(i);while(iA(i))y=A(i);endenddisp('最小元素是：

5、')disp(y)%7、某系统由四个典型环节构成，试根据四阶龙格—库塔公式，求输出量y的阶跃响应，取仿真时间Tf=10s，%步长h=0.02s。r=1;P=[0.110.51;0110;2120;101100];W=[000-1;1000;0100;0010];W0=[1;0;0;0];Wc=[0001];Tf=10;h=0.02;A1=diag(P(:,1));B1=diag(P(:,2));C1=diag(P(:,3));D1=diag(P(:,4));H=B1-D1*W;Q=C1*W-A1;A=inv(H)*Q;B=inv(H)*C1*W0;x

6、=[zeros(length(A),1)];y=[zeros(length(Wc(:,1)),1)];t=0;fori=1:Tf/hK1=A*x+B*r;K2=A*(x+h*K1/2)+B*r;K3=A*(x+h*K2/2)+B*r;K4=A*(x+h*K3)+B*r;x=x+h*(K1+2*K2+2*K3+K4)/6;y=[y,Wc*x];t=[t,t(i)+h];endplot(t,y)%8、系统模型1）判断系统的稳定性2）判定系统是否为最小相位系统。num=[3164128];den=[11411052814942117112];[z,p,k]

7、=tf2zp(num,den);pp=find(real(p)>0);n1=length(pp);if(n1>0)disp('系统不稳定⊙﹏⊙

8、

9、

10、不稳定极点是：')disp(p(ii));elsedisp('我靠，系统竟然稳定!0.0.');endzz=find(real(z)>0);n2=length(zz);if(n2>0)disp('系统是最小相位系统~-。-！')elsedisp('系统为非最小相位系统@_@')endpzmap(p,z)%9、对于典型二阶系统，请编写程序绘制时系统的单位阶跃响应。clear;w=[2:2:10];zeta=

11、0.707;figure(1);holdonforWn=wnum=Wn.^2;den=[1,2*zeta*W