最新matlab答案完美版

1、1已知兩序列 ，計(jì)算兩序列的卷積并繪制其波形。clear;nx=0:6;x1= nx>=0;x2=( nx-5)>=0;x3=.8.A nx; x=x3.*(x1-x2);nh=0:6; h1=nh>=0;h2=(nh-5)>=0; h=h1-h2;y=conv(x,h);ny=0:9;subplot(3,1,1);stem(nx,x);axis(-1,10,-.2,1.1); title('x');grid subplot(3,1,2);stem(nh,h);axis(-1,6,-.2,1.1); title('h');grid sub

2、plot(3,1,3);stem(ny,y(1:10);axis(-1,9,-1,5); title('y');gridP2.已知復(fù)指數(shù)序列，繪制20點(diǎn)該序列的實(shí)部和虛部。n=0:20;x=1.2*exp(1.5+2j)*n);subplot(2,1,1);stem(n,real(x);title('real');gridsubplot(2,1,2);stem(n,imag(x);title('imag');gridP3.編寫長(zhǎng)度為5的中值濾波器程序。原始未受干擾的序列為：sn=3n(0.5) n,加性噪聲信號(hào)dn為隨機(jī)序列，幅度0.4,分別繪制

3、長(zhǎng)度 為 40 的受干擾序列，以及中值濾波器的輸出。>> clear>> n=0:39;s=3*(0.5* n). n;d=0.8*(rand(1,40)-0.5);x=s+d;y=medfilt1(x,5);受干擾 ');grid; 中值輸出 ');grid;subplot(2,1,1);stem(n,x);title(' subplot(2,1,2);stem(n,y);title(' 只有 medfilt1 是中值濾波器函數(shù) 基本語(yǔ)法為Y = MEDFILT1(X,N)% X 為你想平滑的時(shí)間序列， N 為濾波器的階數(shù)，可缺省% Y

4、 為平滑后的序列P4. 已知序列 x1n=2.2,3,1.5,4.2,1.8, x2n= 0.8,1,1.6,0.8 ， xn=x1n?x2n (卷積)，分別繪制序列x1n ,x2n和x n的波形。clear;n1=0:4;n2=0:3;x1=2.2,3,1.5,4.2,1.4;x2=0.8,1,1.6,0.8;x=conv(x1,x2);n=0:7;subplot(3,1,1);stem(n1,x1);title('x1');gridsubplot(3,1,2);stem(n2,x2);title('x2');grid subplot(3,1,3);stem(

5、n,x(1:8);title('x');gridP5編寫4點(diǎn)滑動(dòng)平均濾波器程序。原始未受干擾的序列為：sn=3n(0.8) n,加性噪聲信號(hào)dn為隨機(jī)序列，幅度0.6,受干擾的序列為: xn= sn+ dn ，分別繪制長(zhǎng)度為 40的原始未受干擾的序列，噪聲序列和 受干擾序列，以及滑動(dòng)平均濾波器的輸出。n=0:39; s=3*(0.8* n).A n; d=1.2*(rand(1,40)-0.5);未干擾 ');grid; 噪聲');grid;受干擾 ');grid; 滑動(dòng)輸出 ');grid;x=s+d; windowSize = 4; y=fi

6、lter(ones(1,windowSize)/windowSize,1,x); subplot(4,1,1);stem(n,s);title(' subplot(4,1,2);stem(n,d);title(' subplot(4,1,3);stem(n,x);title(' subplot(4,1,4);stem(n,y);title('P6序列，分別繪制下列序列的波形圖:( 1 ) ； ( 2)( 3)(1) nx=0:7;x1=nx>=0;x2=(nx-8)>=0;x3=.8.A nx;x4=cos(0.75*pi*nx);x=x3.*(x1

7、-x2).*x4;ny=-nx+2; subplot(3,1,1);stem(ny,x);title('yn');grid; h=x;nh=nx+1;nz=0:8;y1=zeros(1,9);y2=y1;y1(find(nz>=0)&(nz<=7)=1)=x;y2(find(nz>=1)&(nz<=8)=1)=h;z=2*y1+y2; subplot(3,1,2);stem(nz,z);title('zn');grid;w=x;nw=2*nx;subplot(3,1,3);stem(nw,w);title('wn&

8、#39;);gridP7序列，分別繪制其偶序列分量 和圓周偶序列分量 的波形圖:nx=0:3;x1=nx>=0;x2=(nx-3)>=0; x=0.75*nx.*(x1-x2);ny=-3:3;y1=zeros(1,7);y2=y1;m=fliplr(x);y1(find(ny>=0)&(ny<=3)=1)=x; y2(find(ny>=-3)&(ny<=0)=1)=m;偶序列 ');grid;圓偶序列 ');grid;y=(y1+y2)/2; subplot(2,1,1);stem(ny,y);title('z=(m

9、+x)/2;subplot(2,1,2);stem(nx,z);title('P8設(shè)原始未受干擾的序列為：sn=3n(0.5) n,加性噪聲信號(hào)dn為隨機(jī)序列, 幅度 0.4，對(duì)該信號(hào)測(cè)量 50次并取樣本的整體平均以去除噪聲，分別繪制某受干擾的長(zhǎng)度為 40的測(cè)量序列，以及整體平均濾波的輸出。n=0:39;s=3*(0.5* n，n;d=0.8*(rand(1,40)-0.5);x=s+d;sub plot(2,1,1);stem(n,x);title(受干擾');grid; np=50;h=1/np;h=h(ones(1,np);y1=conv(h,x);subplot(2,1

10、,2);stem(n,y1);title(除噪信號(hào)');grid;P9繪制指數(shù)序列 的幅度譜和相位譜曲線。n=0:5000;x1=n>=0;x=(0.8.A n).*x1;w=0:0.01*pi:pi;H=freqz(x,1,w); subplot(2,1,1);plot(w/pi,abs(H);grid; subplot(2,1,2);plot(w/pi,angle(H);grid;P10.繪制矩形序列 的幅度譜和相位譜曲線。n=0:10;x1=n>=0;x2=(n-10)>0;h=0.5*(x1-x2);w=0:0.01*pi:pi;H=freqz(h,1,w);

11、sub plot(2,1,1); plot(w/pi,abs(H);grid; sub plot(2,1,2); plot(w/pi,a ngle(H);grid;P11.已知濾波器的差分方程和輸入信號(hào)分別為:yn= - 6.76195xn+13.456335xn-1 - 6.76195xn- 2xn =cos(0.1 n)+2cos(0.4 n)u n繪制該系統(tǒng)的輸入序列和輸出序列。n=0:30;x1=( n>=0);b=-6.76195,13.456335,-6.76195;a=1;x=(cos(0.1* n)+2*cos(0.4* n) ).*x1;y=filter(b,a,x);

12、subplot(2,1,1);stem(n,x);title(輸入序列');grid;subplot(2,1,2);stem(n,y);title(輸出序列');grid;P12.繪制8點(diǎn)滑動(dòng)平均濾波器的幅頻特性和相頻特性。n=0:7;x1= n>=0;x2=( n-7)>0; h=(1/8)*(x1-x2);w=0:0.01* pi:pi;H=freqz(h,1,w);幅頻特性');grid;相頻特性');grid;sub plot(2,1,1); plot(w/pi,abs(H);titleC sub plot(2,1,2); plot(w/pi

13、,a ngle(H);title('P13.設(shè)計(jì)4階巴特沃茲模擬低通濾波器，其3-dB截止頻率為500兀 繪制濾波器的幅頻響應(yīng)曲線。b,a=butter(4,500* pi ,'s');h,w=freqs(b,a);plot(w/pi,20*log10(abs(h);titleC幅頻特性');grid;P14設(shè)計(jì)切比雪夫I型4階模擬低通濾波器，其3-dB截止頻率為500!，繪制 濾波器的幅頻響應(yīng)曲線。b,a=cheby1(4,1,500* pi,'s');h,w=freqs(b,a);plot(w/pi,20*log10(abs(h);title

14、C幅頻特性');grid;P15.設(shè)計(jì)巴特沃茲模擬高通濾波器，其性能指標(biāo)為：歸一化通帶邊緣頻率為Op=8000 n rad/s)s = 2000 n rad/s通帶波紋為0.1dB，最小阻帶衰減為40dB。 繪制所設(shè)計(jì)的濾波器增益響應(yīng)。wp=1;ws=4;afa p=0.1;afas=40;n,wn=buttord(wp,ws,afap,afas,'s');b,a=butter(n,wn,'s');num,den=lp2hp(b,a,2*pi*4000);h,w=freqs(num,den);plot(w/(1000*pi),20*log10(abs(h

15、);title('幅頻特性 ');grid;P16.繪制序列 的16點(diǎn)DFTn=0:10;x1=n>=0;x2=(n-10)>0;x=x1-x2;xk=fft(x,16);m=0:15;stem(m,xk);title('DFT ');grid;P17.已知序列，其中，。繪制該信號(hào)的幅度譜曲線及計(jì)算該序列16點(diǎn)的DFT。n=0:15;x1=n>=0;x2=(n-15)>0;x=cos(2*pi*n/8).*(x1-x2);h,w=freqz(x,1);subplot(3,1,1);plot(w/pi,abs(h);title('x

16、k=fft(x,16);m=0:15;subplot(3,1,2);stem(m,abs(xk);title('DFT subplot(3,1,3);stem(m,angle(xk);title('DFT幅頻曲線 ');grid;幅度');grid; 相位');grid;P18.序列,繪制該序列的8點(diǎn)圓周移位序列。x=5,4,3,2,1,0;x=x zeros(1,8-length(x);n=0:7;y=x(mod(n-4,8)+1);stem(n,y);titleC移位后序列);gridP19.序列,計(jì)算并繪制兩序列的8點(diǎn)圓周卷積。x=5,4,3,2,

17、1,0;h=0.5,0.5,0.5,0.5,0.5,0.5; x=x,zeros(1,8-length(x); h=h,zeros(1,8-length(h); m=0:7;hm=h(mod(-m,8)+1);H=toeplitz(hm,0,h(2:8);y=x*H;stem(m,y);titleC八點(diǎn)圓周卷積序列);gridP20.已知序列，利用DFT計(jì)算并繪制兩序列的卷積。 nx=0:15;x1=nx>=0;x2=(nx-15)>0; x=x1-x2;h=x;j=conv(h,x);x=h,zeros(1,31);h=h,zeros(1,31);X=fft(x,31);H=ff

18、t(h,31);Y=X.*H;y=ifft(Y,31);n=0:30;stem(n,y);title(' 線性卷積和 ');gridP21.已知序列,分別繪制該信號(hào)的16點(diǎn)DFT及32點(diǎn)DFT的幅度譜序列和相位譜序列。n=0:15;x1=n>=0;x2=(n-15)>0;x=(cos(pi*n/8)+2*cos(pi*n/2).*(x1-x2);xk=fft(x,16);m=0:15;sub plot(2,2,1);stem(m,abs(xk);titleCDFT 幅度');grid; sub plot(2,2,2);stem(m,angle(xk);tit

19、leCDFT 相位');grid;yk=fft(x,32);b=0:31;sub plot(2,2,3);stem(b,abs(yk);titleCDFT 幅度');grid;sub plot(2,2,4);stem(b,angle(yk);titleCDFT 相位');grid;P22.分別繪制序列 的16點(diǎn)DFT樣本的實(shí)部和虛部序列分量。n=0:15;x1=n>=0;x2=(n-15)>0;x=(0.5A n).*(x1-x2);xk=fft(x,16);實(shí)部');grid; 虛部');grid;m=0:15;subplot(2,2,1)

20、;stem(m,real(xk);title('DFT subplot(2,2,2);stem(m,imag(xk);title('DFTP23.某信號(hào)的頻譜序列，繪制該序列時(shí)域波形。n=0:15;x1=n>=0;x2=(n-15)>0;x=(n/8).*(x1-x2);y=ifft(x);stem(n,x) title(' 時(shí)域波形 ');gridP24某系統(tǒng)的系統(tǒng)函數(shù)為計(jì)算該系統(tǒng)的直接形式的結(jié)構(gòu)參數(shù)。r = 1,2,1;p = 2,2.5,2.5;k = 1;num, den = residuez(r,p,k);disp( 'Numera

21、tor polynomial coefficients' ); disp(num) disp( 'Denominator polynomial coefficients' ); disp(den)P25.某離散系統(tǒng)的極點(diǎn)為零點(diǎn)為 ，計(jì)算該系統(tǒng)的直接形式的結(jié)構(gòu)參數(shù)。format longzr = -1 2+4*i 2-4*i 0.5;pr = -0.8 0.5+0.4*i 0.5-0.4*i 0.6i; z = zr' p = pr'num, den = zp2tf(z, p, 1);disp( 'H(z)=' );disp(num); d

22、isp(den);P26.計(jì)算序列 的z變換，并在z平面上繪制其零極點(diǎn)圖。P27.繪制IIR系統(tǒng) 的幅頻與相頻響應(yīng)。b=1,-2,4;a=1,2.37,2.7,1.6,0.41;h,w=freqz(b,a);sub plot(2,1,1); plot(w/ pi,abs(h);title(幅頻');gridsub plot(2,1,2); plot(w/ pi,angle(h);title (相頻');gridP28.繪制FIR系統(tǒng) 的幅頻與相頻響應(yīng)。b=1,-2,1;h,w=freqz(b,1);sub plot(2,1,1); plot(w/pi,abs(h);title(

23、幅頻');gridsub plot(2,1,2); plot(w/ pi ,angle(h);title(相頻');gridP29.計(jì)算序列 的z變換，并在z平面上繪制其零極點(diǎn)圖。syms n z;x=2.A n;xz=ztrans(x )；xz =1/2*z/(1/2*z-1)z=0;p=2;zplane(z,p);P30.繪制IIR系統(tǒng) 的幅頻與相頻響應(yīng)。b=1,-2;a=1,0.5,0,-0.25;h,w=freqz(b,a);sub plot(2,1,1); plot(w/ pi,abs(h);title(幅頻');gridsub plot(2,1,2); pl

24、ot(w/ pi,angle(h);title(相頻');gridP31 .繪制FIR系統(tǒng) 的幅頻與相頻響應(yīng)。b=1,-2,0,-2,1;h,w=freqz(b,1);sub plot(2,1,1); plot(w/pi,abs(h);title('幅頻');gridsub plot(2,1,2); plot(w/pi,a ngle(h);title('相頻');gridP32.設(shè)計(jì)切比雪夫I型IIR數(shù)字高通濾波器，其性能指標(biāo)為：通帶波紋g=0.5dB , 最小阻帶衰減 a=43dB，通帶和阻帶邊緣頻率分別為卯=0.75冗rad和ws=0.35冗rad。

25、繪制所設(shè)計(jì)的濾波器增益響應(yīng)。wp=0.75;ws=0.35;Rp=0.5;As=43; N,wc=cheb1ord(wp,ws,Rp,As); b,a=cheby1(N,Rp,wc, 'high' ); h,w=freqz(b,a);P lot(w/ pi,20*log10(abs(h);titleC 增益');grid;P33 設(shè)計(jì)巴特沃茲數(shù)字帶通濾波器，其指標(biāo)要求為：通帶邊緣頻率， ，阻帶邊 緣頻率， ，通帶波紋為 1dB, 阻帶衰減為 40dB, 計(jì)算該濾波器的階數(shù)及 3dB截止頻率。wp=0.45,0.65;ws=0.3,0.75;Rp=1;rs=40;N,wc

26、=buttord(wp,ws,Rp,rs);P34.設(shè)計(jì)I型橢圓IIR數(shù)字低通濾波器，其性能指標(biāo)為：通帶邊緣頻率0.5 n阻帶邊緣頻率0.75 n通帶波紋為1dB，最小阻帶衰減為31dB。計(jì)算該濾 波器的階數(shù)N及3dB截止頻率。wp=0.5,ws=0.75;Rp=1;As=31;N,wc=elli pord(w p, ws,R p,As); disp (N);dis p( wc);P35.設(shè)計(jì)FIR數(shù)字低通濾波器，其指標(biāo)要求為：卯=0.3 n s =0.4 n減為啟=50dB。用凱塞窗設(shè)計(jì)，畫出其幅頻響應(yīng)曲線。阻帶衰fpts=0.3 0.4;As=50;ds=10A(-As/20) mag=1

27、,0;dev=ds ds;N,w n,beta,ft yp e=kaiserord(fpts,mag,dev);b=fir1(N,w n,kaiser(N+1,beta);h,w=freqz(b,1);p lot(w/pi,20*log(abs(h);P36.設(shè)計(jì)FIR數(shù)字高通濾波器，其指標(biāo)要求為：ws = 0.4 np =j0.6 n減為氏=50dB。用凱塞窗設(shè)計(jì)，畫出其幅頻響應(yīng)曲線。阻帶衰fpts=0.4 0.6;As=50;ds=10A(-As/20) mag=0 1;dev=ds ds;N,w n,beta,ft yp e=kaiserord(fpts,mag,dev); b=fir1

28、(N,w n, 'high' ,kaiser(N+1,beta); h,w=freqz(b,1);p lot(w/ pi,20*log(abs( h) );title('幅頻');gridP37.設(shè)計(jì)FIR數(shù)字帶通濾波器，其指標(biāo)要求為：卯1 = 0.3 np2 c= 0.60.2 n， ws2 = 0.75，阻帶衰減為5s = 55dB。用凱塞窗設(shè)計(jì)，畫出其幅頻 響應(yīng)曲線。n s1(=fpts=0.2 0.3 0.6 0.75;As=55;ds=10A(-As/20) mag=0 1 0;dev=ds ds ds;N,w n,beta,ft yp e=kaiserord(fpts,mag,dev); b=fir1(N,w n,ftyp e,kaiser(N+1,beta);h,w=freqz(b,1);p lot(w/ pi,20*log(abs(h);titleC 幅頻')；gridP38.對(duì)10赫茲的正弦信號(hào)