【正文】 。% hamming windowb5=fir1(M,wc/pi,39。high39。,hamming(N+1))。[H5,w]=freqz(b5,wc,512)。H5_db=20*log10(abs(H5))。% hanning windowb6=fir1(M,wc/pi,39。high39。,hanning(N+1))。[H6,w]=freqz(b6,wc,512)。H6_db=20*log10(abs(H6))。figure(3)c=plot(w,H4_db,w,H5_db,39。y39。,w,H6_db,39。r39。)。figure(4)subplot(3,1,1)。stem(n1,real(h4))。axis([0 ])。subplot(3,1,2)。stem(n1,real(h5),39。y39。)。axis([0 ])。subplot(3,1,3)。stem(n1,real(h6),39。r39。)。axis([0 ])。x=1+2*cos(pi/4*n1)+cos(pi/2*n1)。y1=conv(x,h1)。y2=conv(x,h2)。y3=conv(x,h3)。y4=conv(x,h4)。y5=conv(x,h5)。y6=conv(x,h6)。figure(5)subplot(3,2,1)。stem(n2,y1)。axis([0 1])。subplot(3,2,2)。stem(n2,y2)。axis([0 1])。subplot(3,2,3)。stem(n2,y3)。axis([0 1])。subplot(3,2,4)。stem(n2,y4)。axis([0 502 2])。subplot(3,2,5)。stem(n2,y5)。axis([0 502 2])。subplot(3,2,6)。stem(n2,y6)。axis([0 502 2])。***************03040502010.******04050第五章1．某隨機信號由兩余弦信號與噪聲構(gòu)成：x[k]= cos(20πk)+cos(40πk)+ s [k]，s[k]為均值為0，方差為1的高斯白噪聲。(1)繪出此隨機信號的時域波形；(2)試分別用周期圖法、平均周期圖法和Welch法分析該序列的功率譜估計。Fs = 1000。% 抽樣頻率 t = 0:1/Fs:1。% 抽樣時間xn = cos(20*pi*t)+ cos(40*pi*t)+ randn(size(t))。%粗略地估計xn的功率譜，做N=1024點FFT：Pxx = abs(fft(xn,1024)).^2/1001。subplot(3,3,1)。plot(t,xn)。xlabel(39。隨機信號39。)。grid on。subplot(3,3,2)。plot([0:1023]*Fs/1024,10*log10(Pxx))。xlabel(39。利用公式39。)。grid on。window=boxcar(1001)。[Pxx1,F1] = periodogram(xn,window,1024,Fs)。subplot(3,3,3)。plot(F1, 10*log10(Pxx1))。xlabel(39。利用函數(shù)periodogram39。)。grid on。noverlap=500。[Pxx2,F2] = psd(xn, 1024,Fs, window, noverlap)。subplot(3,3,4)。plot(F2, 10*log10(Pxx2))。xlabel(39。利用函數(shù)psd39。)。grid on。noverlap=500。[Pxx3,F3] = pwelch(xn, window39。, noverlap, 1024,Fs)。subplot(3,3,5)。plot(F3, 10*log10(Pxx3))。xlabel(39。利用函數(shù)pwelch39。)。grid on。第五篇：通信原理課程設(shè)計課設(shè)一一、設(shè)計題目信號特性分析（如正弦波信號的波形與頻譜）二、設(shè)計目的通信原理課程設(shè)計是《通信原理》理論課的輔助環(huán)節(jié)。著重體現(xiàn)通信原理教學知識的運用，、擴展通信原理所學知識。其主要通過 matlab 仿真進一步深化對通信原理知識的學習。三、設(shè)計內(nèi)容采用matlab產(chǎn)生不同頻率，不同幅度的兩種正弦波信號，并將這兩個信號疊加為一個信號，觀察這三個信號的波形。對疊加后的信號用FFT作譜分析。要求：繪出正弦信號的時域波形掌握傅立葉變換及其逆變換利用傅立葉變換繪出正弦信號的頻譜四、實驗原理正弦序列x(n)=As2ipfnn/Fs(+j)，在MATLAB中n=0:N1x=A*sin(2*pi*f*n/Fs+fai)將信號源發(fā)出的信號強度按頻率順序展開，使其成為頻率的函數(shù)，并考察變化規(guī)律，稱為頻譜分析。對信號進行頻譜分析，往往對其進行傅里葉變換，觀察其頻譜幅度與頻譜相位。對于信號來說，分模擬信號與數(shù)字信號。對于模擬信號來說，往往對其進行抽樣，然后進行快速傅里葉變換（fft），然后對其幅度（abs）和相位（angle）的圖像進行分析。對于數(shù)字信號，則可直接進行快速傅里葉變換。五、程序截圖六、源程序代碼clear all clc。f1=100。%信號頻率Hz f2=150。%信號頻率Hz fs=1000。%采樣頻率Hz N=20。%采樣點數(shù)t=(0:N1)/fs。%采樣時間s x1=5*sin(2*pi*f1*t)。%信號采樣值 x2=10*sin(2*pi*f2*t)。%信號采樣值 subplot(231)。stem(t,x1,39。.39。)subplot(232)。stem(t,x2,39。.39。)。subplot(233)。stem(t,x1+x2,39。.39。)。y1=fft(x1,512)。subplot(234)。plot(abs(y1).^2)。xlabel(39。頻率(Hz)39。)。ylabel(39。幅值39。)。y2=fft(x1,512)。subplot(235)。plot(abs(y2).^2)。xlabel(39。頻率(Hz)39。)。ylabel(39。幅值39。)。y3=fft(x1+x2,512)。subplot(236)。plot(abs(y3).^2)。xlabel(39。頻率(Hz)39。)。ylabel(39。幅值39。)。課設(shè)二一、設(shè)計題目正弦波信號的波形與頻譜分析AM模擬調(diào)制二、設(shè)計目的熟悉matlab的編程環(huán)境及使用；學會利用matlab進行信號處理及分析；掌握傅立葉變換及其逆變換；學會利用傅立葉變換繪出正弦信號的頻譜；學會用matlab產(chǎn)生特定頻率及功率的正弦信號；學會利用matlab對信號進行載波、解調(diào)處理；三、設(shè)計要求信號特性分析（如正弦波信號的波形與頻譜）采用matlab產(chǎn)生不同頻率，不同幅度的兩種正弦波信號，并將這兩個信號疊加為一個信號，觀察這三個信號的波形。對疊加后的信號用FFT作譜分析。要求：繪出正弦信號的時域波形掌握傅立葉變換及其逆變換利用傅立葉變換繪出正弦信號的頻譜疊加后的正弦信號經(jīng)傅立葉變換后的頻譜8642432012f(KHz)傅立葉逆變換后得到原信號134|S(f)|(V/KHz)a(t)(V)(ms)四、源程序代碼%% waveform and spectrum of sin signal close all k=10。f1=1。f2=2。N=2^k。dt=。%ms df=1/(N*dt)。% KHz T=N*dt。% 截短時間 Bs=N*df/2。% 系統(tǒng)帶寬 f=[Bs:df:Bsdf]。%頻域橫坐標 t=[T/2:dt:T/2dt]。%時域橫坐標 s1=2*sin(2*pi*f1*t)。s2=3*sin(2*pi*f2*t)。s=s1+s2。[f,S]=T2F(t,s)。% S是s的傅氏變換 [t,a]=F2T(f,S)。% a是S的傅氏反變換 a=real(a)。as=abs(S)。f0=max(f1,f2)。figure(1)subplot(3,1,1)。plot(t,s1)。grid axis([2*f0,+2*f0,min(s1),max(s1)])xlabel(39。t39。)。ylabel(39。s139。)title(39。正弦信號s139。)subplot(3,1,2)。plot(t,s2)。grid axis([2*f0,+2*f0,min(s2),max(s2)])xlabel(39。t39。)。ylabel(39。s239。)title(39。正弦信號s239。)subplot(3,1,3)。plot(t,s)。grid axis([2*f0,+2*f0,min(s),max(s)])xlabel(39。t39。)。ylabel(39。s39。)title(39。疊加后的信號s39。)figure(2)subplot(2,1,1)%輸出的頻譜plot(f,as,39。b39。)。grid axis([2*f0,+2*f0,min(as),max(as)])xlabel(39。f(KHz)39。)。ylabel(39。|S(f)|(V/KHz)39。)title(39。疊加后的正弦信號經(jīng)傅立葉變換后的頻譜39。)subplot(2,1,2)plot(t,a,39。black39。)%輸出信號波形畫圖 grid axis([2/f0,+2/f0,5,5])xlabel(39。t(ms)39。)。ylabel(39。a(t)(V)39。)title(39。傅立葉逆變換后得到原信號39。)