聊城大學(xué)通信工程專業(yè)matlab編程題答案

1、1. 階乘的函數(shù)(函數(shù)文件)function result=jiecheng(n)result=1;for i=1:n result=result*i;end(腳本文件)n=input('輸入n:');result=jiecheng(n);disp('階乘結(jié)果是:');result2. 編寫一個計算雙曲正弦的函數(shù)。以及調(diào)用這個函數(shù)的腳本m文件。(函數(shù)文件)function sinhxx=-7:0.1:7;y=(exp(x)-exp(-x)/2;plot(x,y);title('plot of (exp(x)-exp(-x)/2');xlabel(

2、'x');ylabel('y');legend('雙曲正弦');grid on;（腳本文件）disp('雙曲正弦函數(shù)');sinhx3. 成績(函數(shù)文件)。function result=paixu(a)result=a;for i=1:length(result); for j=2:length(result); if result(j-1)>result(j); q=result(j-1); result(j-1)=result(j); result(j)=q; end endend(腳本文件)for i=1:10; a

3、(i)=input('輸入a:'); endazuigao=max(a)zuidi=min(a)pingjun=mean(a)result=paixu(a)4.分段函數(shù)的表示。(腳本文件)disp('輸入x,y:');x=input('輸入x:');y=input('輸入y:');if x >= 0 & y >= 0;fun = x + y;elseif x >= 0 & y < 0;fun = x + y2;elseif x < 0 & y >= 0;fun = x2

4、+ y;elsefun = x2 + y2;endfun5.求一元二次方程的根（討論系數(shù)）。（腳本文件）disp('輸入a,b,c');a=input('輸入a：');b=input('輸入b：');c=input('輸入c：');if(b2-4*a*c)<0; a1=(-b)/(2*a); b1=-sqrt(4*a*c-b2)/(2*a); b2=sqrt(4*a*c-b2)/(2*a); x1=complex(a1,b1); x2=complex(a1,b2);endif(b2-4*a*c)=0; x1=(-b)/(2*

5、a); x2=(-b)/(2*a);endif(b2-4*a*c)>0; x1=(-b)/(2*a)-sqrt(b2-4*a*c)/(2*a); x2=(-b)/(2*a)+sqrt(b2-4*a*c)/(2*a);endx1x26. 如何修改程序能使小紅球上下不停運(yùn)動x=0:0.01:20*pi;p=plot3(sin(x),cos(x),x);h=line(0,1,0,'color','red','marker','.','markersize',20,'erasemode','x

6、or');for k=0:0.01:20*pi; set(h,'xdata',sin(k),'ydata',cos(k),'zdata',k); drawnow;endfor k=20*pi:-0.01:0; set(h,'xdata',sin(k),'ydata',cos(k),'zdata',k); drawnow;endclear close all y = dsolve (' D2y+3*Dy+2*y=2*t+2*t2','y(0)=1,Dy(0)=1'

7、;);yht = dsolve (' D2y+3*Dy+2*y=0');yt = dsolve (' D2y+3*Dy+2*y=2*t+2*t2');yp = yt-yht;yh = y-yp;yzi= dsolve (' D2y+3*Dy+2*y=0','y(0)=1,Dy(0)=1');yzs= dsolve (' D2y+3*Dy+2*y=2*t+2*t2','y(0)=1,Dy(0)=1');t=0:0.01:3;figure(1)ezplot(yzi,0,3);hold on ;ezplo

8、t(yzs,0,3);ezplot(y,0,3);axis(0 3 -1 5);hold off;title('全響應(yīng)，零輸入響應(yīng)，零狀態(tài)響應(yīng) ')；figure(2)ezplot(yh,0,3);hold on ;ezplot(yp,0,3);ezplot(y,0,3);axis(0 3 -1 5);hold off;title（'全響應(yīng)，自由響應(yīng)，強(qiáng)迫響應(yīng)'）clear;close all;N=64;n=0:N-1;A=3 -4 2;B=1 2;x1=(1/2).n;y=filter(B,A,x1);x2=n=0;h=filter(B,A,x2);subpl

9、ot(2,1,1);stem(n,y,'fill');grid on;xlabel('n');title('系統(tǒng)響應(yīng)y(n)');subplot(2,1,2);stem(n,h,'fill');grid on;xlabel('n');title('系統(tǒng)單位取樣響應(yīng)h(n)');clear； close all； y = dsolve (' D2y+3*Dy+2*y=2*t+2*t2','y(0)=1,Dy(0)=1');% 自由響應(yīng)和強(qiáng)迫響應(yīng)yht = dsolve

10、(' D2y+3*Dy+2*y=0');%求齊次通解yt = dsolve (' D2y+3*Dy+2*y=2*t+2*t2'); %求非齊次通解yp = yt-yht;%求特解，即強(qiáng)迫響應(yīng)yh = y-yp;%求齊次解，即自由響應(yīng)yzi= dsolve (' D2y+3*Dy+2*y=0','y(0)=1,Dy(0)=1');%求零輸入相應(yīng)yzs= dsolve (' D2y+3*Dy+2*y=2*t+2*t2','y(0)=1,Dy(0)=1');%求零狀態(tài)響應(yīng)t=0:0.01:3;figure

11、(1)ezplot(yzi,0,3);hold on ;ezplot(yzs,0,3);ezplot(y,0,3);axis(0,3；-1，5);hold off;title('全響應(yīng)，零輸入響應(yīng)，零狀態(tài)響應(yīng) ')；figure(2)ezplot(yh,0,3);hold on ;ezplot(yp,0,3);ezplot(y,0,3);axis(0,3；-1，5);hold off;title（'全響應(yīng)，自由響應(yīng)，強(qiáng)迫響應(yīng)'）clear;close all;N=64;n=0:N-1;A=3 -4 2;B=1 2;x1=(1/2).n;y=filter(B,A,x1);x2=n=0;h=filter(B,A,x2);subplot(2,1,1);s