新·matlab作业(齐).docVIP

作业１ 设,,求使. 程序及运行结果如下: 作业２ 利用公式计算的值. 程序： i = 0; pi = 0; while (2*i+1)=21 pi = pi+(-1)^i/(2*i+1); i = i+1; end pi = 4*pi; fprintf(pi = %f\n,pi); 运行结果为： pi = 3.232316 作业３ 在同一坐标系下绘制,并标注. 程序： t = 0:.1:2*pi; y1 = t.^2; y2 = -y1; y3 = y1.*sin(t); plot(t,y1,-r,t,y2,:b,t,y3,-.k); grid legend(t^2,-t^2,t^2*sin(t)); 运行结果： 作业４ 在一个图形窗口绘制半径为1的球面,的平面以及马鞍面. 程序: % 半径为1的球面 subplot(2,2,1) x0 = 0; y0 = 0; z0 = 0; r = 1; [x,y,z] = sphere; surf(r*x+x0,r*y+y0,r*z+z0) axis equal % z=4 的平面 subplot(2,2,2) a = [-10,10]; b = [-10,10]; [x,y] = meshgrid(a,b); z = 4*ones(size(x)); surf(x,y,z) % 马鞍面z=2*x^2-y^2 subplot(2,2,3) [x,y] = meshgrid(-10:0.2:10); z = 2*x.^2-y.^2; surf(x,y,z) 运行结果： 作业5 设某信号,其中为白噪声.求信号的DFT及其反DFT变换,画出原始信号,幅频特性和反变换后信号的曲线. 程序： N = 8; f = 2; fs = 10; n = 0:N-1; xn = 5*sin(2*pi*f*n/fs)+2*randn(size(n)); WN = exp(-i*2*pi/N); k = 0:N-1; nk = n*k; WNnk = WN.^nk; XK = xn*WNnk; WNkn = WN.^(-nk); xn1 = XK*WNkn/N; [h,w] = freqz(xn); [p,q] = freqz(xn1); subplot(2,2,1);stem(n,xn);title(原始信号) subplot(2,2,2);plot(w/pi,abs(h));title(原始幅频特性) subplot(2,2,3);stem(n,xn1);title(IDFT变换后的信号) subplot(2,2,4);plot(q/pi,abs(p));title(IDFT变换后得幅频特性); 运行结果: 作业6 设计IIR滤波器（带通，三种方法，fs=2000HZ,通带频率300~500HZ，阶数自选，画频率特性并分析比较）. (1)第一种方法：采用butterworth滤波器，程序如下： n=10; wn=[300 500]/1000; [b,a]=butter(n,wn); freqz(b,a,512,2000);title(10阶带通butterworth滤波器); 运行结果: (2)第二种方法：采用椭圆滤波器，程序如下： n=10; wn=[300 500]/1000; [b,a]=ellip(n,0.5,20,wn); freqz(b,a,128,2000);title(10阶带通椭圆滤波器); 运行结果: (3)第三种方法:采用切比雪夫滤波器，程序如下： n=10; wn=[300 500]/1000; [b,a]=cheby1(n,0.5,wn); freqz(b,a,128,2000);title(10阶带通cheby1滤波器); 运行结果: 分析比较：Butterworth带通滤波器的通带内的幅值响应比较平坦，但是截止频率处的截止特性不是非常理想。Chebyshev I型带通滤波器通带内的特性也很好，缺点也是截止特性不是很好。与前两种滤波器相比，椭圆带通滤波器在截止频率处的截止特性比较好. 作业7 对一幅分别添加了椒盐噪声和高斯白噪声的图象进行均值和中值滤波. (1)对添加了椒盐噪声得图像进行均值和中值滤波,程序如下: a = imread(cameraman.tif); b = imnoise(a,salt pepper); c = filter2(fspecial(average,5),b)/255; d = medfilt2(b); subplot(2,2,1); imshow(a); title(原始图像); subplot(2,2,2); imshow(b); title(椒