通信系統(tǒng)計(jì)算機(jī)仿真-通信信號(hào)與系統(tǒng)分析課件

第3章通信信號(hào)與系統(tǒng)分析

第3章通信信號(hào)與系統(tǒng)分析

本章內(nèi)容3.1信號(hào)和系統(tǒng)的時(shí)域分析3.2信號(hào)和系統(tǒng)的頻域分析3.3帶通信號(hào)的低通等效3.4隨機(jī)信號(hào)分析本章內(nèi)容3.1信號(hào)和系統(tǒng)的時(shí)域分析3.1.1信號(hào)的時(shí)域分析連續(xù)信號(hào)的處理：MATLAB并不能處理連續(xù)信號(hào)，而是用等時(shí)間間隔點(diǎn)的樣值來近似表示連續(xù)信號(hào)。（取樣時(shí)間間隔足夠?。o論是連續(xù)時(shí)間信號(hào)還是離散時(shí)間信號(hào)，MATLAB都是用一個(gè)數(shù)字序列來表示信號(hào)，這個(gè)數(shù)字序列在MATLAB中叫做向量(vector)。通常的情況下，需要與時(shí)間變量相對(duì)應(yīng)。在MATLAB可視化繪圖中，對(duì)于以t為自變量的連續(xù)信號(hào)，在繪圖時(shí)統(tǒng)一用plot函數(shù)；而對(duì)n為自變量的離散序列，在繪圖時(shí)統(tǒng)一用stem函數(shù)。4.1.1、常用信號(hào)的MATLAB表示3.1.1信號(hào)的時(shí)域分析連續(xù)信號(hào)的處理：MATLAB并不能一、幾種常用信號(hào)的表示1．單位沖擊函數(shù)、序列一、幾種常用信號(hào)的表示1．單位沖擊函數(shù)、序列單位沖激函數(shù)、單位沖激序列的實(shí)現(xiàn)案例t=-5:0.01:5;

y=(t==0);

subplot(1,2,1);

plot(t,y,'r');

n=-5:5;

x=(n==0);

subplot(1,2,2);

stem(n,x);單位沖激函數(shù)、單位沖激序列的實(shí)現(xiàn)案例t=-5:0.01:2．單位階躍函數(shù)、序列

