1、傳播優(yōu)秀word版文檔 ，希望對您有幫助，可雙擊去除！matlab上機模擬試題1、 計算當時，的值。2、 計算3、 求解矩陣方程： 4、 求解矩陣方程：5、 計算序列-2 0 1 1 3和序列1 2 0 -1的離散卷積。 6、 用matlab計算差分方程 p138 157當輸入序列為時的輸出結果。n=41;a=0.8 -0.44 0.36 0.22; b=1 0.7 -0.45 -0.6; x=1 zeros(1,n-1); k=0:1:n-1; y=filter(a,b,x); stem(k,y) xlabel('n');ylabel('幅度') 7、用mat

2、lab計算差分方程 p156所對應的系統(tǒng)函數(shù)的頻率響應。b=0.8 -0.44 0.36 0.02;a=1 0.7 -0.45 -0.6;freqz(b,a);8、 用fft計算序列的頻譜。xn=1 3 -1 1 2 3 3 1;xn=fft(xn);figure(1);plot(abs(xn);9、求下列直接型系統(tǒng)函數(shù)的零、極點。b=1,-0.1,-0.3,-0.3,-0.2;a=1,0.1,0.2,0.2,0.5;disp('h(z)');printsys(b,a,'z');z,p,k=tf2zp(b,a)10、畫出序列 的圖形并補零，序列長度為20.xn=

3、1 2 3 4 5;b=zeros(1,15);xn=xn,b;stem(xn)11、 產(chǎn)生一個1024點的正弦加白噪聲序列，并繪其頻譜。x=0:1023;y=sin(x);y1=awgn(y,10); freqz(y1)12、 已知信號，求其fft變換的幅值和相位值，并繪出結果。y=2.*(cos(2.*x).*(cos(100.*x);y1=fft(y);freqz(y1)13、給定，。分別求兩序列相加、相卷積。用4個子圖顯示。n1=0:4;x1=1 2 3 4 5;subplot(4,1,1),stem(n1,x1);ylabel('x1');n2=0:3;x2= 1 2

4、 1 2;subplot(4,1,2),stem(n2,x2);ylabel('x2');n=0:4;x1=x1 ,zeros(1,5-length(n1);x2= x2,zeros(1,5-length(n2);x=x1+x2;y=conv(x1,x2);subplot(4,1,3),stem(n,x);ylabel('x');subplot(4,1,4),stem(y);ylabel('y')14、 生成一個長度為50的漢寧窗，并觀察期頻率特性。n=51;window=hanning(n);h,w=freqz(window,1);subplo

5、t(1,2,1)stem(window);subplot(1,2,2);plot(w/pi,20*log(abs(h)/abs(h(1);15、已知有限長序列x(n)=7,6,5,4,3,2,求dft和idft，要求：畫出序列傅立葉變換對應的幅度譜和相位譜；畫出原信號與傅立葉逆變換idftx(k)的圖形進行比較。title('|x(k)|');subplot(2,2,4),stem(k,angle(xk);title('arg|x(k)|');xn=7,6,5,4,3,2;n=length(xn);n=0:(n-1);k=0:(n-1);xk=xn*exp(-j

6、*2*pi/n).(n'*k);x=(xk*exp(j*2*pi/n).(n'*k)/n;subplot(2,2,1),stem(n,xn);title('x(n)');subplot(2,2,2),stem(n,abs(x);title('idft|x(k)|');subplot(2,2,3),stem(k,abs(xk);16、 已知周期序列的主值x(n)=7,6,5,4,3,2,求x(n)周期重復次數(shù)為3次時的dfs和idfs。要求：畫出原信號序列的主值和周期序列的圖形；畫出離散傅立葉變換對應的幅度譜和相位譜. xn=7,6,5,4,3,2

7、; n=length(xn); n=0:3*n-1;k=0:3*n-1;xn1=xn(mod(n,n)+1);xk=xn1*exp(-j*2*pi/n).(n'*k);subplot(2,2,1),stem(xn);title('原主值信號x(n)');subplot(2,2,2),stem(n,xn1);title('周期序列信號');subplot(2,2,3),stem(k,abs(xk);title('|x(k)|');subplot(2,2,4),stem(k,angle(xk);title('arg|x(k)|'

8、;);17、 求x(n)=7,6,5,4,3,2, 0=<n<=5的dtft，將（-2*pi，2*pi）區(qū)間分成500份。要求：畫出原信號；畫出離散傅立葉變換對應的幅度譜和相位譜；求有限長序列x(n)=7,6,5,4,3,2在n=100時的dft，并與dtft進行對比。xn=7,6,5,4,3,2;n=length(xn);n=0:n-1;w=linspace(-2*pi,2*pi,500);x=xn*exp(-j*n'*w);subplot(3,1,1),stem(n,xn,'k');ylabel('x(n)');subplot(3,1,2

9、),plot(w,abs(x),'k');axis(-2*pi,2*pi,1.1*min(abs(x),1.1*max(abs(x);ylabel('幅度譜');subplot(3,1,3),plot(w,angle(x),'k');axis(-2*pi,2*pi,1.1*min(angle(x),1.1*max(angle(x);ylabel('相位譜');（2）n=100;xn=7,6,5,4,3,2,zeros(1,n-6);n=0:(n-1);k=0:(n-1);xk=xn*exp(-j*2*pi/n).(n'*k

10、);x=(xk*exp(j*2*pi/n).(n'*k)/n;subplot(2,1,1),stem(k,abs(xk);title('|x(k)|');subplot(2,1,2),stem(k,angle(xk);title('arg|x(k)|');18、繪出下列直接型系統(tǒng)函數(shù)的頻率響應。b=1 -0.1 -0.3 -0.3 -0.2;a=1 0.1 0.2 0.2 0.5;freqz(b,a)19、 設采樣周期t=250s（采樣頻率fs =4khz），用脈沖響應不變法和雙線性變換法設計一個三階巴特沃茲濾波器，并繪出頻率響應，其3db邊界頻率為fc

11、 =1khz。b,a=butter(3,2*pi*1000,'s'); num1,den1=impinvar(b,a,4000); h1,w=freqz(num1,den1); b,a=butter(3,2/0.00025,'s'); num2,den2=bilinear(b,a,4000); h2,w=freqz(num2,den2); f=w/pi*2000; plot(f,abs(h1),'-.',f,abs(h2),'-'); grid; xlabel('頻率/hz ') ylabel('幅值/db

12、') 程序中第一個butter的邊界頻率2×1000，為脈沖響應不變法原型低通濾波器的邊界頻率；第二個butter的邊界頻率2/t=2/0.00025，為雙線性變換法原型低通濾波器的邊界頻率.圖3.1給出了這兩種設計方法所得到的頻響，虛線為脈沖響應不變法的結果；實線為雙線性變換法的結果。脈沖響應不變法由于混疊效應，使得過渡帶和阻帶的衰減特性變差，并且不存在傳輸零點。同時，也看到雙線性變換法，在z=-1即=或f=2000hz處有一個三階傳輸零點，這個三階零點正是模擬濾波器在=處的三階傳輸零點通過映射形成的。20、 用漢寧窗設計一個線性相位fir濾波器，已知，n=37. 并繪出頻

13、率響應。wc=0.3*pi;n=37;b=fir1(n,wc,hanning(n+1);freqz(b);ylabel('頻率響應');b=fir1(37,0.3);a=1;freqz(b,a)21、 用頻率采樣法設計一個線性相位fir濾波器，已知，n=37. 并繪出頻率響應。n=37;f=0:1/37:1; (2)a=ones(1,18),zeros(1,n-17);b=fir2(n,f,a);freqz(b)22、 生成周期為50hz和200hz的合成正弦信號，設計模擬濾波器濾除50hz的正弦信號，給出濾波前后信號的對比圖，并繪出濾波器頻率響應。wp=40 60/100;w

14、s=49 51/100;rp=3;rs=20;fs=200;n,wn=buttord(wp,ws,rp,rs);h=butter(n,wn,'stop');figure(1),freqz(h,512,fs);t=0:0.001:0.05;f1=50;f2=200;x1=sin(2*pi*f1*t);x2=sin(2*pi*f2*t)s=x1+x2;sf=filter(h,1,s);figure(2),subplot(2,1,1),plot(t,s);ylabel('幅值');xlabel('hz');title('混合正弦信號')

15、;figure(2),subplot(2,1,2),plot(t,sf);ylabel('幅值');xlabel('hz');title('濾波后的信號');y=fft(s,512);pyy=y.*conj(y)/512;f=1000/512*(0:255);figure(3),subplot(2,1,1),plot(f/5,pyy(1:256);ylabel('幅值');xlabel('hz');title('濾波前的頻譜');y=fft(sf,512);pyy=y.*conj(y)/512;f=

16、1000/512*(0:255);figure(3),subplot(2,1,2),plot(f/5,pyy(1:256);ylabel('幅值');xlabel('hz');title('濾波后的頻譜')23、 生成周期為50hz和200hz的合成正弦信號，設計fir濾波器濾除50hz的正弦信號，采樣頻率為1000hz，給出濾波前后信號的對比圖。t=0:0.001:0.05fs=1000;f1=50;f2=200;x1=sin(2*pi*f1*t);x2=sin(2*pi*f2*t);x0=x1+x2;figure(1);plot(t,x0);xlabel('t/s');ylabel('幅值');fm=fs/2;wc=75/fm;b=fir1(24,wc,'high');y0=filter(b,1,x0);figure(2);plot(t,y0);xlabel('t/s');ylabel('幅值');24、 生成周期為50hz、100 hz和200hz的合成正弦信號，設計fir濾波器濾除100hz的正弦信號，采樣頻率為1000hz，給出濾波前后信號的對比圖。t=0:0.001:0.05fs=1000;f1=50;f2=100;f3=20