资源预览内容
第1页 / 共4页
第2页 / 共4页
第3页 / 共4页
第4页 / 共4页
亲,该文档总共4页全部预览完了,如果喜欢就下载吧!
资源描述
1.阶乘的函数( (函数文件)函数文件)function result=jiecheng(n)result=1;for i=1:nresult=result*i;end(脚本文件)(脚本文件)n=input(输入n:);result=jiecheng(n);disp(阶乘结果是:);result2. 编写一个计算双曲正弦的函数。以及调用这个函数的脚本m文件。(函数文件(函数文件) )function sinhxx=-7:0。1:7;y=(exp(x)-exp(x)/2;plot(x,y);title(plot of (exp(x)exp(-x))/2);xlabel(x);ylabel(y);legend(双曲正弦);grid on;( (脚本文件脚本文件) )disp(双曲正弦函数);sinhx3。 成绩(函数文件(函数文件) )。function result=paixu(a)result=a;for i=1:length(result);for j=2:length(result);if result(j-1)result(j);q=result(j1);result(j1)=result(j);result(j)=q;endendend( (脚本文件脚本文件) )for i=1:10;a(i)=input(输入a:);endazuigao=max(a)zuidi=min(a)pingjun=mean(a)result=paixu(a)4.4.分段函数的表示。分段函数的表示。( (脚本文件)脚本文件)disp(输入x,y:);x=input(输入x:);y=input(输入y:);if x = 0 y = 0;fun = x + y;elseif x = 0 & y 0;x1=(-b)/(2a)-sqrt(b2-4*ac)/(2a);x2=(-b)/(2a)+sqrt(b2-4*a*c)/(2a);endx1x26。 如何修改程序能使小红球上下不停运动x=0:0。01:20pi;p=plot3(sin(x),cos(x),x);h=line(0,1,0,color,red,marker,.,markersize,20,erasemode,xor);for k=0:0。01:20pi;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 ally = dsolve ( D2y+3Dy+2y=2*t+2*t2,y(0)=1,Dy(0)=1);yht = dsolve ( D2y+3Dy+2*y=0);yt = dsolve ( D2y+3Dy+2y=2t+2*t2);yp = yt-yht;yh = yyp;yzi= dsolve ( D2y+3Dy+2y=0,y(0)=1,Dy(0)=1);yzs= dsolve( D2y+3Dy+2y=2t+2t2, y(0)=1,Dy(0)=1);clear; close all;t=0:0。01:3;figure(1)ezplot(yzi,0,3 ) ;hold on ;ezplot(yzs,0,3 ) ;ezplot(y,0,3);axis(0 3 1 5 );hold off;title(全响应,零输入响应,零状态响应 );figure(2)ezplot(yh,0,3);hold on ;ezplot(yp,0,3) ;ezplot(y, 0,3);axis( 0 3 1 5) ;hold off;title(全响应,自由响应,强迫响应)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) ;stem(n,y,fill) ;grid on;xlabel(n);title(系统响应y(n) );subplot(2,1,2);stem(n,h,fill);grid on;xlabel ( n); title ( 系统单位取样响应h(n);y = dsolve ( D2y+3Dy+2y=2t+2t2,y(0)=1,Dy(0)=1); 自由响应和强迫响应yht = dsolve ( D2y+3*Dy+2*y=0) ;求齐次通解yt = dsolve ( D2y+3*Dy+2*y=2*t+2t2) ; 求非齐次通解yp = ytyht;%求特解,即强迫响应yh = yyp;%求齐次解,即自由响应yzi= dsolve ( D2y+3Dy+2y=0,y(0)=1,Dy(0)=1);求零输入相应yzs= dsolve ( D2y+3Dy+2*y=2t+2t2,y(0)=1,Dy(0)=1);求零状态响应t=0:0。01:3;figure(1)ezplot(yzi,0,3);hold on ;ezplot(yzs,0,3);ezplot(y,0,3);axis(0,3;1,5);hold off;title(全响应,零输入响应,零状态响应 );figure(2)ezplot(yh,0,3);hold on ;ezplot(yp,0,3);ezplot(y,0,3);axis(0,3;1,5);hold off;title(全响应,自由响应,强迫响应)clear;close all;N=64;n=0:N1;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);stem(n,y,fill);grid on;xlabel(n);title(系统响应y(n) );subplot(2,1,2);stem(n,h,fill) ;grid on;xlabel(n);title(系统单位取样响应h(n);
收藏 下载该资源
网站客服QQ:2055934822
金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号