下载文档原格式
第四章1.a=input('请输入一个4位数:');while (a<1000|a>9999)a=input('输入错误，请重新输入一个4位数:'); endb=fix(a/1000);c=rem(fix(a/100),10);d=rem(fix(a/10),10);e=rem(a,10);b=b+7;c=c+7;d=d+7;e=e+7;b=rem(b,10);c=rem(c,10);d=rem(c,10);e=rem(e,10);g=b;b=d;d=g;g=c;c=e;e=g;a=1000*b+100*c+10*d+e;disp(['加密后：',num2str(a)])2.逻辑表达式法：a=input('请输入a: ');b=input('请输入b: ');c=input('请输入c: ');x=0.5:1:5.5;x1=(x>=0.5&x<1.5);x2=(x>=1.5&x<3.5);x3=(x>=3.5&x<=5.5);y1=a.*(x.^2)+b.*x+c;y2=a*(sin(b)^c)+x;y3=log(abs(b+c./x));y=y1.*x1+y1.*x2+y3.*x3; disp(y)if语句法：a=input('请输入a: ');b=input('请输入b: ');c=input('请输入c: ');for x=0.5:1:5.5if x>=0.5 & x<1.5y=a.*(x.^2)+b.*x+celseif x>=1.5 & x<3.5y=a*(sin(b)^c)+xelseif x>=3.5 & x<5.5y=log(abs(b+c./x))endendswitch语句法：a=input('请输入a: ');b=input('请输入b: ');c=input('请输入c: ');for x=0.5:1:5.5switch floor(x/0.5)case {1,2}y=a.*(x.^2)+b.*x+c;case {3,4,5,6}y=a*(sin(b)^c)+x;case {7,8,9,10}y=log(abs(b+c./x)); enddisp(y)end3.x=fix(rand(1,20)*89)+10;x1=mean(x);n=find(rem(x,2)==0 & x<x1);disp(['小于平均数的偶数是：',num2str(x(n))]);4.(1)A=input('请输入20个数的一个行向量：');a=A(1);b=A(1);for m=Aif a>=ma=m;elseif b<=mb=m;endenddisp(['最小数是：',num2str(a)])disp(['最大数是：',num2str(b)])(2)A=input('请输入20个数的一个行向量：'); maxval=max(A)minval=min(A)5.s=0;for a=0:63c=2^a;s=s+c;enddisp(['2的0次方到63次方的和是：',num2str(s)])k=0:63n=2.^ks=sum(n)6.（1）sum1=0;for n=1:100x=(-1)^(n+1)*(1/n);sum1=sum1+x;enddisp(['当n取100时： sum=',num2str(sum1)])sum2=0;for n=1:1000x=(-1)^(n+1)*(1/n);sum2=sum2+x;enddisp(['当n取1000时： sum=',num2str(sum2)])sum3=0;for n=1:10000x=(-1)^(n+1)*(1/n);sum3=sum3+x;enddisp(['当n取10000时：sum=',num2str(sum3)])（2）sum1=0;n1=0;for n=1:2:100x=(-1)^n1*(1/n);sum1=sum1+x;n1=n1+1;enddisp(['当n取100时： sum=',num2str(sum1)])sum2=0;n2=0;for n=1:2:1000x=(-1)^n2*(1/n);sum2=sum2+x;n2=n2+1;enddisp(['当n取1000时： sum=',num2str(sum2)])sum3=0;n3=0;for n=1:2:10000x=(-1)^n3*(1/n);sum3=sum3+x;n3=n3+1;enddisp(['当n取10000时：sum=',num2str(sum3)])（3）sum1=0;for n=1:100x=1/(4^n);sum1=sum1+x;enddisp(['当n取100时： sum=',num2str(sum1)])sum2=0;for n=1:1000x=1/(4^n);sum2=sum2+x;enddisp(['当n取1000时： sum=',num2str(sum2)])sum3=0;for n=1:10000x=1/(4^n);sum3=sum3+x;enddisp(['当n取10000时：sum=',num2str(sum3)])（4）sum1=1;for n=1:100x=4*n*n/(2*n-1)/(2*n+1);sum1=sum1*x;enddisp(['当n取100时： sum=',num2str(sum1)])sum2=1;for n=1:1000x=4*n*n/(2*n-1)/(2*n+1);sum2=sum2*x;enddisp(['当n取1000时： sum=',num2str(sum2)])sum3=1;for n=1:10000x=4*n*n/(2*n-1)/(2*n+1);sum3=sum3*x;enddisp(['当n取10000时：sum=',num2str(sum3)])7.函数文件function f=fibnacci(n)if n==1 | n==2f=1;elsef=fibnacci(n-1)+fibnacci(n-2);end命令文件：shulie=[];for k=1:nshulie=[shulie fibnacci(k)];endshulie8.function [f1,f2]=juzhenji(x1,x2)f1=x1*x2;f2=x1.*x2;命令文件：clear alla=input('请输入一个矩阵：');b=input('请再输入一个矩阵：(注意：两矩阵要可以相乘)'); [f1,f2]=juzhenji(a,b);disp(f1)disp(f2)9.function sum=qiuhe(n,m)if n<=1sum=0;elsesum=n^m+qiuhe(n-1,m);end命令文件：clear ally=qiuhe(100,1)+qiuhe(50,2)+qiuhe(10,-1);disp(y)10.s=0；a=[12,13,14;15,16,17;18,19,20;21,22,23];for k=afor j=1:4if rem(k(j),2)~=0 s=s+k(j);endendendss =108（2）global xx=1:2:5;y=2:2:6;sub(y);xyfunction fun=sub(z) global xz=3*x;x=x+z;x =4 12 20 y =2 4 6。
第4章数值运算习题 4 及解答1 根据题给的模拟实际测量数据的一组t和)(t y试用数值差分diff或数值梯度gradient指令计算)(t y'，然后把)(t y和)(t y'曲线绘制在同一张图上，观察数值求导的后果。
（模拟数据从prob_data401.mat获得）〖目的〗●强调：要非常慎用数值导数计算。
●练习mat数据文件中数据的获取。
●实验数据求导的后果●把两条曲线绘制在同一图上的一种方法。
〖解答〗（1）从数据文件获得数据的指令假如prob_data401.mat文件在当前目录或搜索路径上clearload prob_data401.mat（2）用diff求导的指令dt=t(2)-t(1);yc=diff(y)/dt; %注意yc的长度将比y短1plot(t,y,'b',t(2:end),yc,'r')（3）用gradent 求导的指令（图形与上相似）dt=t(2)-t(1);yc=gradient(y)/dt;plot(t,y,'b',t,yc,'r')grid on〖说明〗● 不到万不得已，不要进行数值求导。
● 假若一定要计算数值导数，自变量增量dt 要取得比原有数据相对误差高1、2个量级以上。
● 求导会使数据中原有的噪声放大。
2 采用数值计算方法，画出dt tt x y x ⎰=0sin )(在]10 ,0[区间曲线，并计算)5.4(y 。
〖提示〗● 指定区间内的积分函数可用cumtrapz 指令给出。
● )5.4(y 在计算要求不太高的地方可用find 指令算得。
〖目的〗● 指定区间内的积分函数的数值计算法和cumtrapz 指令。
● find 指令的应用。
〖解答〗dt=1e-4;t=0:dt:10;t=t+(t==0)*eps;f=sin(t)./t;s=cumtrapz(f)*dt;plot(t,s,'LineWidth',3)ii=find(t==4.5);s45=s(ii)s45 =1.65413 求函数x ex f 3sin )(=的数值积分⎰=π0 )(dx x f s ，并请采用符号计算尝试复算。
第二章1、x=[2,4];y=x八3+(x-0.98)八2./(x+1.35)八3-5*(x+1./x) 2、y=cos(pi/3)-(9-sqrt(2))A(1/3) 3、a=3;A=4; b=a.A2; B=b.A2-1; c=a+A-2*B;C=a+2*B+c4、x=[1 2 3;4 5 6;7 8 9]Desktop->Workspace, 双击变量x 5、clearx=magic(3) y=randn(3,3) xy=[x,y] yx=[x,y] z=xy(:,1:2) 6、clear x=eye(4,4)y=triu(x)7、clear x=rand(4,5) y=x([1,2],:)z=(y>=0.3).*y8、clear x=randn(5,5) y=inv(x)9、clear x=randn(5,5)z=x A510、clearA=[1,4,8;-3,6,-5;2,-7,-12];B=[5,4,3;6,-2,3;-1,3,-9];C=A*BD=A.*B11、clear x=linspace(0,2*pi,125);y=cos(x).*(0.5+3*sin(x)./(1+x.A2)); plot(x,y)12、z=-45:1:45;x=z.*sin(3*z);y=z.*cos(3*z);plot3(x,y,z);13、x=-2:0.1:2,y=x;[x,y]=meshgrid(x,y);z=x.A2.*exp(-x.A2-y.A2)surf(x,y,z);14、x=-2:0.1:2,y=x;[x,y]=meshgrid(x,y);z=x.A2.*exp(-x.A2-y.A2);surf(x,y,z);z1=0.05*x-0.05*y+0.1;hold on,mesh(x,y,z1);15、(1)n=2;alfa=0;t=0:0.1:10;x=cos(t);y=sin(n*t+alfa); subplot(2,2,1);plot(t,x,t,y);(2)n=2;alfa=0;t=0:0.1:10;x=cos(t);y=sin(n*t+alfa); subplot(2,2,1);plot(t,x,t,y);n=2;alfa=pi/3;t=0:0.1:10;x=cos(t);y=sin(n*t+alfa); subplot(2,2,2);plot(t,x,t,y);n=2;alfa=pi/2;t=0:0.1:10;x=cos(t);y=sin(n*t+alfa);subplot(2,2,3);plot(t,x,t,y);n=2;alfa=pi;t=0:0.1:10;x=cos(t);y=sin(n*t+alfa);subplot(2,2,4);plot(t,x,t,y);球体clearfor i=1:100for j=1:100x(i,j)=i/100*cos(j*2*pi/100);y(i,j)=i/100*sin(j*2*pi/100);z(i,j)=sqrt(1.001-x(i,jF2-y(i,j)八2);endendsurf(x,y,z);hold on;surf(x,y,-z);axisequal棱锥clearfor i=1:100for j=1:100x(i,j)=i/100*cos(j*2*pi/4);y(i,j)=i/100*sin(j*2*pi/4);z(i,j)=sqrt(x(i,j)A2+y(i,j)A2);endendsurf(x,y,-z);求最大值x=[67 87 56 99 43] max=0;for i=1:5 if max>=x(i) max=max; elsemax=x(i) endendmax求最小值x=[67 87 56 99 43] min=inf;for i=1:5if min<=x(i) min=min;else min=x(i) end end min求和x=[67 87 56 99 43]; sum=0; for i=1:5 sum=sum+x(i) end sum第三章1、h0=[446,714,950,1422,1634];t0=[7.04,4.28,3.40,2.54,2.13];t1=interp1(h0,t0,500,'linear') 2、x0=[1,0,-1];y0=[0,1,0]; p=polyfit(x0,y0,3); x=-1:0.1:1; y=polyval(p,x);plot(x,y,-x,-y),axis equal3、clear x0=[1.0,1.1,1.2,1.3,1.4];y0=[0.25,0.2268,0.2066,0.1890,0.1736];p=polyfit(x0,y0,3);p1=polyder(p); y=polyval(p1,[1.0,1.2])4、p=[3 4 7 2 9 12]; r=roots(p)5、r=[-3 -5 -8 -9]; p=poly(r)6、function ypie=fun1(x,y)ypie=x42./y-x.*cos(y);[x,y]=ode45( 'fun1' ,[0,5],1)plot(x,y)7、function y=fun2(x)y=x.A4-3*x A3+5*cos(x)+8;fplot( 'fun2' ,[1,5])hold on;x=0:5;y=0*x;plot(x,y)a=fzero(fun2,2)b=fzero(fun2,3)8、a=[2 4 9;4 2 4;9 4 18];[v,gama]=eig(a)9、function xpie=fun3(t,x)f=exp(-t)xpie=[0 1 0 0;1 0 0 -1;0 0 0 1;0 -1 -10]*x+[0;1;0;0]*f[t,x]=ode45( 'fun3' ,[0,5],[0 0 0 0]'); plot(t,x(:,1),'b' ,t,x(:,3), 'r' )10、function f=fun4(x,y)f=4*(x-y)-x八2-y八2;v=-2:0.2:2;[x,y]=meshgrid(v)z=fun4(x,y);surf(x,y,z)求极小值 (有误) function f=fun5(x,y)x=v(1);y=v(2)f=-(4*(x-y)-x.A2-y.A2);v=-50:2:50;[x,y]=meshgrid(v)z=fun4(x,y);surf(x,y,z)min=fminsearch('fun4',[0,0])11、function f=fun6(x,y)f=sin(y)+exp(x)-x.*y.A2ezplot( 'fun6' )12、function f=fun7(x,y)f=(x-y).A2.*(sin(x+y)).A2s=dblquad( 'fun7' ,pi,2*pi,0,pi)14、a=randn(4,4)[l,u]=lu(a)a=randn(4,4)[u,gama,v]=svd(a)15、b=[3 0.5 4];a=[3 4 5 6]; [r,p,k]=residue(b,a)b=[3 0.5 4];a=[3 4 5 6];[r,p,k]=residue(b,a)t=0:0.1:100;y=0*t;for i=1:3 y=y+r(i)*exp(p(i)*t) end plot(t,y) b=[3 0.5 4];a=[3 4 5 6]; impulse(b,a) b=[3 0.5 4];a=[3 4 5 6]; step(b,a)第四章1、n=1:11;x=cos(pi/6*n); subplot(2,1,1);stem(x); y=abs(fft(x)); subplot(2,1,2);stem(y);clear n=0:11;x=cos(pi/6*n); subplot(2,1,1);stem(x); y=abs(fft(x));subplot(2,1,2);stem(y);2、(有误)n=0:19;x=5*0.6.八n;subplot(3,1,1);stem(n,x);for i=-20:39xl(i+21)=x(mod(i+40,20)+1); end n1=-20:39; subplot(3,1,2);stem(n1,x1);x2=x1(10:29); subplot(3,1,3);stem(n,x2)新方法m=10; e=0:19;c=0.6;k=5;a=k*c.A e;a=a';b=circshift(a,m);L=length(a)-1;n=0:L;subplot(2,1,1);stem(n,a);axis([0,L,min(a),max(a)]);subplot(2,1,2);stem(n,b);axis([0,l,min(a),max(a)])3、x=[0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0 0 0 0 0 0 0 0];h=[1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0]y=conv(x,h);stem(y)4、x=[0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0 0 0 0 0 0 0 0];h=[1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0]y=conv(x,h);subplot(2,1,1);stem(y)xx=fft(x);hh=fft(h);yy=ifft(xx.*hh);subplot(2,1,2);stem(yy)5、b=[2 3];a=[1 0.4 1];[z p k]=tf2zp(b,a) 6、b=[4 15.6 6 2.4 -6.4];a=[3 2.4 6.3 -11.4 6];[z p k]=tf2zp(b,a)7、b=[18];a=[18 3 -4 -1];[r p k]=residue(b,a)8、b=[0.2 0.3 1];a=[1 0.4 1]; freqs(b,a)9、Fs=2000;wc=[100 200]/(Fs/2);N=10;[b,a]=butter(N,wc, 'z' ); %freqz(b,a)dimpulse(b,a,100)10、wp=100;ws=200;rp=2;rs=15;Fs=500;wp=wp/(Fs/2);ws=ws/(Fs/2);[N,wc]=buttord(wp,ws,rp,rs); [b,a]=butter(N,wc, 'z' ) freqz(b,a)b=fir1(48,[0.35 0.65]);a=1; freqz(b,a)12、b=fir1(37,0.3);a=1; freqz(b,a)13、低通：F=[0:1/56:1];A=[ones(1,50),zeros(1,56-49)]; b=fir2(56,F,A);a=1; freqz(b,a)11、带通：F=[0:1/56:1];A=[ones(1,25),zeros(1,25),ones(1,7)];b=fir2(56,F,A);a=1;freqz(b,a)高通F=[0:1/56:1];A=[zeros(1,50),ones(1,56-49)]; b=fir2(56,F,A);a=1; freqz(b,a)%滤波器设计(椭圆滤波器)wp1=300;wp2=700;ws1=301;ws2=699; rp=0.1;rs=50;Fs=2000; wp=[wp1]/(Fs/2);ws=[ws1]/(Fs/2);[N,wc]=ellipord(wp,ws,rp,rs,'z')[b,a]=ellip(N,rp,rs,wc,'high')%滤波器特性分析freqz(b,a)%信号采样及时频变换Fs=2000;t=0:1/Fs:10;yt=sin(2*pi*200*t)+sin(2*pi*500*t)+sin(2*pi*800*t); figure(2);subplot(3,1,1);plot(t,yt)figure(2);subplot(3,1,2);plot(t,yt);axis([0,0.1,-5,5]) yf=abs(fft(yt));f=(1:length(yf))/length(yf)*Fs;figure(2);subplot(3,1,3);plot(f,yf)%信号滤波及时频分析yt1=filter(b,a,yt);figure(3);subplot(3,1,1);plot(t,yt1)figure(3);subplot(3,1,2);plot(t,yt1);axis([0,0.1,-5,5]) yf1=abs(fft(yt1));f=(1:length(yf1))/length(yf1)*Fs; figure(3);subplot(3,1,3);plot(f,yf1)。
习题二1.如何理解“矩阵是MATLAB最基本的数据对象”？答：因为向量可以看成是仅有一行或一列的矩阵，单个数据（标量）可以看成是仅含一个元素的矩阵，故向量和单个数据都可以作为矩阵的特例来处理。
因此，矩阵是MATLAB最基本、最重要的数据对象。
2.设A和B是两个同维同大小的矩阵，问：(1)A*B和A.*B的值是否相等？答：(2)(3)答：(4)答：即逆左乘3.(1)答：5])(2)答：(3)答：(4)答：(5)答：(6)将含有12个元素的向量x转换成34⨯矩阵。
答：reshape(x,3,4);(7)求一个字符串的ASCII码。
答：abs(‘123’); 或double(‘123’);(8)求一个ASCII码所对应的字符。
答：char(49);4.下列命令执行后，L1、L2、L3、L4的值分别是多少？A=1:9;B=10-A;...L1=A==B;L2=A<=5;L3=A>3&A<7;L4=find(A>3&A<7);答：L1的值为[0, 0, 0, 0, 1, 0, 0, 0, 0]L2的值为[1, 1, 1, 1, 1, 0, 0, 0,0]L3的值为[0, 0, 0, 1, 1, 1, 0, 0,0]L4的值为[4, 5, 6]6.当A=[34, NaN, Inf, -Inf, -pi, eps, 0]时，分析下列函数的执行结果：all(A)，any(A)，isnan(A)，isinf(A)，isfinite(A)。
答：all(A)的值为0any(A) 的值为1isnan(A) 的值为[ 0, 1, 0, 0, 0, 0, 0]isinf(A) 的值为[ 0, 0, 1, 1, 0, 0, 0]isfinite(A) 的值为[1, 0, 0, 0, 1, 1, 1]7.用结构体矩阵来存储5名学生的基本情况数据，每名学生的数据包括学号、姓名、专业和6门课程的成绩。
《MATLAB 原理及应用》实验报告三.课后练习题答案1.为 ⎪⎩⎪⎨⎧<-=>+=01001x x x xx x y 编写赋值程序。
程序如下：①建立如下的M 文件：x=input('x=');%让用户通过键盘输入数值、字符串或表达式if x>0y=x+1;elseif x==0y=x;else x<0y=x-1;e nd程序执行结果如下>> kh1 %在当前工作目录下，文件名为“kh.1.m ” x=1>> yy =22.使用for ... end循环的array向量编程求出1+3+5...+100 的值程序如下：sum=0;>> for k=1:2:100sum=sum+k;end>> sumsum =25003.计算1+3+5...+100 的值，当和大于1000时终止计算。
程序如下：sum=0;for m=1:2:100; %建立1 3 5….100的向量if sum<=1000 %如果sum小于1000则可以继续加sum=sum+m; %累加elsebreak; %若sum的结果不符合条件就跳出整个循环endend结果为：sum =1024k =653.1计算从1开始多少个自然数之和超过100。
程序如下：>> sum=0;n=0;>> while sum<=100n=n+1;sum=sum+n;end结果为：n =14sum =1054.求1！+2！+3！+……+8！的值程序如下：n=1;sum=1;for m=2:8; %循环7次使得得到各次阶乘n=n*m;sum=sum+n; %累加end结果为:sum =462335.写程序，判断一年是否为闰年，符合下面两条件之一：（1990~2014）A、能被4整除，不能被100整除B、能被400整除程序如下count=0;for y=1990:2014;if((rem(y,4)==0&rem(y,100)~=0)|(rem(y,4)==0&rem(y,400)~=0));count=count+1;endend结果为:count =5。
欢迎共阅习题二1.如何理解“矩阵是MATLAB 最基本的数据对象”？ 答：因为向量可以看成是仅有一行或一列的矩阵，单个数据（标量）可以看成是仅含一个元素的矩阵，故向量和单个数据都可以作为矩阵的特例来处理。
(2) 删除矩阵A 的第7号元素。
答：A(7)=[](3) 将矩阵A 的每个元素值加30。
答：A=A+30;(4) 求矩阵A 的大小和维数。
答：size(A);ndims(A);(5) 将向量 t 的0元素用机器零来代替。
答：t(find(t==0))=eps; (6) 将含有12个元素的向量 x 转换成34⨯矩阵。
答：reshape(x,3,4); (7) 求一个字符串的ASCII码。
答：abs(‘123’); 或E 。
答：B=A(1:3,:); C=A(:,1:2); D=A(2:4,3:4);E=B*C;(2) 分别求E<D 、E&D 、E|D 、~E|~D 和find(A>=10&A<25)。
答：E<D=010001⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦，E&D=110111⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦，11⎡⎤⎢⎥答：student(1).id='0001';student(1).name='Tom'; student(1).major='computer';student(1).grade=[89,78,67,90,86,85]; 8.建立单元矩阵B 并回答有关问题。
B{1,1}=1;B{1,2}='Brenden';B{2,1}=reshape(1:9,3,3); B{2,2}={12,34,2;54,21,3;4 ,23,67};(1)size(B)和ndims(B)的值分别是多少？答：size(B) 的值为2, 2。
ndims(B) 的值为2。
(2)B(2)和B(4)的值分别是(2)建立5×6随机矩阵A，其元素为[100,200]范围内的随机整数。
1. 求下列联立方程的解81025695832475412743-=+-+-=-+-=++-=--+w z y x w z x w z y x w z y x程序A=[3,4,-7,-12;5,-7,4,2;1,0,8,-5;-6,5,-2,10];B=[4;-3;9;-8];X=A\B解：X =-1.4841-0.68160.5337-1.24292. 设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡------=81272956313841A ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----=793183262345B 求C1=AB’;C2=A’B;C3=A.*B,并求它们的逆阵。
程序 A=[1,4,8,13;-3,6,-5,-9;2,-7,-12,-8];B=[5,4,3,-2;6,-2,3,-8;-1,3,-9,7];C1= A*B',C2 = A'*B, C3 = A.*Binv(C1),inv(C2),inv(C3)3. a. 列出2×2阶的单位矩阵I, 4×4阶魔方矩阵M 和4×2阶的全幺矩阵A,全零矩阵Bb. 将这些矩阵拼接为6×6阶的矩阵C:⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=M B A I C ' c. 求出C 的第2,4,6行,组成3×6阶的矩阵C1,及第2,4,6,裂,组成6×3阶的矩阵C2,d. 求D=C1C2及D1=C2C1.程序 >>I=eye(2),A=ones(4,2),B=zeros(4,2),M=magic(4),C=[I,A';B,M]>> C1=C([2,4,6,],:),C2=C(:,[2,4,6,])>> D=C1*C2, D1=C2*C14．设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y把x=0~2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线解：程序x=linspace(0,2*pi,101)y=cos(x).*(0.5+3*sin(x)./(1+x.^2));plot(x,y),grid5．求代数方程3x 5+4x 4+7x 3+2x 2+9x+12=0的所有根。
习题4.1绘制图形)y t--=8t)(tsin(21e≤t并在横轴标注Time纵0≤轴标注“amplitude”，图形的标题为Decaying-oscillating Exponential t=0:0.1:8;y=1-2*exp(-t).*sin(t);plot(t,y,'k')grid ontitle('Decaying-oscillating Exponential')xlabel('Time');ylabel('amplitude')4.2绘制如下函数图形并加上适当的修饰3008.0)309.0cos(5)(22.0≤≤+︒-=--t e t e t y ttt=0:0.1:30;y=5*exp(-0.2*t).*cos(0.9*t-pi/6); plot(t,y,'r:')4.3在同一张图中绘制下列函数曲线100625.0)24083.2cos(23.1)(625.0)(≤≤+︒+==t t t z t yt=0:0.1:10;y=0.625;z=1.23*cos(2.83*t+4*pi/3)+0.625;plot(t,y,'r:')hold onplot(t,z,'b:')text(0,(0.625-1.23*0.5),'这是z(t=0)的位置')text(10,4*pi/3,'这是z(t=10)的位置')4.4对应于250≤≤t 区间内，在同一图中绘制下列函数曲线tt t te t e t y e t y et y ----+︒-===25.1)128554.0cos(02.2)(02.2)(25.1)(3.033.021t=0:0.1:25;y1=1.25*exp(-t);y2=2.02*exp(-0.3*t);y3=2.02*exp(-0.3*t).*cos(0.554*t-128)+1.25*exp(-t);plot(t,y1,'k')hold onplot(t,y2,'b')hold on plot(t,y3,'r')4.5已知椭圆的长、短轴分别为a=4，b=2，用“小红点线”画椭圆 {)cos()sin(t a x t b y ==t=0:pi/100:2*pi;a=4;b=2;x=a.*cos(t);y=b.*sin(t);plot(x,y,'r.')4.6在同一图形框中绘制两个子图，分别显示下列曲线x=0:0.01:10;y1=sin(2*x).*sin(3*x);y2=0.4*x;subplot(1,2,1),plot(x,y1,'r'),title('y1')subplot(1,2,2),plot(x,y2,'b'),title('y2')4.7试将图形分割成3个子图，并分别绘制lgx在100≤区间内对数坐0≤t标、x半对数坐标和y半对数坐标的函数曲线，并加上适当的图形修饰。
（完整版）matlab经典习题及解答第1章 MATLAB 概论1.1 与其他计算机语⾔相⽐较，MATLAB 语⾔突出的特点是什么？MATLAB 具有功能强⼤、使⽤⽅便、输⼊简捷、库函数丰富、开放性强等特点。
1.2 MATLAB 系统由那些部分组成？MATLAB 系统主要由开发环境、MATLAB 数学函数库、MATLAB 语⾔、图形功能和应⽤程序接⼝五个部分组成。
1.4 MATLAB 操作桌⾯有⼏个窗⼝？如何使某个窗⼝脱离桌⾯成为独⽴窗⼝？⼜如何将脱离出去的窗⼝重新放置到桌⾯上？在MATLAB 操作桌⾯上有五个窗⼝，在每个窗⼝的右上⾓有两个⼩按钮，⼀个是关闭窗⼝的Close 按钮，⼀个是可以使窗⼝成为独⽴窗⼝的Undock 按钮，点击Undock 按钮就可以使该窗⼝脱离桌⾯成为独⽴窗⼝，在独⽴窗⼝的view 菜单中选择Dock ……菜单项就可以将独⽴的窗⼝重新防⽌的桌⾯上。
1.5 如何启动M ⽂件编辑/调试器？在操作桌⾯上选择“建⽴新⽂件”或“打开⽂件”操作时，M ⽂件编辑/调试器将被启动。
在命令窗⼝中键⼊edit 命令时也可以启动M ⽂件编辑/调试器。
1.6 存储在⼯作空间中的数组能编辑吗？如何操作？存储在⼯作空间的数组可以通过数组编辑器进⾏编辑：在⼯作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输⼊修改内容即可。
1.7 命令历史窗⼝除了可以观察前⾯键⼊的命令外，还有什么⽤途？命令历史窗⼝除了⽤于查询以前键⼊的命令外，还可以直接执⾏命令历史窗⼝中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M ⽂件中。
1.8 如何设置当前⽬录和搜索路径，在当前⽬录上的⽂件和在搜索路径上的⽂件有什么区别？当前⽬录可以在当前⽬录浏览器窗⼝左上⽅的输⼊栏中设置，搜索路径可以通过选择操作桌⾯的file 菜单中的Set Path 菜单项来完成。
在没有特别说明的情况下，只有当前⽬录和搜索路径上的函数和⽂件能够被MATLAB 运⾏和调⽤，如果在当前⽬录上有与搜索路径上相同⽂件名的⽂件时则优先执⾏当前⽬录上的⽂件，如果没有特别说明，数据⽂件将存储在当前⽬录上。
matlab课后习题及答案详解第1章MATLAB概论与其余计算机语言对比较，MATLAB语言突出的特色是什么？MATLAB拥有功能强盛、使用方便、输入简捷、库函数丰富、开放性强等特色。
MATLAB系统由那些部分构成？MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。
安装MATLAB时，在选择组件窗口中哪些部分一定勾选，没有勾选的部分此后怎样补安装？在安装MATLAB时，安装内容由选择组件窗口中个复选框能否被勾选来决定，能够依据自己的需要选择安装内容，但基本平台（即MATLAB选项）一定安装。
第一次安装没有选择的内容在补安装时只要依据安装的过程进行，不过在选择组件时只勾选要补装的组件或工具箱即可。
MATLAB操作桌面有几个窗口？怎样使某个窗口离开桌面成为独立窗口？又怎样将离开出去的窗口从头搁置到桌面上？在MATLAB操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是封闭窗口的Close按钮，一个是能够使窗口成为独立窗口的Undock 按钮，点击Undock按钮就能够使该窗口离开桌面成为独立窗口，在独立窗口的view菜单中选择Dock,,菜单项就能够将独立的窗口从头防备的桌面上。
怎样启动M文件编写/调试器？在操作桌面上选择“成立新文件”或“翻开文件”操作时，M文件编写/调试器将被启动。
在命令窗口中键入edit命令时也能够启动M文件编写/调试器。
储存在工作空间中的数组能编写吗？怎样操作？储存在工作空间的数组能够经过数组编写器进行编写：在工作空间阅读器中双击要编写的数组名翻开数组编写器，再选中要改正的数据单元，输入改正内容即可。
命令历史窗口除了能够察看前方键入的命令外，还有什么用途？命令历史窗口除了用于查问从前键入的命令外，还能够直接履行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。
怎样设置目前目录和搜寻路径，在目前目录上的文件和在搜寻路径上的文件有什么差别？目前目录能够在目前目录阅读器窗口左上方的输入栏中设置，搜寻路径能够经过选择操作桌面的file菜单中的SetPath菜单项来达成。
第4章数值运算习题 4 及解答1 根据题给的模拟实际测量数据的一组t和)(t y试用数值差分diff或数值梯度gradient指令计算)(t y'，然后把)(t y和)(t y'曲线绘制在同一张图上，观察数值求导的后果。
（模拟数据从prob_data401.mat 获得）〖目的〗●强调：要非常慎用数值导数计算。
●练习mat数据文件中数据的获取。
●实验数据求导的后果●把两条曲线绘制在同一图上的一种方法。
〖解答〗（1）从数据文件获得数据的指令假如prob_data401.mat文件在当前目录或搜索路径上clearload prob_data401.mat（2）用diff求导的指令dt=t(2)-t(1);yc=diff(y)/dt; %注意yc的长度将比y短1plot(t,y,'b',t(2:end),yc,'r')grid on（3）用gradent 求导的指令（图形与上相似）dt=t(2)-t(1);yc=gradient(y)/dt;plot(t,y,'b',t,yc,'r')grid on〖说明〗● 不到万不得已，不要进行数值求导。
● 假若一定要计算数值导数，自变量增量dt 要取得比原有数据相对误差高1、2个量级以上。
● 求导会使数据中原有的噪声放大。
2 采用数值计算方法，画出dt tt x y x ⎰=0sin )(在]10 ,0[区间曲线，并计算)5.4(y 。
〖提示〗● 指定区间内的积分函数可用cumtrapz 指令给出。
● )5.4(y 在计算要求不太高的地方可用find 指令算得。
〖目的〗● 指定区间内的积分函数的数值计算法和cumtrapz 指令。
● find 指令的应用。
〖解答〗dt=1e-4;t=0:dt:10;t=t+(t==0)*eps;f=sin(t)./t;s=cumtrapz(f)*dt;plot(t,s,'LineWidth',3)ii=find(t==4.5);s45=s(ii)s45 =1.65413 求函数x ex f 3sin )(=的数值积分⎰=π0 )(dx x f s ，并请采用符号计算尝试复算。
clc,clear,clf,close all;disp('题解4.5,4.6,4.7');pause;echo on;a=[3 4 7 2 9 12],roots(a),pause,b=[1 0 0 0 0 -1],roots(b),pause,c=[-3 -5 -8 -9],poly(c),echo off;disp('题解4.8');pause;echo on;a=[1 2 5 4 3],b=[1],[r,p,k]=residue(b,a)pause,t=0:0.2:10;ys=r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t)+r(3)*exp(p(3)*t)+r(4)*exp(p(4)*t); h=plot(t,ys);set(h,'linewidth',2),echo off;disp('题解4.9,4.10');pause;echo on;R1=randn(8,6),mean(R1),std(R1),pause,mean(R1(:)),std(R1(:)),pause,R=ceil(16*rand(4,6)),a=R(:,[1,2,3,4]),inv(a),echo off;disp('题解4.11');echo on;t=[0:0.05:10];x1=2.*cos(t)+3.*t;y1=2.*sin(t)+3;plot(x1,y1),shg,hold on;x2=3.*cos(t)+3.*t;y2=3.*sin(t)+3;plot(x2,y2),shgx3=4.*cos(t)+3.*t;y3=4.*sin(t)+3;plot(x3,y3),shghold off;pause;figure;plot(x1,y1,x2,y2,x3,y3)echo off;disp('题解4.12');echo on;t=[-2*pi:0.05:2*pi];a=0;x=sin(t);y1=sin(t+a);y2=sin(2.*t+a);y3=sin(3.*t+a);y4=sin(4.*t+a);subplot(2,2,1),plot(x,y1), subplot(2,2,2),plot(x,y2), subplot(2,2,3),plot(x,y3), subplot(2,2,4),plot(x,y4), pause;y11=sin(2.*t);y22=sin(2.*t+pi/3);y33=sin(2.*t+pi/2);y44=sin(2.*t+pi);figure;subplot(2,2,1),plot(x,y11), subplot(2,2,2),plot(x,y22), subplot(2,2,3),plot(x,y33), subplot(2,2,4),plot(x,y44), echo off;disp('题解4.13');pause; echo on;p=[1 -4 0 3 -2 6],x=linspace(-2,8,100);y=polyval(p,x);plot(x,y),pauseaxis([-2 8 -20 20])pause,roots(p),echo off;disp('题解4.14')echo on;z=[0:0.05:10];x=z.*sin(3.*z);y=z.*cos(3.*z);plot3(x,y,z);pause;figure,plot3(x,y,z,'r:')pause;xlabel('x'),ylabel('y'),zlabel('z'), echo off;disp('题解4.15 4.16');echo on;[x,y]=meshgrid(-2:0.1:2,-2:0.1:2); z=x.^2.*exp(-(x.^2+y.^2)); mesh(z);hold on;pause;z1=0.05.*x-0.05.*y+0.1;mesh(z1);echo off;for i=1:41for j=1:41if z(i,j)==z1(i,j)a=z(i,j);plot3(i,j,a,'r+');endendendhold off;disp('题解4.18');pause;type f1.mi=1;a(1)=fzero('f1',-2);for m=-2:1:8b=fzero('f1',m);if abs(a(i)-b)>0.00001i=i+1;a(i)=b;endendapause,echo on,fplot('f1',[-2 8]);shg;hold on,plot(a,0,'ro'),shg,pause,axis([-2 8 -20 20]),grid,hold off;echo off;disp('题解4.18第二问');pause; type f11.mi=1;a(1)=fzero('f11',-2);for m=-2:1:8b=fzero('f11',m);if abs(a(i)-b)>0.00001i=i+1;a(i)=b;endendapause,echo on,fplot('f11',[-2 8]);shg;hold on,plot(a,0,'ro'),shg,pause,axis([-2 8 -20 20]),grid,hold off;echo off;disp('题解4.19');echo on;x=[0:0.05:4];y=f31(x);plot(x,y)echo off;disp('题解420');pause;type f420.mi=1;a(1)=fzero('f420',0.1);for m=0.1:0.1:4b=fzero('f420',m);if abs(a(i)-b)>0.0001i=i+1;a(i)=b;endendapause,echo on,fplot('f420',[0 4]);shg;hold on,plot(a,f420(a),'ro'),shg,pause, axis([0 4 -20 20]),grid,hold off;pause,quad('f420',a(1),a(2)),echo off;disp('题解421');pause;type f421.mecho on[x,y]=ode23('f421',[0 5],1);plot(x,y),shg,pause,p=polyfit(x,y,9),xi=0:0.1:5;yi=polyval(p,xi);plot(x,y,'ro',xi,yi,'b')echo off,disp('题解422');pause;echo onpause,s='magic(3)',pause,eval(s);pause,y(3,:)='magic(3)';y(4,:)='magic(4)';y(5,:)='magic(5)';y,pausefor n=3:5eval(y(n,:));endecho offdisp('题解424');pause;echo one=exp(1);pause,st=sprintf('自然对数底数e=%21.20f',e);disp(st) echo offdisp('题解425');pause;echo on='Jone';student.birthplace='Shanghai';student(2).name='David';student(2).birthplace='Nanjing';student(3).name='Tom';student(3).birthplace='Hangzhou';student,pause,name={'Jone';'David';'Tom'};birthplace={'Shanghai';'Nanjing';'Hangzhou'}; studentt={name;birthplace};studentt,pause,student1=['Jone ';'Shanhai'],student2=['David ';'Nanjing'],student3=['Tom ';'Hongzhou'],echo off调用子程序f1.mfunction y=f1(x)y=x.^5-4*x.^4+3*x.^2-2*x+6;f11.mfunction y=f11(x)y=x.^5-4*x.^4+3*x.^2-2*x+6+x.*sin(x);f31.mfunction f1=f31(x)f1=1./((x-2).^2+0.1)+1./((x-3).^4+0.01);f420.mfunction y=f420(x)y=x.^3-2*(x.^2).*sin(x)+5*x.*cos(x)+1./(x); f421.mfunction yp=f421(x,y)yp=(x.^2)./y-x.*cos(y);。