(一）幾種常用信號(hào)的表示2．單位階躍函數(shù)、序列(一）幾種常用信號(hào)的表示單位階躍函數(shù)、單位階躍序列的實(shí)現(xiàn)案例t=-5:0.01:5;y=(t>=0);subplot(1,2,1);plot(t,y,'r')

n=-5:5;

x=(n>=0);

subplot(1,2,2);

stem(n,x);單位階躍函數(shù)、單位階躍序列的實(shí)現(xiàn)案例t=-5:0.01:53．實(shí)指數(shù)序列直接實(shí)現(xiàn)：n=[ns:nf];x=a.^n;4．復(fù)指數(shù)序列直接實(shí)現(xiàn)：n=[ns:nf];x=exp((sigema+jw)*n);5．正(余)弦序列直接實(shí)現(xiàn)：n=[ns:nf];x=cos(w*n+sita);3．實(shí)指數(shù)序列直接實(shí)現(xiàn)：n=[ns:nf];x=a.^n;生成上述三種信號(hào)的實(shí)現(xiàn)案例t=-5:0.01:5;subplot(2,2,1);a=2y1=2.^tplot(t,y1,'r')subplot(2,2,2);a=2;theat=pi/3;y2=sin(2*pi*t+theat)plot(t,y2)subplot(2,1,2);w=4;y3=exp((a+j*w)*t);plot(t,y3,‘y')生成上述三種信號(hào)的實(shí)現(xiàn)案例t=-5:0.01:5;工具箱中的信號(hào)產(chǎn)生函數(shù)函數(shù)名功能函數(shù)名功能sawtooth產(chǎn)生鋸齒波或三角波信號(hào)pulstran產(chǎn)生沖激串square產(chǎn)生方波信號(hào)rectpule產(chǎn)生非周期的方波信號(hào)sinc產(chǎn)生sinc函數(shù)波形tripuls產(chǎn)生非周期的三角波信號(hào)chirp產(chǎn)生調(diào)頻余弦信號(hào)diric產(chǎn)生Dirichlet或周期sinc函數(shù)gauspuls產(chǎn)生高斯正弦脈沖信號(hào)gmonopuls產(chǎn)生高斯單脈沖信號(hào)vco電壓控制振蕩器工具箱中的信號(hào)產(chǎn)生函數(shù)函數(shù)名功能函數(shù)名功能sawtooth產(chǎn)1、sawtooth函數(shù)產(chǎn)生鋸齒波或三角波例：產(chǎn)生f=50Hz的鋸齒波、三角波Fs=10000;%采樣頻率t=0:1/Fs:0.1;%采樣間隔1/Fsf=50;%50Hzx1=sawtooth(2*pi*50*t,0);x2=sawtooth(2*pi*50*t,1);x3=sawtooth(2*pi*50*t,0.5);subplot(311);plot(t,x1);subplot(312);plot(t,x2);subplot(313);plot(t,x3);1、sawtooth函數(shù)2、square函數(shù)產(chǎn)生矩形波例：產(chǎn)生50Hz占空比分別為20％和50％的矩形波。Fs=10000;%采樣頻率t=0:1/Fs:0.1;%采樣間隔1/Fsf=50;%50Hzx1=square(2*pi*50*t,20);x2=square(2*pi*50*t,50);subplot(211);plot(t,x1);subplot(212);plot(t,x2);2、square函數(shù)3、sinc函數(shù)產(chǎn)生sinc波形或sin(πt)/(πt)波形例：t=linspace(-10,+10,200);x=sinc(t);plot(t,x);3、sinc函數(shù)2信號(hào)運(yùn)算1信號(hào)的相加與相乘

y(n)=x1(n)+x2(n)MATLAB實(shí)現(xiàn)：y=x1+x2y(n)=x1(n)×x2(n)MATLAB實(shí)現(xiàn)：y=x1+x2;y=x1.*x2例4-10已知兩序列[0.5,1.5,1,-0.5]與[2,2,2,2,2,2]，試求該兩序列的和與積example4_10.m2信號(hào)運(yùn)算1信號(hào)的相加與相乘例4-10已知兩序列function[f,k]=sigadd(f1,k1,f2,k2)%實(shí)現(xiàn)序列f1,f2的相加,相減,相乘,可據(jù)實(shí)際需要作選擇%f1,k1;f2,k2是參加運(yùn)算的序列向量及其時(shí)間向量%f,k作為返回的和(差,積)序列及其時(shí)間向量%將f1,f2轉(zhuǎn)換成等長序列s1,s2k=min(min(k1),min(k2)):max(max(k1),max(k2));s1=zeros(1,length(k));s2=s1;%初始化序列s1(find((k>=min(k1))&(k<=max(k1))==1))=f1;s2(find((k>=min(k2))&(k<=max(k2))==1))=f2;f=s1+s2;%序列相加%f=s1-s2;%序列相減%f=s1.*s2;%序列相乘%stem(k,f,'fill')；%axis([(min(min(k1),min(k2))-1),(max(max(k1),max(k2))+1),(min(f)-0.5),(max(f)+0.5)])function[f,k]=sigadd(f1,k1,f22序列移位與周期延拓運(yùn)算序列移位：y(n)=x(n-m)。MATLAB實(shí)現(xiàn)：y=x;ny=nx+m序列圓周移位：y(n)=x((n))M，MATLAB實(shí)現(xiàn)：ny=nxs:nxf；y=x(mod(ny,M)+1)3序列翻褶序列翻褶：y(n)=x(-n)。MATLAB可實(shí)現(xiàn):y=fliplr(x)2序列移位與周期延拓運(yùn)算序列移位：y(n)=x(n-m)。3例：實(shí)現(xiàn)f(t+2)，f(-t)f1=f;t1=t-2;subplot(2,2,2);plot(t1,f1)axis([-5502]);f2=fliplr(f);subplot(2,2,3);plot(t,f2)例：實(shí)現(xiàn)f(t+2)，f(-t)f1=f;.4兩序列的卷積運(yùn)算兩序列卷積運(yùn)算：

MATLAB實(shí)現(xiàn)：y=conv(x1,x2)。序列x1(n)和x2(n)必須長度有限。設(shè)x1的非零區(qū)間為n1~n2，長度為L1=n2-n1+1；x2的非零區(qū)間為m1~m2，長度為L2=m2-m1+1；則Y的長度為L=L1+L2-1；Y對(duì)應(yīng)的橫軸序號(hào)

為n1+m1~n1+m1+L-1.4兩序列的卷積運(yùn)算兩序列卷積運(yùn)算：MATLAB實(shí)連續(xù)卷積和離散卷積的關(guān)系：所以，可以用離散卷積和CONV（）求連續(xù)卷積，只需足夠小以及在卷積和的基礎(chǔ)上乘以。連續(xù)卷積和離散卷積的關(guān)系：所以，可以用離散卷積和CONV（）兩序列的相關(guān)運(yùn)算兩序列相關(guān)運(yùn)算：

MATLAB實(shí)現(xiàn)：y=xcorr(x1,x2)。序列累加的數(shù)學(xué)描述為：

MATLAB實(shí)現(xiàn)：y=cumsum(x)序列累加運(yùn)算兩序列的相關(guān)運(yùn)算兩序列相關(guān)運(yùn)算：MATLAB實(shí)現(xiàn)：y=xc信號(hào)的能量和功率1.信號(hào)能量數(shù)字定義：MATLAB實(shí)現(xiàn):E=sum(x.*conj(x));或E=sum(abs(x).^2);數(shù)字定義：2.信號(hào)功率MATLAB實(shí)現(xiàn):P=sum(x.*conj(x))/N;或E=sum(abs(x).^2)/N;信號(hào)的能量和功率1.信號(hào)能量數(shù)字定義：MATLAB實(shí)現(xiàn):E3.1.2線性系統(tǒng)時(shí)域分析

線性時(shí)(移)不變系統(tǒng)特性線性時(shí)(移)不變系統(tǒng)表示方法線性時(shí)(移)不變系統(tǒng)的時(shí)域響應(yīng)3.1.2線性系統(tǒng)時(shí)域分析

線性時(shí)(移)不變系統(tǒng)特性線性時(shí)(移)不變系統(tǒng)特性

假設(shè)系統(tǒng)在輸入信號(hào)x1(t)作用時(shí)的響應(yīng)信號(hào)為y1(t)，在輸入信號(hào)x2(t)作用時(shí)的響應(yīng)信號(hào)為y2(t)，給定兩個(gè)常數(shù)a和b，如果當(dāng)輸入信號(hào)為x(t)時(shí)系統(tǒng)的響應(yīng)信號(hào)為y(t)，且滿足

x(t)=x1(t)+x2(t)（4-15）

y(t)=y1(t)+y2(t)（4-16）則該系統(tǒng)具有疊加性（Additivity）。如果滿足

x(t)=a*x1(t)（4-17）

y(t)=a*y1(t)（4-18）則該系統(tǒng)具有齊次性（Homogeneity）。一個(gè)系統(tǒng)如果是線性系統(tǒng)的話，那么這個(gè)系統(tǒng)必須同時(shí)具有疊加性和齊次性。線性時(shí)(移)不變系統(tǒng)特性

假設(shè)系統(tǒng)在輸入信號(hào)x1(t)作用時(shí)線性時(shí)(移)不變系統(tǒng)表示方法1．常系數(shù)線性微分/差分方程2．系統(tǒng)傳遞函數(shù)：tf3．零－極點(diǎn)增益模型zpk連續(xù)系統(tǒng)：

連續(xù)系統(tǒng)：

離散系統(tǒng)：

離散系統(tǒng)：

線性時(shí)(移)不變系統(tǒng)表示方法1．常系數(shù)線性微分/差分方程24．極點(diǎn)留數(shù)模型離散系統(tǒng)：

連續(xù)系統(tǒng)：

5．二次分式模型sos連續(xù)系統(tǒng)：

離散系統(tǒng)：

6．狀態(tài)空間模型ss連續(xù)系統(tǒng)：

離散系統(tǒng)：

4．極點(diǎn)留數(shù)模型離散系統(tǒng)：連續(xù)系統(tǒng)：5．二次分式模型so系統(tǒng)模型的轉(zhuǎn)換函數(shù)在MATLAB中，用sos、ss、tf、zp分別表示二次分式模型、狀態(tài)空間模型、傳遞函數(shù)模型和零－極點(diǎn)增益模型。其中sos表示二次分式，g為比例系數(shù)，sos為L×6的矩陣，即

（4－15）

1．ss2tf函數(shù)

格式：[num,den]=ss2tf(A,B,C,D,iu)功能：將指定輸入量iu的線性系統(tǒng)(A,B,C,D)轉(zhuǎn)換為傳遞函數(shù)模型[num,den]。2．zp2tf函數(shù)格式：[num,den]=zp2tf(z,p,k)功能：將給定系統(tǒng)的零－極點(diǎn)增益模型轉(zhuǎn)換為傳遞函數(shù)模型，z、p、k分別為零點(diǎn)列向量、極點(diǎn)列向量和增益系數(shù)。系統(tǒng)模型的轉(zhuǎn)換函數(shù)在MATLAB中，用sos線性系統(tǒng)模型的變換函數(shù)函數(shù)名功能說明函數(shù)名功能說明ss2tf狀態(tài)空間模型轉(zhuǎn)換為傳遞函數(shù)模型zp2tf零－極點(diǎn)增益模型轉(zhuǎn)換為傳遞函數(shù)模型ss2zp狀態(tài)空間模型轉(zhuǎn)換為零－極點(diǎn)增益模型zp2ss零－極點(diǎn)增益模型轉(zhuǎn)換為狀態(tài)空間模型ss2sos狀態(tài)空間模型轉(zhuǎn)換為二次分式模型zp2sos零－極點(diǎn)增益模型轉(zhuǎn)換為二次分式模型tf2ss傳遞函數(shù)模型轉(zhuǎn)換為狀態(tài)空間模型sos2tf二次分式模型轉(zhuǎn)換為傳遞函數(shù)模型tf2zp傳遞函數(shù)模型轉(zhuǎn)換為零－極點(diǎn)增益模型sos2zp二次分式模型轉(zhuǎn)換為零－極點(diǎn)增益模型tf2sos傳遞函數(shù)模型轉(zhuǎn)換為二次分式模型sos2ss二次分式模型轉(zhuǎn)換為狀態(tài)空間模型線性系統(tǒng)模型的變換函數(shù)函數(shù)名功能說明函數(shù)名功能說明ss2tf一、連續(xù)系統(tǒng)的時(shí)域響應(yīng)及MATLAB實(shí)現(xiàn)二、離散時(shí)間系統(tǒng)的時(shí)域響應(yīng)及MATLAB實(shí)現(xiàn)線性時(shí)(移)不變系統(tǒng)的時(shí)域響應(yīng)

一、連續(xù)系統(tǒng)的時(shí)域響應(yīng)及MATLAB實(shí)現(xiàn)線性時(shí)(移)不變連續(xù)系統(tǒng)的時(shí)域響應(yīng)及MATLAB實(shí)現(xiàn)1、impulse函數(shù)功能：求連續(xù)系統(tǒng)的單位沖激響應(yīng)。impulse（）函數(shù)有如下幾種調(diào)用格式：impulse(b,a)：該調(diào)用格式以默認(rèn)方式由向量a和b表示的連續(xù)系統(tǒng)在指定時(shí)間范圍內(nèi)的沖激響應(yīng)h(t)的時(shí)域波形圖，并能求出指定時(shí)間范圍內(nèi)沖激響應(yīng)的數(shù)值解。impulse(b,a,t)：繪出系統(tǒng)在0～t時(shí)間范圍內(nèi)沖激響應(yīng)的時(shí)域波形。impulse(b,a,t1:p:t2)：繪出在t1~t2時(shí)間范圍內(nèi)，且以時(shí)間間隔p均勻取樣的沖激響應(yīng)波形。y=impulse(b,a,t1:p:t2)：不繪出波形，而是求出系統(tǒng)沖激響應(yīng)的數(shù)值解：連續(xù)系統(tǒng)的時(shí)域響應(yīng)及MATLAB實(shí)現(xiàn)1、impulse函連續(xù)系統(tǒng)的時(shí)域響應(yīng)2、step函數(shù)功能：求連續(xù)系統(tǒng)的單位階躍響應(yīng)，與impulse（）函數(shù)一樣，也有四種調(diào)用格式。3、initial函數(shù)功能：求連續(xù)系統(tǒng)的零輸入響應(yīng)。4、lsim函數(shù)功能：求任意輸入的連續(xù)系統(tǒng)進(jìn)行仿真?！纠?-20】已知二階系統(tǒng)y″(t)+y(t)=f(t);y′(0+)=-1，y(0+)=0，f(t)=cost,求系統(tǒng)的單位階躍響應(yīng)、零輸入響應(yīng)及全響應(yīng)。

連續(xù)系統(tǒng)的時(shí)域響應(yīng)2、step函數(shù)3、initial函數(shù)【例*1、impz函數(shù)求離散系統(tǒng)（數(shù)字濾波器）的單位沖擊響應(yīng)。離散時(shí)間系統(tǒng)響應(yīng)·impz(b,a)以默認(rèn)方式繪出向量

a和

b定義的離散系統(tǒng)

h(k)的時(shí)域波形；·impz(b,a,n)繪出向量

a和

b定義的離散系統(tǒng)在

0～n（n必須為整數(shù)）離散時(shí)間范圍內(nèi)單位脈沖響應(yīng)

h(k)的時(shí)域波形；·impz(b,a,n1:n2)繪出向量

a和

b定義的離散系統(tǒng)在

n1～n2（n1、n2必須為整數(shù)，且n1<n2）離散時(shí)間范圍內(nèi)單位脈沖響應(yīng)h(k)的時(shí)域波形；·y=impz(b,a,n1:n2)不繪出系統(tǒng)的

h(k)的時(shí)域波形，而是求出向量

a和

b定義的離散系統(tǒng)在

n1～n2（n1、n2必須為整數(shù)，且

n1<n2）離散時(shí)間范圍內(nèi)單位脈沖響應(yīng)

h(k)的數(shù)值解。*1、impz函數(shù)離散時(shí)間系統(tǒng)響應(yīng)·impz(b,a)以默2、stepz函數(shù)功能：求離散系統(tǒng)的單位階躍響應(yīng)格式：stepz(a,b,c,d)、stepz(num,den)3、dinitial函數(shù)功能：求離散系統(tǒng)的零輸入響應(yīng)格式：dinitial(a,b,c,d,x0)、dinitial(num,den,x0)4、dlsim函數(shù)功能：求任意輸入的離散系統(tǒng)進(jìn)行仿真格式：dlsim(a,b,c,d,u)、dlsim(num,den,u)5、fliter函數(shù)功能：求任意輸入的離散系統(tǒng)進(jìn)行仿真格式：filter(a,b,c,d,u)、filter(num,den,u)2、stepz函數(shù)設(shè)系統(tǒng)差分方程為,求該系統(tǒng)對(duì)信號(hào)的響應(yīng)

B=1;A=[1,-0.8];N=0:31;x=0.8.^n;y=filter(B,A,x);subplot(2,1,1);stem(x)subplot(2,1,2);stem(y)例3-22：設(shè)系統(tǒng)差分方程為3.1.3信號(hào)與系統(tǒng)的頻域分析信號(hào)的傅里葉分析連續(xù)時(shí)間系統(tǒng)的S域分析離散時(shí)間系統(tǒng)的Z域分析3.1.3信號(hào)與系統(tǒng)的頻域分析信號(hào)的傅里葉分析信號(hào)的傅里葉分析Matlab的符號(hào)運(yùn)算工具箱（SymbolicMathToolbox）提供了能直接求解傅里葉變換和逆變換的符號(hào)運(yùn)算函數(shù)fourier()和ifourier()。兩函數(shù)的調(diào)用格式如下。（1）傅里葉變換在Matlab中，傅里變換變換由函數(shù)fourier()實(shí)現(xiàn)。fourier()有三種調(diào)用格式：①F=fourier(f)求時(shí)間函數(shù)f(t)的傅里葉變換，返回函數(shù)F的自變量默認(rèn)為w，②F=fourier(f,v)求時(shí)間函數(shù)f(t)的傅里葉變換，返回函數(shù)F的自變量為v，③F=fourier(f,u,v)對(duì)自變量為u的函數(shù)f(u)求傅里葉變換，返回函數(shù)F的自變量為v，1、非周期連續(xù)信號(hào)的頻域分析-傅氏變換信號(hào)的傅里葉分析Matlab的符號(hào)運(yùn)算工具箱（Symboli（2）傅里葉逆變換在Matlab中，傅里變換逆變換由函數(shù)ifourier()實(shí)現(xiàn)。與函數(shù)fourier()相類似，ifourier()也有三種調(diào)用格式：①f=ifourier(F)②f=ifourier(F,u)③f=ifourier(F,v,u)（2）傅里葉逆變換例3-23：求單邊指數(shù)函數(shù)的傅里葉變換，畫出其幅頻特性和相頻特性圖symstwff=exp(-2*t)*sym('Heaviside(t)');F=fourier(f)subplot(3,1,1);ezplot(f,[0:2,0:1.2]);subplot(3,1,2);ezplot(abs(F),[-10:10]);subplot(3,1,3);ezplot(angle(F),[-10:10])例3-23：求單邊指數(shù)函數(shù)用數(shù)值計(jì)算的方法計(jì)算連續(xù)時(shí)間信號(hào)的傅里葉變換需要有個(gè)限定條件，即信號(hào)是時(shí)限信號(hào)（Timelimitedsignal），也就是當(dāng)時(shí)間|t|大于某個(gè)給定時(shí)間時(shí)其值衰減為零或接近于零。采用數(shù)值計(jì)算算法的理論依據(jù)是：若信號(hào)為時(shí)限信號(hào)，當(dāng)時(shí)間間隔T取得足夠小時(shí)，上式可演變?yōu)椋河肕ATLAB表示為：

X=x*exp(-j*w*t’)*T用數(shù)值計(jì)算的方法計(jì)算連續(xù)時(shí)間信號(hào)的傅里葉變換需要有個(gè)限定條件例3-25已知非周期連續(xù)信號(hào)求它的傅里葉變換example4_25.m例3-25已知非周期連續(xù)信號(hào)2、周期連續(xù)信號(hào)的頻域分析-傅里葉級(jí)數(shù)周期T1基頻0=2/T1假設(shè)需要計(jì)算的諧波次數(shù)為N，則總的系數(shù)個(gè)數(shù)為2N+1，對(duì)于某個(gè)系數(shù)：

傅里葉級(jí)數(shù)的系數(shù)為2、周期連續(xù)信號(hào)的頻域分析-傅里葉級(jí)數(shù)周期T1假設(shè)需要計(jì)算的例4-26：給定一個(gè)周期為T1=2s的連續(xù)時(shí)間周期方波信號(hào)，其一個(gè)周期內(nèi)的數(shù)學(xué)表達(dá)式為

試求其傅里葉級(jí)數(shù)。例4-26：給定一個(gè)周期為T1=2s的連續(xù)時(shí)間周期方波信3時(shí)間離散、連續(xù)頻率－序列傅里葉變換4離散時(shí)間、離散頻率－離散傅里葉級(jí)數(shù)5離散時(shí)間、離散頻率－離散傅里葉變換(DFT)正變換：

逆變換：

正變換：

逆變換：

正變換：

逆變換：

3時(shí)間離散、連續(xù)頻率－序列傅里葉變換4離散時(shí)間、離散頻率1．一維快速正傅里葉變換函數(shù)fft格式：X=fft(x,N)功能：采用FFT算法計(jì)算序列向量x的N點(diǎn)DFT變換，當(dāng)N缺省時(shí)，fft函數(shù)自動(dòng)按x的長度計(jì)算DFT。當(dāng)N為2整數(shù)次冪時(shí)，fft按基-2算法計(jì)算，否則用混合算法。2．一維快速逆傅里葉變換函數(shù)ifft格式：x=ifft(X,N)功能：采用FFT算法計(jì)算序列向量X的N點(diǎn)IDFT變換。[例4-36]用快速傅里葉變換FFT計(jì)算下面兩個(gè)序列的卷積。,并測(cè)試直接卷積和快速卷積的時(shí)間。圖4.35快速卷積框圖6、快速傅里葉變換1．一維快速正傅里葉變換函數(shù)fft[例4-36]用快速傅里MATLAB程序(部分)：%線性卷積xn=sin(0.4*[1:15]);%對(duì)序列x(n)賦值,M=15hn=0.9.^(1:20);%對(duì)序列h(n)賦值,N=20yn=conv(xn,hn); %直接調(diào)用函數(shù)conv計(jì)算卷積%圓周卷積L=pow2(nextpow2(M+N-1));

Xk=fft(xn,L);

Hk=fft(hn,L);

Yk=Xk.*Hk;

yn=ifft(Yk,L);

圖4.36x(n),h(n)及其線性卷積波形MATLAB程序(部分)：圖4.36x(n),h(n)及連續(xù)時(shí)間系統(tǒng)的S域分析Matlab的符號(hào)數(shù)學(xué)工具箱（SymbolicMathToolbox）提供了能直接求解拉普拉斯變換和逆變換的符號(hào)運(yùn)算函數(shù)laplace()和ilaplace()。1.連續(xù)時(shí)間信號(hào)的拉普拉斯變換連續(xù)時(shí)間系統(tǒng)的S域分析Matlab的符號(hào)數(shù)學(xué)工具箱（Symb例3-30：求以下函數(shù)的拉普拉斯變換。symstf1=sym('exp(-2*t)*Heaviside(t)');F1=laplace(f1)%求f1(t)的拉普拉斯變換f2=sym('t*exp(-t)*Heaviside(t)');F2=laplace(f2)例3-30：求以下函數(shù)的拉普拉斯變換。symstMatlab工具箱中提供了專門用于分析連續(xù)系統(tǒng)頻率響應(yīng)的freqs()函數(shù)。它可求出系統(tǒng)頻率響應(yīng)的數(shù)值解，也可繪出其幅頻及相頻特性曲線。freqs()函數(shù)的調(diào)用形式有如下四種。（1）H＝freqs(b,a,w1:dw:w2)可求得指定頻率范圍（w1~w2）內(nèi)相應(yīng)頻點(diǎn)處系統(tǒng)頻率響應(yīng)的樣值。其中，w1、w2分別為頻率起始值和終止值，dw為頻率取樣間隔。（2）[H,w]＝freqs(b,a)將計(jì)算默認(rèn)頻率范圍內(nèi)200個(gè)頻點(diǎn)上系統(tǒng)頻率響應(yīng)的樣值，并賦值給返回向量H，200個(gè)頻點(diǎn)則記錄在向量w中。（3）[H,w]＝freqs(b,a,n)將計(jì)算默認(rèn)頻率范圍內(nèi)n個(gè)頻點(diǎn)上系統(tǒng)頻率響應(yīng)的樣值，并賦值給返回向量H，n個(gè)頻點(diǎn)則記錄在向量w中。（4）freqs(b,a)該調(diào)用方式將繪出系統(tǒng)的幅頻特性和相頻特性曲線。2.連續(xù)系統(tǒng)的復(fù)頻域分析

Matlab工具箱中提供了專門用于分析連續(xù)系統(tǒng)頻率響應(yīng)的fr例3-32：已知一因果系統(tǒng)的系統(tǒng)函數(shù)為

試分析其頻率特性。a=[12321];b=[10-4];p=roots(a)pxm=max(real(p));ifpxm>=0'系統(tǒng)不穩(wěn)定'elsefreqs(b,a)end例3-32：已知一因果系統(tǒng)的系統(tǒng)函數(shù)為

試分析其頻率特性。離散時(shí)間系統(tǒng)的Z域分析1.離散時(shí)間信號(hào)的Z變換MATLAB的符號(hào)數(shù)學(xué)工具箱(SymbolicMathTools)提供了計(jì)算z正變換的函數(shù)ztrans和計(jì)算逆z變換的函數(shù)iztrans。其調(diào)用形式為：F=ztrans(f)%求符號(hào)函數(shù)f的z變換，返回函數(shù)的自變量為z；F=ztrans(f,w)%求符號(hào)函數(shù)f的z變換，返回函數(shù)的自變量為w；F=ztrans(f,k,w)%對(duì)自變量為k的符號(hào)函數(shù)f求z變換，返回函數(shù)的自變量為w。f=iztrans(F)%對(duì)自變量為z的符號(hào)函數(shù)F求逆z變換，返回函數(shù)的自變量為n；f=iztrans(F,k)%對(duì)自變量為z的符號(hào)函數(shù)F求逆z變換，返回函數(shù)的自變量為k；f=iztrans(F,w,k)%對(duì)自變量為w的符號(hào)函數(shù)F求逆z變換，返回函數(shù)的自變量為k。離散時(shí)間系統(tǒng)的Z域分析1.離散時(shí)間信號(hào)的Z變換2、離散時(shí)間系統(tǒng)的Z域分析利用freqz()函數(shù)可求得系統(tǒng)的頻率響應(yīng)。調(diào)用格式如下：freqz(b，a)該調(diào)用方式將繪制系統(tǒng)在0~范圍內(nèi)的幅頻特性和相頻特性圖，其中，b、a分別為系統(tǒng)函數(shù)分子、分母多項(xiàng)式的系數(shù)向量。freqz(b，a，’whole’)該調(diào)用方式將繪制系統(tǒng)在0~2范圍內(nèi)的幅頻特性和相頻特性圖。freqz(b，a，N)該調(diào)用方式將繪制系統(tǒng)在0~范圍內(nèi)N個(gè)頻率等分點(diǎn)的幅頻特性和相頻特性圖，N的缺省值為512；freqz(b，a，N，’whole’)該調(diào)用方式將繪制系統(tǒng)在0~2范圍內(nèi)N個(gè)頻率等分點(diǎn)的幅頻特性和相頻特性圖。2、離散時(shí)間系統(tǒng)的Z域分析a.畫出零極點(diǎn)圖

b.求系統(tǒng)響應(yīng)

c.求系統(tǒng)的幅頻特性和相頻特性例4-37：設(shè)a.畫出零極點(diǎn)圖

b.求系統(tǒng)響應(yīng)

c.通信系統(tǒng)計(jì)算機(jī)仿真-通信信號(hào)與系統(tǒng)分析課件第3章通信信號(hào)與系統(tǒng)分析

第3章通信信號(hào)與系統(tǒng)分析

本章內(nèi)容3.1信號(hào)和系統(tǒng)的時(shí)域分析3.2信號(hào)和系統(tǒng)的頻域分析3.3帶通信號(hào)的低通等效3.4隨機(jī)信號(hào)分析本章內(nèi)容3.1信號(hào)和系統(tǒng)的時(shí)域分析3.1.1信號(hào)的時(shí)域分析連續(xù)信號(hào)的處理：MATLAB并不能處理連續(xù)信號(hào)，而是用等時(shí)間間隔點(diǎn)的樣值來近似表示連續(xù)信號(hào)。（取樣時(shí)間間隔足夠小）無論是連續(xù)時(shí)間信號(hào)還是離散時(shí)間信號(hào)，MATLAB都是用一個(gè)數(shù)字序列來表示信號(hào)，這個(gè)數(shù)字序列在MATLAB中叫做向量(vector)。通常的情況下，需要與時(shí)間變量相對(duì)應(yīng)。在MATLAB可視化繪圖中，對(duì)于以t為自變量的連續(xù)信號(hào)，在繪圖時(shí)統(tǒng)一用plot函數(shù)；而對(duì)n為自變量的離散序列，在繪圖時(shí)統(tǒng)一用stem函數(shù)。4.1.1、常用信號(hào)的MATLAB表示3.1.1信號(hào)的時(shí)域分析連續(xù)信號(hào)的處理：MATLAB并不能一、幾種常用信號(hào)的表示1．單位沖擊函數(shù)、序列一、幾種常用信號(hào)的表示1．單位沖擊函數(shù)、序列單位沖激函數(shù)、單位沖激序列的實(shí)現(xiàn)案例t=-5:0.01:5;

y=(t==0);

subplot(1,2,1);

plot(t,y,'r');

n=-5:5;

x=(n==0);

subplot(1,2,2);

stem(n,x);單位沖激函數(shù)、單位沖激序列的實(shí)現(xiàn)案例t=-5:0.01:2．單位階躍函數(shù)、序列

(一）幾種常用信號(hào)的表示2．單位階躍函數(shù)、序列(一）幾種常用信號(hào)的表示單位階躍函數(shù)、單位階躍序列的實(shí)現(xiàn)案例t=-5:0.01:5;y=(t>=0);subplot(1,2,1);plot(t,y,'r')

n=-5:5;

x=(n>=0);

subplot(1,2,2);

stem(n,x);單位階躍函數(shù)、單位階躍序列的實(shí)現(xiàn)案例t=-5:0.01:53．實(shí)指數(shù)序列直接實(shí)現(xiàn)：n=[ns:nf];x=a.^n;4．復(fù)指數(shù)序列直接實(shí)現(xiàn)：n=[ns:nf];x=exp((sigema+jw)*n);5．正(余)弦序列直接實(shí)現(xiàn)：n=[ns:nf];x=cos(w*n+sita);3．實(shí)指數(shù)序列直接實(shí)現(xiàn)：n=[ns:nf];x=a.^n;生成上述三種信號(hào)的實(shí)現(xiàn)案例t=-5:0.01:5;subplot(2,2,1);a=2y1=2.^tplot(t,y1,'r')subplot(2,2,2);a=2;theat=pi/3;y2=sin(2*pi*t+theat)plot(t,y2)subplot(2,1,2);w=4;y3=exp((a+j*w)*t);plot(t,y3,‘y')生成上述三種信號(hào)的實(shí)現(xiàn)案例t=-5:0.01:5;工具箱中的信號(hào)產(chǎn)生函數(shù)函數(shù)名功能函數(shù)名功能sawtooth產(chǎn)生鋸齒波或三角波信號(hào)pulstran產(chǎn)生沖激串square產(chǎn)生方波信號(hào)rectpule產(chǎn)生非周期的方波信號(hào)sinc產(chǎn)生sinc函數(shù)波形tripuls產(chǎn)生非周期的三角波信號(hào)chirp產(chǎn)生調(diào)頻余弦信號(hào)diric產(chǎn)生Dirichlet或周期sinc函數(shù)gauspuls產(chǎn)生高斯正弦脈沖信號(hào)gmonopuls產(chǎn)生高斯單脈沖信號(hào)vco電壓控制振蕩器工具箱中的信號(hào)產(chǎn)生函數(shù)函數(shù)名功能函數(shù)名功能sawtooth產(chǎn)1、sawtooth函數(shù)產(chǎn)生鋸齒波或三角波例：產(chǎn)生f=50Hz的鋸齒波、三角波Fs=10000;%采樣頻率t=0:1/Fs:0.1;%采樣間隔1/Fsf=50;%50Hzx1=sawtooth(2*pi*50*t,0);x2=sawtooth(2*pi*50*t,1);x3=sawtooth(2*pi*50*t,0.5);subplot(311);plot(t,x1);subplot(312);plot(t,x2);subplot(313);plot(t,x3);1、sawtooth函數(shù)2、square函數(shù)產(chǎn)生矩形波例：產(chǎn)生50Hz占空比分別為20％和50％的矩形波。Fs=10000;%采樣頻率t=0:1/Fs:0.1;%采樣間隔1/Fsf=50;%50Hzx1=square(2*pi*50*t,20);x2=square(2*pi*50*t,50);subplot(211);plot(t,x1);subplot(212);plot(t,x2);2、square函數(shù)3、sinc函數(shù)產(chǎn)生sinc波形或sin(πt)/(πt)波形例：t=linspace(-10,+10,200);x=sinc(t);plot(t,x);3、sinc函數(shù)2信號(hào)運(yùn)算1信號(hào)的相加與相乘

y(n)=x1(n)+x2(n)MATLAB實(shí)現(xiàn)：y=x1+x2y(n)=x1(n)×x2(n)MATLAB實(shí)現(xiàn)：y=x1+x2;y=x1.*x2例4-10已知兩序列[0.5,1.5,1,-0.5]與[2,2,2,2,2,2]，試求該兩序列的和與積example4_10.m2信號(hào)運(yùn)算1信號(hào)的相加與相乘例4-10已知兩序列function[f,k]=sigadd(f1,k1,f2,k2)%實(shí)現(xiàn)序列f1,f2的相加,相減,相乘,可據(jù)實(shí)際需要作選擇%f1,k1;f2,k2是參加運(yùn)算的序列向量及其時(shí)間向量%f,k作為返回的和(差,積)序列及其時(shí)間向量%將f1,f2轉(zhuǎn)換成等長序列s1,s2k=min(min(k1),min(k2)):max(max(k1),max(k2));s1=zeros(1,length(k));s2=s1;%初始化序列s1(find((k>=min(k1))&(k<=max(k1))==1))=f1;s2(find((k>=min(k2))&(k<=max(k2))==1))=f2;f=s1+s2;%序列相加%f=s1-s2;%序列相減%f=s1.*s2;%序列相乘%stem(k,f,'fill')；%axis([(min(min(k1),min(k2))-1),(max(max(k1),max(k2))+1),(min(f)-0.5),(max(f)+0.5)])function[f,k]=sigadd(f1,k1,f22序列移位與周期延拓運(yùn)算序列移位：y(n)=x(n-m)。MATLAB實(shí)現(xiàn)：y=x;ny=nx+m序列圓周移位：y(n)=x((n))M，MATLAB實(shí)現(xiàn)：ny=nxs:nxf；y=x(mod(ny,M)+1)3序列翻褶序列翻褶：y(n)=x(-n)。MATLAB可實(shí)現(xiàn):y=fliplr(x)2序列移位與周期延拓運(yùn)算序列移位：y(n)=x(n-m)。3例：實(shí)現(xiàn)f(t+2)，f(-t)f1=f;t1=t-2;subplot(2,2,2);plot(t1,f1)axis([-5502]);f2=fliplr(f);subplot(2,2,3);plot(t,f2)例：實(shí)現(xiàn)f(t+2)，f(-t)f1=f;.4兩序列的卷積運(yùn)算兩序列卷積運(yùn)算：

MATLAB實(shí)現(xiàn)：y=conv(x1,x2)。序列x1(n)和x2(n)必須長度有限。設(shè)x1的非零區(qū)間為n1~n2，長度為L1=n2-n1+1；x2的非零區(qū)間為m1~m2，長度為L2=m2-m1+1；則Y的長度為L=L1+L2-1；Y對(duì)應(yīng)的橫軸序號(hào)

為n1+m1~n1+m1+L-1.4兩序列的卷積運(yùn)算兩序列卷積運(yùn)算：MATLAB實(shí)連續(xù)卷積和離散卷積的關(guān)系：所以，可以用離散卷積和CONV（）求連續(xù)卷積，只需足夠小以及在卷積和的基礎(chǔ)上乘以。連續(xù)卷積和離散卷積的關(guān)系：所以，可以用離散卷積和CONV（）兩序列的相關(guān)運(yùn)算兩序列相關(guān)運(yùn)算：

MATLAB實(shí)現(xiàn)：y=xcorr(x1,x2)。序列累加的數(shù)學(xué)描述為：

MATLAB實(shí)現(xiàn)：y=cumsum(x)序列累加運(yùn)算兩序列的相關(guān)運(yùn)算兩序列相關(guān)運(yùn)算：MATLAB實(shí)現(xiàn)：y=xc信號(hào)的能量和功率1.信號(hào)能量數(shù)字定義：MATLAB實(shí)現(xiàn):E=sum(x.*conj(x));或E=sum(abs(x).^2);數(shù)字定義：2.信號(hào)功率MATLAB實(shí)現(xiàn):P=sum(x.*conj(x))/N;或E=sum(abs(x).^2)/N;信號(hào)的能量和功率1.信號(hào)能量數(shù)字定義：MATLAB實(shí)現(xiàn):E3.1.2線性系統(tǒng)時(shí)域分析

線性時(shí)(移)不變系統(tǒng)特性線性時(shí)(移)不變系統(tǒng)表示方法線性時(shí)(移)不變系統(tǒng)的時(shí)域響應(yīng)3.1.2線性系統(tǒng)時(shí)域分析

線性時(shí)(移)不變系統(tǒng)特性線性時(shí)(移)不變系統(tǒng)特性

假設(shè)系統(tǒng)在輸入信號(hào)x1(t)作用時(shí)的響應(yīng)信號(hào)為y1(t)，在輸入信號(hào)x2(t)作用時(shí)的響應(yīng)信號(hào)為y2(t)，給定兩個(gè)常數(shù)a和b，如果當(dāng)輸入信號(hào)為x(t)時(shí)系統(tǒng)的響應(yīng)信號(hào)為y(t)，且滿足

x(t)=x1(t)+x2(t)（4-15）

y(t)=y1(t)+y2(t)（4-16）則該系統(tǒng)具有疊加性（Additivity）。如果滿足

x(t)=a*x1(t)（4-17）

y(t)=a*y1(t)（4-18）則該系統(tǒng)具有齊次性（Homogeneity）。一個(gè)系統(tǒng)如果是線性系統(tǒng)的話，那么這個(gè)系統(tǒng)必須同時(shí)具有疊加性和齊次性。線性時(shí)(移)不變系統(tǒng)特性

假設(shè)系統(tǒng)在輸入信號(hào)x1(t)作用時(shí)線性時(shí)(移)不變系統(tǒng)表示方法1．常系數(shù)線性微分/差分方程2．系統(tǒng)傳遞函數(shù)：tf3．零－極點(diǎn)增益模型zpk連續(xù)系統(tǒng)：

連續(xù)系統(tǒng)：

離散系統(tǒng)：

離散系統(tǒng)：

線性時(shí)(移)不變系統(tǒng)表示方法1．常系數(shù)線性微分/差分方程24．極點(diǎn)留數(shù)模型離散系統(tǒng)：

連續(xù)系統(tǒng)：

5．二次分式模型sos連續(xù)系統(tǒng)：

離散系統(tǒng)：

6．狀態(tài)空間模型ss連續(xù)系統(tǒng)：

離散系統(tǒng)：

4．極點(diǎn)留數(shù)模型離散系統(tǒng)：連續(xù)系統(tǒng)：5．二次分式模型so系統(tǒng)模型的轉(zhuǎn)換函數(shù)在MATLAB中，用sos、ss、tf、zp分別表示二次分式模型、狀態(tài)空間模型、傳遞函數(shù)模型和零－極點(diǎn)增益模型。其中sos表示二次分式，g為比例系數(shù)，sos為L×6的矩陣，即

（4－15）

1．ss2tf函數(shù)

格式：[num,den]=ss2tf(A,B,C,D,iu)功能：將指定輸入量iu的線性系統(tǒng)(A,B,C,D)轉(zhuǎn)換為傳遞函數(shù)模型[num,den]。2．zp2tf函數(shù)格式：[num,den]=zp2tf(z,p,k)功能：將給定系統(tǒng)的零－極點(diǎn)增益模型轉(zhuǎn)換為傳遞函數(shù)模型，z、p、k分別為零點(diǎn)列向量、極點(diǎn)列向量和增益系數(shù)。系統(tǒng)模型的轉(zhuǎn)換函數(shù)在MATLAB中，用sos線性系統(tǒng)模型的變換函數(shù)函數(shù)名功能說明函數(shù)名功能說明ss2tf狀態(tài)空間模型轉(zhuǎn)換為傳遞函數(shù)模型zp2tf零－極點(diǎn)增益模型轉(zhuǎn)換為傳遞函數(shù)模型ss2zp狀態(tài)空間模型轉(zhuǎn)換為零－極點(diǎn)增益模型zp2ss零－極點(diǎn)增益模型轉(zhuǎn)換為狀態(tài)空間模型ss2sos狀態(tài)空間模型轉(zhuǎn)換為二次分式模型zp2sos零－極點(diǎn)增益模型轉(zhuǎn)換為二次分式模型tf2ss傳遞函數(shù)模型轉(zhuǎn)換為狀態(tài)空間模型sos2tf二次分式模型轉(zhuǎn)換為傳遞函數(shù)模型tf2zp傳遞函數(shù)模型轉(zhuǎn)換為零－極點(diǎn)增益模型sos2zp二次分式模型轉(zhuǎn)換為零－極點(diǎn)增益模型tf2sos傳遞函數(shù)模型轉(zhuǎn)換為二次分式模型sos2ss二次分式模型轉(zhuǎn)換為狀態(tài)空間模型線性系統(tǒng)模型的變換函數(shù)函數(shù)名功能說明函數(shù)名功能說明ss2tf一、連續(xù)系統(tǒng)的時(shí)域響應(yīng)及MATLAB實(shí)現(xiàn)二、離散時(shí)間系統(tǒng)的時(shí)域響應(yīng)及MATLAB實(shí)現(xiàn)線性時(shí)(移)不變系統(tǒng)的時(shí)域響應(yīng)

一、連續(xù)系統(tǒng)的時(shí)域響應(yīng)及MATLAB實(shí)現(xiàn)線性時(shí)(移)不變連續(xù)系統(tǒng)的時(shí)域響應(yīng)及MATLAB實(shí)現(xiàn)1、impulse函數(shù)功能：求連續(xù)系統(tǒng)的單位沖激響應(yīng)。impulse（）函數(shù)有如下幾種調(diào)用格式：impulse(b,a)：該調(diào)用格式以默認(rèn)方式由向量a和b表示的連續(xù)系統(tǒng)在指定時(shí)間范圍內(nèi)的沖激響應(yīng)h(t)的時(shí)域波形圖，并能求出指定時(shí)間范圍內(nèi)沖激響應(yīng)的數(shù)值解。impulse(b,a,t)：繪出系統(tǒng)在0～t時(shí)間范圍內(nèi)沖激響應(yīng)的時(shí)域波形。impulse(b,a,t1:p:t2)：繪出在t1~t2時(shí)間范圍內(nèi)，且以時(shí)間間隔p均勻取樣的沖激響應(yīng)波形。y=impulse(b,a,t1:p:t2)：不繪出波形，而是求出系統(tǒng)沖激響應(yīng)的數(shù)值解：連續(xù)系統(tǒng)的時(shí)域響應(yīng)及MATLAB實(shí)現(xiàn)1、impulse函連續(xù)系統(tǒng)的時(shí)域響應(yīng)2、step函數(shù)功能：求連續(xù)系統(tǒng)的單位階躍響應(yīng)，與impulse（）函數(shù)一樣，也有四種調(diào)用格式。3、initial函數(shù)功能：求連續(xù)系統(tǒng)的零輸入響應(yīng)。4、lsim函數(shù)功能：求任意輸入的連續(xù)系統(tǒng)進(jìn)行仿真。【例4-20】已知二階系統(tǒng)y″(t)+y(t)=f(t);y′(0+)=-1，y(0+)=0，f(t)=cost,求系統(tǒng)的單位階躍響應(yīng)、零輸入響應(yīng)及全響應(yīng)。

連續(xù)系統(tǒng)的時(shí)域響應(yīng)2、step函數(shù)3、initial函數(shù)【例*1、impz函數(shù)求離散系統(tǒng)（數(shù)字濾波器）的單位沖擊響應(yīng)。離散時(shí)間系統(tǒng)響應(yīng)·impz(b,a)以默認(rèn)方式繪出向量

a和

b定義的離散系統(tǒng)

h(k)的時(shí)域波形；·impz(b,a,n)繪出向量

a和

b定義的離散系統(tǒng)在

0～n（n必須為整數(shù)）離散時(shí)間范圍內(nèi)單位脈沖響應(yīng)

h(k)的時(shí)域波形；·impz(b,a,n1:n2)繪出向量

a和

b定義的離散系統(tǒng)在

n1～n2（n1、n2必須為整數(shù)，且n1<n2）離散時(shí)間范圍內(nèi)單位脈沖響應(yīng)h(k)的時(shí)域波形；·y=impz(b,a,n1:n2)不繪出系統(tǒng)的

h(k)的時(shí)域波形，而是求出向量

a和

b定義的離散系統(tǒng)在

n1～n2（n1、n2必須為整數(shù)，且

n1<n2）離散時(shí)間范圍內(nèi)單位脈沖響應(yīng)

h(k)的數(shù)值解。*1、impz函數(shù)離散時(shí)間系統(tǒng)響應(yīng)·impz(b,a)以默2、stepz函數(shù)功能：求離散系統(tǒng)的單位階躍響應(yīng)格式：stepz(a,b,c,d)、stepz(num,den)3、dinitial函數(shù)功能：求離散系統(tǒng)的零輸入響應(yīng)格式：dinitial(a,b,c,d,x0)、dinitial(num,den,x0)4、dlsim函數(shù)功能：求任意輸入的離散系統(tǒng)進(jìn)行仿真格式：dlsim(a,b,c,d,u)、dlsim(num,den,u)5、fliter函數(shù)功能：求任意輸入的離散系統(tǒng)進(jìn)行仿真格式：filter(a,b,c,d,u)、filter(num,den,u)2、stepz函數(shù)設(shè)系統(tǒng)差分方程為,求該系統(tǒng)對(duì)信號(hào)的響應(yīng)

B=1;A=[1,-0.8];N=0:31;x=0.8.^n;y=filter(B,A,x);subplot(2,1,1);stem(x)subplot(2,1,2);stem(y)例3-22：設(shè)系統(tǒng)差分方程為3.1.3信號(hào)與系統(tǒng)的頻域分析信號(hào)的傅里葉分析連續(xù)時(shí)間系統(tǒng)的S域分析離散時(shí)間系統(tǒng)的Z域分析3.1.3信號(hào)與系統(tǒng)的頻域分析信號(hào)的傅里葉分析信號(hào)的傅里葉分析Matlab的符號(hào)運(yùn)算工具箱（SymbolicMathToolbox）提供了能直接求解傅里葉變換和逆變換的符號(hào)運(yùn)算函數(shù)fourier()和ifourier()。兩函數(shù)的調(diào)用格式如下。（1）傅里葉變換在Matlab中，傅里變換變換由函數(shù)fourier()實(shí)現(xiàn)。fourier()有三種調(diào)用格式：①F=fourier(f)求時(shí)間函數(shù)f(t)的傅里葉變換，返回函數(shù)F的自變量默認(rèn)為w，②F=fourier(f,v)求時(shí)間函數(shù)f(t)的傅里葉變換，返回函數(shù)F的自變量為v，③F=fourier(f,u,v)對(duì)自變量為u的函數(shù)f(u)求傅里葉變換，返回函數(shù)F的自變量為v，1、非周期連續(xù)信號(hào)的頻域分析-傅氏變換信號(hào)的傅里葉分析Matlab的符號(hào)運(yùn)算工具箱（Symboli（2）傅里葉逆變換在Matlab中，傅里變換逆變換由函數(shù)ifourier()實(shí)現(xiàn)。與函數(shù)fourier()相類似，ifourier()也有三種調(diào)用格式：①f=ifourier(F)②f=ifourier(F,u)③f=ifourier(F,v,u)（2）傅里葉逆變換例3-23：求單邊指數(shù)函數(shù)的傅里葉變換，畫出其幅頻特性和相頻特性圖symstwff=exp(-2*t)*sym('Heaviside(t)');F=fourier(f)subplot(3,1,1);ezplot(f,[0:2,0:1.2]);subplot(3,1,2);ezplot(abs(F),[-10:10]);subplot(3,1,3);ezplot(angle(F),[-10:10])例3-23：求單邊指數(shù)函數(shù)用數(shù)值計(jì)算的方法計(jì)算連續(xù)時(shí)間信號(hào)的傅里葉變換需要有個(gè)限定條件，即信號(hào)是時(shí)限信號(hào)（Timelimitedsignal），也就是當(dāng)時(shí)間|t|大于某個(gè)給定時(shí)間時(shí)其值衰減為零或接近于零。采用數(shù)值計(jì)算算法的理論依據(jù)是：若信號(hào)為時(shí)限信號(hào)，當(dāng)時(shí)間間隔T取得足夠小時(shí)，上式可演變?yōu)椋河肕ATLAB表示為：

X=x*exp(-j*w*t’)*T用數(shù)值計(jì)算的方法計(jì)算連續(xù)時(shí)間信號(hào)的傅里葉變換需要有個(gè)限定條件例3-25已知非周期連續(xù)信號(hào)求它的傅里葉變換example4_25.m例3-25已知非周期連續(xù)信號(hào)2、周期連續(xù)信號(hào)的頻域分析-傅里葉級(jí)數(shù)周期T1基頻0=2/T1假設(shè)需要計(jì)算的諧波次數(shù)為N，則總的系數(shù)個(gè)數(shù)為2N+1，對(duì)于某個(gè)系數(shù)：

傅里葉級(jí)數(shù)的系數(shù)為2、周期連續(xù)信號(hào)的頻域分析-傅里葉級(jí)數(shù)周期T1假設(shè)需要計(jì)算的例4-26：給定一個(gè)周期為T1=2s的連續(xù)時(shí)間周期方波信號(hào)，其一個(gè)周期內(nèi)的數(shù)學(xué)表達(dá)式為

試求其傅里葉級(jí)數(shù)。例4-26：給定一個(gè)周期為T1=2s的連續(xù)時(shí)間周期方波信3時(shí)間離散、連續(xù)頻率－序列傅里葉變換4離散時(shí)間、離散頻率－離散傅里葉級(jí)數(shù)5離散時(shí)間、離散頻率－離散傅里葉變換(DFT)正變換：

逆變換：

正變換：

逆變換：

正變換：

逆變換：

3時(shí)間離散、連續(xù)頻率－序列傅里葉變換4離散時(shí)間、離散頻率1．一維快速正傅里葉變換函數(shù)fft格式：X=fft(x,N)功能：采用FFT算法計(jì)算序列向量x的N點(diǎn)DFT變換，當(dāng)N缺省時(shí)，fft函數(shù)自動(dòng)按x的長度計(jì)算DFT。當(dāng)N為2整數(shù)次冪時(shí)，fft按基-2算法計(jì)算，否則用混合算法。2．一維快速逆傅里葉變換函數(shù)ifft格式：x