信號處理分析課程設計報告(編輯修改稿)
2025-05-08 22:34 本頁面
【文章內(nèi)容簡介】 gth(z0)/2+1))。frq0 = [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。輸出圖:圖31 原始信號波形及其頻譜圖32 已調(diào)信號波形及其頻譜 噪聲部分代碼:%噪聲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。)。輸出圖:圖33 加入噪聲后信號波形圖及其頻譜 帶通濾波部分代碼：%帶通濾波器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。)。輸出圖：圖34 帶通濾波后信號波形及其頻譜 解調(diào)部分代碼：%解調(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。)。輸出圖：圖35 解調(diào)后信號波形圖及其頻譜 低通濾波部分代碼：%低通濾波器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(B
環(huán)評公示相關推薦
鄂ICP備17016276號-1