MATLAB经典例题.doc

Matlab参考试题主文件main.mclear;while 1fprintf(请输入你要调用的文件的序号.n);fprintf(0:退出程序n );fprintf(1: :绘制n); fprintf(2:一组数据n); fprintf(3:设计一个n) fprintf(4:编写m程序n) a=input( ); while (a4|a0) disp(输入有误,重新输入) a=input( ); end if a=0 break; end switch a case 1 t1; case 2 t2; case 3 t3; case 4 t4; end end1.绘制图形在0 2范围内绘制二维曲线图y=sin(x)*cos(5x)在0 1范围内绘制z = f (x,y) =/sqrt(1-x*x+y*y)+1/sqrt(1+x*x+y*y) 的三维曲面图和三视图function t1()x=0:0.1:2*pi;y=sin(x).*cos(5*x);plot(x,y);grid on;title(y=sin(x)*cos(5x)曲线);xlabel(自变量x);ylabel(函数值y);function t2()x=0:1:20;y=x;X,Y=meshgrid(x,y);Z=1./sqrt(1-X.2+Y2)+1./sqrt(1+X.2+Y.2);subplot(2,2,1);surf(X,Y,Z);title(3D曲面);subplot(2,2,2);plot(x,y);title(xy曲线);subplot(2,2,3);zx=2./sqrt(1+y.2);plot(y,zx);title(yz曲线);subplot(2,2,4);zy=1./sqrt(1-x.2)+1./sqrt(1+x.2);plot(x,zy);title(xz曲线);绘制z=sin(x)*cos(y)的三维网格和三维曲面图，x,y变化范围均为 0 2在0 2范围内绘制以Y轴为对数的二维曲线图。y=|1000sin(4x)|+1function t1()x=0:0.1:2*pi;y=0:0.1:2*pi;X,Y=meshgrid(x,y);Z=sin(X).*cos(Y);figure(1);mesh(X,Y,Z);figure(2);surf(X,Y,Z);function t1()x = 0:0.1:2*pi;y = abs(1000*sin(4*x)+1;semilogy(x,y);set(gca,xtick,0:0.5:6.5);title(y = |1000*sin(4*x)|+1的函数图象);xlabel(x),ylabel(y)legend(y = |1000*sin(4*x)|+1,4)gridv=0,6.5,0,10000;axis(v);在0 2范围内,在同一坐标系下绘制t*t, -t*t, t*t*sin(t)的曲线图function t1()t=0:0.01*pi:2*pi;y1=t.*t;y2=-1*y1;y3=y1.*sin(t);plot(t,y1,r,t,y2,b,t,y3,k);grid on;title(第一小题作图t*t,-t*t,t*t*sin(t)曲线);xlabel(自变量t);ylabel(函数值y);legend(y1=t*t,y2=-t*t,y3=t*t*sin(t);2.曲线拟合与插值一组数据 .447 2.978 1.11 4.25 3.02 6.66 6.01 4.58 8.45 10.35 .22，进行数据拟合并绘制原数据以及拟合曲线一组数据 -.447 1.978 3.11 5.25 5.02 4.66 4.01 4.58 3.45 5.35 9.22，利用样条插值的办法计算其他中间点并绘制图形function t2()x=0:10;y=.447 2.978 1.11 4.25 3.02 6.66 6.01 4.58 8.45 10.35 8.22;p=polyfit(x,y,8); t=0:0.1:10;y1=polyval(p,t);plot(x,y,ro,t,y1,b);grid on;xlabel(自变量x);ylabel(函数值y);title(第二小题曲线拟合并作图);legend(原始数据,拟合曲线);function t2()x0=0:1:10;y0=-.447 1.978 3.11 5.25 5.02 4.66 4.01 4.58 3.45 5.35 9.22;x1=0:0.1:10;y1=spline(x0,y0,x1);plot(x0,y0,bs,x1,y1, k);grid on;xlabel(自变量x);ylabel(函数值y);title(第二小题样条插值并作图);legend(原始数据,插值曲线);hold on;x0=input(请输入需要计算的中间点（010） x0=)y0 = interp1(x,y,x0,spline);plot(x0,y0,ro);legend(插值曲线,原始数据点,需要计算的其他中间点,4)disp(y0=);disp(y0);3.信号处理读入文本文件testdat.txt内的数据(采样频率500Hz)，滤掉70Hz的频率，计算功率谱，显示滤波前后的原始数据波形以及频谱。function t3()fid=fopen(testdat.txt,r);signal,number=fscanf(fid,%f,501);fclose(fid);t=0:0.002:1;b,a=butter(8,0.2,high);signal_y=filter(b,a,signal);subplot(2,1,1);plot(t,signal,t,signal_y);title(滤波前后波形);Fy=fft(signal_y,512);P=Fy.*conj(Fy)/512; f=1000*(0:255)/512;subplot(2,1,2);plot(f,P(1:256);title(频谱);function t3()fid=fopen(test.dat,r);signal,number=fread(fid,501,float32);fclose(fid);t=0:0.002:1;b,a=butter(3,0.14);signal_y=filter(b,a,signal);subplot(2,1,1);plot(t,signal,t,signal_y);title(滤波前后波形);Fy=fft(signal_y,512);P=Fy.*conj(Fy)/512; f=1000*(0:255)/512;subplot(2,1,2);plot(f,P(1:256);title(频谱);设计一个合适的FIR滤波器，100Hz-200Hz带通，要求滤波性能要好。显示其频率响应生成一个采样频率2000，由50Hz、120Hz、200Hz正弦信号以及噪声信号组成。计算其功率谱密度并显示原始数据波形以及频谱function t3()b=fir1(30,0.2 0.4);figure(1);freqz(b,1,512);title(频率响应);t=0:0.001:1;x=sin(t*2*pi*350)+sin(t*2*pi*10);y=fftfilt(b,x);figure(2)subplot(2,1,1);plot(x)title(滤波前波形);subplot(2,1,2);plot(y)title(滤波后波形);function t3()t=0:1/2000:1; x=sin(2*pi*50*t)+sin(2*pi*120*t)+sin(2*pi*200*t);y=x+1.5*randn(1,length(t);Y=fft(y,512);P=Y.*conj(Y)/512;%功率谱f=2000*(0:255)/512;subplot(2,1,1);plot(f,P(1:256);%功率谱图形set(gca,ytick,0:25:150);title(频谱); xlabel(采样点);ylabel(频率);gridsubplot(2,1,2);plot(y);title(原数据波形); xlabel(时间t);ylabel(y);axis tight生成y=5sin(2pi*50)+10sin(2pi*120)信号，滤掉100Hz的频率，计算功率谱，显示滤波前后的原始数据波形以及频谱设计一个合适的IIR滤波器，100Hz-200Hz带阻，要求滤波性能要好。显示其频率响应function t3()t=0:0.001:1;y=5*sin(t*2*pi*50)+10*sin(t*2*pi*120);b,a=butter(8,0.1);y1=filter(b,a,y);subplot(2,1,1);plot(t,y,r);xlabel(时间t);ylabel(信号y);title(原始波形);subplot(2,1,2);plot(t,y1,b);axis(0 1 -20 20);xlabel(时间t);ylabel(信号y1);title(滤波后波形);y2=fft(y,1024);p=y2.*conj(y2)/1024;f=1000*(0:511)/1024;y3=fft(y1,1024);p1=y3.*conj(y3)/1024;figure(2);subplot(2,1,1);plot(f,p(1:512);xlabel(频率f);ylabel(功率密度p);title(滤波前功率谱密度);subplot(2,1,2);plot(f,p1(1:512);xlabel(频率f);ylabel(功率密度p);title(滤波后功率谱密度);function t3()b=fir1(30,0.2 0.4, stop);figure(1);freqz(b,1,512);title(频率响应);t=0:0.001:1;x=sin(t*2*pi*350)+sin(t*2*pi*10);y=fftfilt(b,x);figure(2)subplot(2,1,1);plot(x)title(滤波前波形);subplot(2,1,2);plot(y)title(滤波后波形);4.编程验证哥德巴赫猜想：使用input函数读入一个大于的偶数，并输出一个哥德巴赫素数等式例如:6=3+3 100=41+59要求读入一个，计算一个，能够连续执行，直到读入才结束，退出程序计算Fibonacci数列的值 F(n)=F(n-1)+F(n-2), F(1)=F(2)=1使用input函数读入一个大于正整数，并输出Fibonacci数列的值例如：F(6)=8F(10)=55 要求读入一个，计算一个，能够连续执行，直到读入才结束，退出程序。function t4()n=input(输入一个偶数(输入0退出):);while n0 a=fix(n/2); b=n/2; if(ba) disp(输入不为偶数！); n=input(输入一个偶数(输入0退出):); else for m = 1:n/2; if(isprime(m) k = n - m; if(isprime(k) sprintf(%d=%d+%dn,n,m,k) end end end n=input(输入一个偶数(输入0退出):); end endfunction t4()n = input(输入Fibonacci数列的n值：);while n=1 fn = fibon(n); disp(fn); n = input(输入Fibonacci数列的n值：);endfunction fn = fibon(n)if n=1 | n=2 fn = 1;else fn = fibon(n-1) + fibon(n-2);end编写m程序，找出300以内所有满足下两个条件的数(1)各位数字之和为奇数 (2)该数为素数function t4()n=0;for m=2:30