信號處理分析課程設(shè)計報告-資料下載頁
2025-04-11 22:34本頁面
【正文】 [0:length(z0)1]*Fs/length(z0)/2。figure(1)。subplot(2,1,1)。 plot(x)。 grid。title(39。原始信號39。)。subplot(2,1,2)。 plot(frq0,z0)。 grid。title(39。原始信號的頻譜39。)。%雙邊帶調(diào)制y1 = ammod(x,Fc,Fs)。z1 = fft(y1)。z1 = abs(z1(1:length(z1)/2+1))。frq1 = [0:length(z1)1]*Fs/length(z1)/2。figure(2)。subplot(2,1,1)。 plot(y1)。 grid。title(39。雙邊帶信號39。)。subplot(2,1,2)。 plot(frq1,z1)。title(39。雙邊帶信號的頻譜39。)。grid。%噪聲noisy=randn(1,N)。 y2=y1+noisy/2。z2=fft(y2)。%frq2=[0:N1]*Fs/N。frq2 = [0:length(z2)1]*Fs/length(z2)。figure(3)。subplot(2,1,1)。 plot(y2)。 grid。title(39。加入噪聲后的信號39。)。subplot(2,1,2)。 plot(frq2,abs(z2))。 grid。title(39。加入噪聲后信號的頻譜39。)。%帶通濾波器rp=1。rs=10。 wp=2*pi*[45,55]。ws=2*pi*[40,60]。[N,wc]=buttord(wp,ws,rp,rs,39。s39。)。[B,A]=butter(N,wc,39。s39。)。[Bz,Az]=impinvar(B,A,Fs)。y3=filter(Bz,Az,y2)。z3=fft(y3)。figure(4)。subplot(2,1,1)。 plot(y3)。 grid。title(39。帶通濾波后的的信號39。)。subplot(2,1,2)。 plot(frq2,abs(z3))。 grid。title(39。帶通濾波后信號的頻譜39。)。%解調(diào)y4 = amdemod(y3,Fc,Fs,0,0)。z4=fft(y4)。z4 = abs(z4(1:length(z4)/2+1))。frq4 = [0:length(z4)1]*Fs/length(z4)/2。figure(5)。subplot(2,1,1)。 plot(y4)。 grid。title(39。解調(diào)后的信號39。)。subplot(2,1,2)。 plot(frq4,z4)。grid。title(39。解調(diào)后的信號的頻譜39。)。%低通濾波器fp1=1。fs1=4。rp1=1。rs1=10。 wp1=2*pi*fp1。ws1=2*pi*fs1。 [N1,wc1]=buttord(wp1,ws1,rp1,rs1,39。s39。)。 [B1,A1]=butter(N1,wc1,39。s39。)。 [Bz1,Az1]=impinvar(B1,A1,Fs)。 y5=filter(Bz1,Az1,y4)。 z5=fft(y5)。z5 = abs(z5(1:length(z5)/2+1))。frq5 = [0:length(z5)1]*Fs/length(z5)/2。figure(6)。subplot(2,1,1)。 plot(y5)。 grid。title(39。低通濾波后的信號39。)。subplot(2,1,2)。 plot(frq5,z5)。grid。title(39。低通濾波后信號的頻譜39。)。20
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1