信號(hào)與系統(tǒng)Matlab實(shí)驗(yàn)報(bào)告

1、實(shí)驗(yàn)一matlab 程序入門(mén)和基礎(chǔ)應(yīng)用一、實(shí)驗(yàn)名稱 matlab 程序入門(mén)和基礎(chǔ)應(yīng)用二、實(shí)驗(yàn)?zāi)康?.學(xué)習(xí)matlab軟件的基本使用方法；2.了解matlab的數(shù)值計(jì)算，符號(hào)運(yùn)算，可視化功能；3. matlab程序設(shè)計(jì)入門(mén)四、實(shí)驗(yàn)設(shè)備計(jì)算機(jī) matlab軟件六、實(shí)驗(yàn)內(nèi)容及具體步驟1、打開(kāi)matlab的系統(tǒng)界面，對(duì)其功能做一個(gè)大致了解；2、學(xué)習(xí)變量的描述方法，掌握幾個(gè)固定變量：i，j，pi，inf的使用。注意，變量描述以字母開(kāi)頭，可以由字母、數(shù)字和下劃線混合組成，區(qū)分字母大，小寫(xiě)字符長(zhǎng)度不超過(guò)31個(gè)。3、學(xué)習(xí)數(shù)值，矩陣，運(yùn)算符，向量的矩陣運(yùn)算，數(shù)組運(yùn)算的描述方法。（1）用一個(gè)簡(jiǎn)單命令求解線性系統(tǒng)3

2、x1+ x2 - x3 =3.6 x1+2x2+4x3 = 2.1 -x1+4x2+5x3 = -1.4a=3 1 -1;1 2 4;-1 4 5;b=3.6;2.1;-1.4;x=ab結(jié)果：x = 1.4818 -0.4606 0.3848（2）用簡(jiǎn)短命令計(jì)算并繪制在0x6范圍內(nèi)的sin(2x)、sinx2、sin2x。x=linspace(0,6)y1=sin(2*x),y2=sin(x.2),y3=(sin(x).2;plot(x,y1,x, y2,x, y3)4、matlab符號(hào)運(yùn)算功能（1）符號(hào)運(yùn)算的過(guò)程在符號(hào)運(yùn)算的整個(gè)過(guò)程中，所有的運(yùn)算均是以符號(hào)進(jìn)行的，即使以數(shù)字形式出現(xiàn)的量也是字

3、符量。做一個(gè)對(duì)sin(x/2)求導(dǎo)的過(guò)程。在命令窗口中輸入如下符號(hào)表達(dá)式按回車： f=sin(x/2);dfdx=diff(f)顯示結(jié)果如下： dfdx = 1/2*cos(1/2*x)整個(gè)求導(dǎo)的過(guò)程都是由符號(hào)變量和符號(hào)表達(dá)式完成，沒(méi)有涉及到具體的數(shù)值運(yùn)算，其中1/2也被當(dāng)作是字符量 。注意：符號(hào)變量前先要進(jìn)行定義，定義語(yǔ)句是：sym 或syms 變量名列表。前者定義一個(gè)單一的符號(hào)變量，后者可以一次定義多個(gè)符號(hào)變量。如： sym a 定義a為符號(hào)變量 sym a b c 定義a,b,c均為符號(hào)變量符號(hào)表達(dá)式是由符號(hào)變量組成的一個(gè)表達(dá)式，符號(hào)方程是將一個(gè)符號(hào)表達(dá)式通過(guò)等號(hào)給一個(gè)符號(hào)變量。凡是用到

4、sym命令的時(shí)候，表達(dá)式和方程式對(duì)空格都是敏感的，因此不要隨意添加空格符到式中 （2）符號(hào)表達(dá)式的創(chuàng)建：用sym命令直接創(chuàng)建符號(hào)表達(dá)式，這種創(chuàng)建方式不需要在前面有任何說(shuō)明，使用快捷方便。例如： f=sym(a+b+c)在命令窗口中輸入上述語(yǔ)句按回車鍵，出現(xiàn)以下結(jié)果:f= a+b+c。說(shuō)明已成功將符號(hào)表達(dá)式a+b+c賦給變量f5、matlab語(yǔ)言的繪圖功能(1)plot 最基本的二維圖形指令分別完成plot(x)， plot(x,y) ，plot(x1,y1,x2,y2) ，三種格式的調(diào)用，得出結(jié)論。調(diào)用plot(x,y,s)，注意開(kāi)關(guān)量字符串s設(shè)定的曲線顏色和繪圖方式。（2）曲線繪圖觀察 以下

5、各段語(yǔ)句的繪圖結(jié)果：a: x=0, 0.48,0.84,1,0.91,0.6,0.14 x1, x2, x3, x4, x5, x6, x7,plot (x)b: t=0:pi/100:2*pi;y=sin(t);y1=sin(t+0.25);y2=sin(t+0.5);plot(t,y,t,y1,t,y2)c: subplot(1,3,1); plot(t,y)subplot(1,3,2); plot(t,y3)subplot(1,3,3); plot(t,y2)d: t=0:pi/100:2*pi;y=sin(t);y1=sin(t+0.25);y2=sin(t+0.5);plot(t,y

6、,t,y1,t,y2)e: subplot(3,1,1);plot(t,y)subplot(3,1,2);plot(t,y3)subplot(3,1,3);plot(t,y2)f: ezplot(sin(x)ezplot(sin(x),cos(y),-4*pi 4*pi,figure(2)思考題：1、 說(shuō)明產(chǎn)生一個(gè)matlab行向量的兩種方法，其中行向量從0開(kāi)始，并在pi結(jié)束，共有13個(gè)等間隔點(diǎn)。 答：a=0:pi/12:pi b=linspace(0,pi,13)2、假定行向量x=3 6 9,y=5 3 0,下列matlab運(yùn)算的結(jié)果是什么？a x+y 結(jié)果：ans =8 9 9 b, x.

7、*y 結(jié)果： 15 18 0c, x.y 結(jié)果： 243 216 1 d, x./y 結(jié)果： 0.6000 2.0000 inf實(shí)驗(yàn)二 matlab描述常用信號(hào)一、實(shí)驗(yàn)名稱： matlab描述常用信號(hào)二、實(shí)驗(yàn)?zāi)康模?. 在了解matlab這個(gè)軟件的基本應(yīng)用之后，學(xué)習(xí)用matlab描述常用信號(hào)的方法2. 掌握連續(xù)時(shí)間信號(hào)和離散時(shí)間信號(hào)的描述3. 學(xué)會(huì)用向量和符號(hào)表示法4. 掌握信號(hào)的時(shí)域變換方法三、實(shí)驗(yàn)原理：matlab強(qiáng)大的圖形處理功能及符號(hào)運(yùn)算功能，為我們實(shí)現(xiàn)信號(hào)的可視化提供了強(qiáng)有力的工具。在matlab中通常有兩種方法來(lái)表示信號(hào)，一種是用向量來(lái)表示信號(hào)，另一種則是用符號(hào)運(yùn)算的方法來(lái)表示信號(hào)

8、。用適當(dāng)?shù)膍atlab語(yǔ)句表示出信號(hào)后，我們就可以利用matlab的繪圖命令繪制出直觀的信號(hào)波形。1、連續(xù)時(shí)間信號(hào)所謂連續(xù)時(shí)間信號(hào)，是指自變量的取值范圍是連續(xù)的，且對(duì)于一切自變量的取值，除了有若干不連續(xù)點(diǎn)以外，信號(hào)都有確定的值與之對(duì)應(yīng)的信號(hào)。從嚴(yán)格意義上來(lái)講，matlab并不能處理連續(xù)信號(hào)，在matlab中，是用連續(xù)信號(hào)在等時(shí)間間隔點(diǎn)的樣值來(lái)近似地表示連續(xù)信號(hào)的，當(dāng)取樣時(shí)間間隔足夠小時(shí)，這些離散的樣值就能較好地近似出連續(xù)信號(hào)。在matlab中連續(xù)信號(hào)可用向量或符號(hào)運(yùn)算功能來(lái)表示。（1）向量表示法對(duì)于連續(xù)時(shí)間信號(hào)f(t)，我們可以用兩個(gè)行向量f和t來(lái)表示，其中向量t是行如tt1：p：t2的mat

9、lab命令定義的時(shí)間范圍向量，t1為信號(hào)起始時(shí)間，t2為中止時(shí)間，p為時(shí)間間隔。向量f為連續(xù)信號(hào)f(t)在向量t所定義的時(shí)間點(diǎn)上的樣值。例如對(duì)于連續(xù)信號(hào)f(t)=sin(t),我們可以用如下兩個(gè)向量來(lái)表示：t=-10:1.5:10； f=sin(t)用上述向量對(duì)連續(xù)信號(hào)表示后，就可以用plot命令來(lái)繪出該信號(hào)的時(shí)域波形。plot命令可將點(diǎn)與點(diǎn)間用直線連接，當(dāng)點(diǎn)與點(diǎn)間的距離很小時(shí)，繪出的圖形就成了光滑的曲線。命令如下：plot(t,f)title(f(t)=sint)xlabel(t)axis(-10,10,-1.1,1.1)繪制的信號(hào)波形如圖3.1所示，當(dāng)把時(shí)間間隔p取得更小(如0.01)時(shí)，

10、就可得到sint較好的近似波形，如圖3.2所示。圖3.1 p=1.5的sint近似波形 圖3.2 p=0.01的sint近似波形（2）符號(hào)運(yùn)算表示法如果信號(hào)可以用一個(gè)符號(hào)表達(dá)式來(lái)表示它，則我們可用ezplot命令繪制出信號(hào)的波形。例如對(duì)于連續(xù)信號(hào)f(t)= ，我們可以用符號(hào)表達(dá)式表示為：syms tfsym(exp(-t/2)f= exp(-t/2)然后用ezplot命令繪制其波形：ezplot(f,-6, 6)該命令繪制的信號(hào)波形如圖3.3所示圖3.3 指數(shù)信號(hào)波形 利用上面兩種表達(dá)方式我們可以描述出很多種連續(xù)信號(hào)，常用的信號(hào)有：階躍信號(hào)，門(mén)信號(hào)，斜線信號(hào)，指數(shù)信號(hào)等等。二、離散時(shí)間信號(hào)一般

11、說(shuō)來(lái)，離散時(shí)間信號(hào)用f(k)表示，其中變量k為整數(shù)，代表離散的采樣時(shí)間點(diǎn)。f(k)可表示為：f(k)=.f(-2),f(-1),f(0),f(1),f(2) k0在matlab中，用一個(gè)向量f即可表示一個(gè)有限長(zhǎng)度的序列。但是，這樣的向量并沒(méi)有包含其對(duì)應(yīng)的時(shí)間序號(hào)信息。所以，要完整地表示一個(gè)離散信號(hào)需要用兩個(gè)向量。如序列： f(k)=1,2,-1,3,2,4,-1 k=0在matlab中應(yīng)表示為：k=-3,-2,-1,0,1,2,3或是k-3：3；f= 1, 2,-1, 3, 2, 4,-1在用matlab表示離散序列并將其可視化時(shí)，我們要注意以下幾點(diǎn)：第一，與連續(xù)時(shí)間信號(hào)不同，離散時(shí)間信號(hào)無(wú)法

12、用符號(hào)運(yùn)算來(lái)表示；第二，由于在matlab中，矩陣的元素個(gè)數(shù)是有限的，因此，matlab無(wú)法表示無(wú)限序列；第三，在繪制離散信號(hào)波形時(shí)，要使用專門(mén)繪制離散數(shù)據(jù)的stem命令，而不是plot命令。如對(duì)于上面定義的二向量f和k，可用如下stem命令繪圖：stem (k, f)，得到對(duì)應(yīng)序列波形圖，如圖3.4所示。圖3.4 隨機(jī)序列的波形同樣的，單位階躍序列，正弦序列，離散時(shí)間指數(shù)序列等的離散信號(hào)我們都可以用類似的方法描述并繪制出圖形。本節(jié)使用了plot、stem和ezplot這幾個(gè)繪圖命令(顯示連續(xù)信號(hào)中的不連續(xù)點(diǎn)用stairs命令繪圖)，在繪制連續(xù)信號(hào)時(shí)要得到光滑曲線用plot命令；繪制離散信號(hào)則

13、要用stem命令；繪制用matlab符號(hào)表達(dá)式表示的信號(hào)用ezplot命令。 三、信號(hào)的時(shí)域變換信號(hào)的時(shí)域變換包括信號(hào)的平移、反折、倒相及信號(hào)的尺度變換。我們就分別介紹連續(xù)時(shí)間信號(hào)和離散時(shí)間信號(hào)的各種時(shí)域變換11。1、連續(xù)信號(hào)的時(shí)域變換如前所述，matlab可以有兩種方法來(lái)表示連續(xù)信號(hào)。用這兩種方法均可實(shí)現(xiàn)連續(xù)信號(hào)的時(shí)域變換，但用符號(hào)運(yùn)算的方法則較為簡(jiǎn)便。（1）移位對(duì)于連續(xù)信號(hào)f(t),若有常數(shù)t00，延時(shí)信號(hào)f(t-t0)是將原信號(hào)沿正t軸方向平移時(shí)間t0，而f(t+t0)是將原信號(hào)沿負(fù)t軸方向移動(dòng)時(shí)間t0。我們可用下面的命令來(lái)實(shí)現(xiàn)連續(xù)信號(hào)的平移及其結(jié)果可視化，其中f是用符號(hào)表達(dá)式表示的連續(xù)

14、時(shí)間信號(hào)，t是符號(hào)變量，subs命令則將連續(xù)信號(hào)中的時(shí)間變量t用tt0替換：y=subs(f,t,t-t0); ezplot(y)（2）反折連續(xù)信號(hào)的反折，是指將信號(hào)以縱坐標(biāo)為軸反折，即將信號(hào)f(t)中的自變量t換為t。實(shí)現(xiàn)如下：y=subs(f,t,-t); ezplot(y)（3）尺度變換連續(xù)信號(hào)的尺度變換，是指將信號(hào)的橫坐標(biāo)進(jìn)行展寬或壓縮變換，即將信號(hào)f(t)中的自變量t換為at，當(dāng)a1時(shí)，信號(hào)f(at)以原點(diǎn)為基準(zhǔn)，沿橫軸壓縮到原來(lái)的1/a;當(dāng)0a0，延時(shí)信號(hào)f(t-t0)是將原信號(hào)沿正t軸方向平移時(shí)間t0，而f(t+t0)是將原信號(hào)沿負(fù)t軸方向移動(dòng)時(shí)間t0。我們可用下面的命令來(lái)實(shí)現(xiàn)連

15、續(xù)信號(hào)的平移及其結(jié)果可視化，其中f是用符號(hào)表達(dá)式表示的連續(xù)時(shí)間信號(hào)，t是符號(hào)變量，subs命令則將連續(xù)信號(hào)中的時(shí)間變量t用tt0替換： syms t f=sin(t)y=subs(f,t,t-t0); ezplot(y)反折連續(xù)信號(hào)的反折，是指將信號(hào)以縱坐標(biāo)為軸反折，即將信號(hào)f(t)中的自變量t換為t。實(shí)現(xiàn)如下：y=subs(f,t,-t); ezplot(y)3尺度變換和倒相連續(xù)信號(hào)的尺度變換，是指將信號(hào)的橫坐標(biāo)進(jìn)行展寬或壓縮變換，即將信號(hào)f(t)中的自變量t換為at，當(dāng)a1時(shí)，信號(hào)f(at)以原點(diǎn)為基準(zhǔn)，沿橫軸壓縮到原來(lái)的1/a;當(dāng)0a0）的頻譜為： 3.分別介紹以下周期信號(hào)和非周期信號(hào)的

16、頻譜。周期信號(hào)的頻譜周期信號(hào)可以分解成一系列正弦信號(hào)或虛指數(shù)信號(hào)之和。為了直觀地表示信號(hào)所含各分量的振幅，以角頻率為橫坐標(biāo)，以各次諧波的振幅或?yàn)榭v坐標(biāo)，畫(huà)一條直線到x軸，這就形成信號(hào)的幅度譜。類似的可以畫(huà)出角頻率和初相角的線圖，即是信號(hào)的相位譜。周期信號(hào)是功率信號(hào)，研究周期信號(hào)在1電阻上的平均功率，畫(huà)出角頻率與功率的線圖，即是信號(hào)的功率譜。而我們通常所說(shuō)的信號(hào)的頻譜就是指信號(hào)的幅度譜。各譜線頂點(diǎn)連線稱為頻譜包絡(luò)線，它反應(yīng)了各次諧波分量的幅度隨頻率變化的情況。以為縱坐標(biāo)繪制的頻譜是信號(hào)的單邊頻譜，以為縱坐標(biāo)繪制的頻譜是信號(hào)的雙邊頻譜。周期信號(hào)的頻譜只出現(xiàn)在原周期信號(hào)頻率整數(shù)倍的頻率上，它的譜線是

17、離散的，具有以下三個(gè)特點(diǎn)：a離散性，譜線是離散的，只出現(xiàn)在頻率為0，的位置上。b諧波性，譜線只出現(xiàn)在基波頻率的整數(shù)倍上。c收斂性，諧波幅度隨諧波次數(shù)增高而減小。非周期信號(hào)的頻譜的逆變換為： 因?yàn)?所以 上式中第二項(xiàng)的被積函數(shù)是的奇函數(shù)，所以積分為0，因此可以進(jìn)一步化簡(jiǎn)為下式： 上式表明，非周期信號(hào)可以看成是由不同頻率的余弦“分量”疊加而成的，它包含了頻率從零到無(wú)限大的一切分量。由式可見(jiàn)，相當(dāng)于各“分量”的振幅，它是無(wú)窮小量。所以非周期信號(hào)頻譜不能再用幅度表示，而改用密度函數(shù)表示。類似于物質(zhì)的密度是單位體積的質(zhì)量，函數(shù)看作是單位頻率的振幅，故又稱為信號(hào)的頻譜密度函數(shù)，即非周期信號(hào)的頻譜函數(shù)。我們

18、只要求出周期信號(hào)傅立葉級(jí)數(shù)，就可以畫(huà)出信號(hào)的頻譜圖。那么，我們?nèi)绾卫胢atlab畫(huà)出周期信號(hào)的頻譜圖呢，我們?cè)谥耙颜f(shuō)過(guò)，matlab有強(qiáng)大的數(shù)值計(jì)算和繪圖功能，很顯然，繪制頻譜圖之前必須先求出周期信號(hào)的傅立葉系數(shù)，而這一點(diǎn)通過(guò)matlab的積分函數(shù)是很容易實(shí)現(xiàn)的。傅立葉系數(shù)求出以后通過(guò)簡(jiǎn)單的計(jì)算就可以得到或，然后再利用matlab的繪圖函數(shù)就可以畫(huà)出周期信號(hào)的頻譜圖，由于周期信號(hào)的頻譜圖是離散的，所以我們用stem( )函數(shù)即可。下面給出用matlab軟件繪制周期信號(hào)頻譜圖的整個(gè)實(shí)現(xiàn)流程：a編寫(xiě)子函數(shù)y，用符號(hào)表達(dá)式表示出周期信號(hào)在第一個(gè)周期內(nèi)的符號(hào)表達(dá)式，并賦值給返回符號(hào)變量y。b編寫(xiě)子

19、函數(shù)x，求出該周期信號(hào)在繪圖區(qū)間內(nèi)的信號(hào)樣值，并賦值給返回變量x。編寫(xiě)求解信號(hào)傅立葉復(fù)指數(shù)系數(shù)及繪制頻譜圖的通用函數(shù)，該函數(shù)的流程如下： (1) 用函數(shù)y，獲取周期信號(hào)的符號(hào)表達(dá)式。 (2) 求出信號(hào)三角形的傅立葉級(jí)數(shù)展開(kāi)系數(shù)和。 (3) 求出信號(hào)指數(shù)形的傅立葉級(jí)數(shù)展開(kāi)系數(shù)。 (4) 繪制的振幅頻譜圖。 (5) 調(diào)用函數(shù)x，繪制信號(hào)波形圖。主要的matlab算法如下：采用符號(hào)積分int求0,t內(nèi)時(shí)間信號(hào)三角級(jí)數(shù)展開(kāi)系數(shù)，求出，。用forend循環(huán)語(yǔ)句(語(yǔ)法和功能與c語(yǔ)言中的該語(yǔ)句一致)求出三角級(jí)數(shù)展開(kāi)系數(shù)，的具體數(shù)值(n的最大值為諧波最大次數(shù)的值)。求出復(fù)指數(shù)展開(kāi)系數(shù)，因?yàn)?，中n的取值范圍是(

20、0，n)，而中n的取值范圍是(n，n)。所以我們需要用到第一章中提到的反折函數(shù)求出。4.上面的部分用matlab繪制信號(hào)頻譜的算法與實(shí)現(xiàn)流程已經(jīng)全部清楚，余下我們需要學(xué)習(xí)的就是如何編寫(xiě)程序，畫(huà)出信號(hào)的頻譜圖，下面以周期方波脈沖頻譜的matlab實(shí)現(xiàn)作為例子來(lái)講解。給定周期方波脈沖信號(hào)，其幅度為1，脈沖寬度占空比為1/2，周期t5，用matlab編程繪出該信號(hào)頻譜。根據(jù)上節(jié)的實(shí)現(xiàn)流程及算法提示，編寫(xiě)程序如下：function an,bn=fbsbpp% t 輸入信號(hào)周期% 函數(shù)輸入輸出都是數(shù)值量% nf 諧波最大次數(shù)% nn 輸出數(shù)據(jù)的準(zhǔn)確位數(shù)% an cos項(xiàng)系數(shù)，第1項(xiàng)是直流項(xiàng)，以后依次為1

21、，2，3次諧波cos項(xiàng)系數(shù)% bn sin項(xiàng)系數(shù)，第2，3，4元素依次是1，2，3次諧波sin項(xiàng)系數(shù)syms t n k y 定義符號(hào)變量t=5;if nargin2;nf=input(pleas input 所需展開(kāi)的最高諧波次數(shù));endif nargin3;nn=32;end nargin表示輸入?yún)?shù)個(gè)數(shù)y=fb_symfun; 調(diào)用方波的符號(hào)表達(dá)式子函數(shù)，具體函數(shù)見(jiàn)后面的子函數(shù) a0=2*int(y,t,0,t)/t; 根據(jù)定義求a0a=int(2*y*cos(2*pi*n*t/t)/t,t,0,t)；求cos項(xiàng)系數(shù)的符號(hào)積分表達(dá)式b=int(2*y*sin(2*pi*n*t/t)/t

22、,t,0,t); 求sin項(xiàng)系數(shù)的符號(hào)積分表達(dá)式an(1)=double(vpa(a0,nn); 將a0的值以雙精度32位形式傳給an(1)for k=1:nf 循環(huán)語(yǔ)句，每執(zhí)行一次，k加1直到等于nf an(k+1)=double(vpa(subs(a,n,k),nn); 根據(jù)定義求an(n=2nf+1)bn(k+1)=double(vpa(subs(b,n,k),nn);end 根據(jù)定義求bn(n=2,3,nf+1)if nargout=0nargout表示輸出參數(shù)，輸出參數(shù)為零表示所有參數(shù)都已計(jì)算結(jié)束s1=fliplr(an) 對(duì)an陣左右對(duì)稱變換，ana2a(nf+1)，所以s1a(n

23、f+1), a2s1(1,k+1)=an(1) s1陣列擴(kuò)展為k1項(xiàng)，第k1項(xiàng)定為a1，s1a(nf+1), a(nf), a1 s2=fliplr(1/2*s1) s1取半然后反折后給s2，s2a1/2,a2/2,a(nf+1)/2s3=fliplr(1/2*bn) bn取半然后反折后給s3，s3b(nf+1)/2,bnf/2,b2/2 s3(1,k+1)=0 s3陣列擴(kuò)展為k1項(xiàng)，第k1項(xiàng)定為0，s3b(nf+1)/2,bnf/2,b2/2,0 s4=fliplr(s3) 將s3反折后賦給s4， s40,b2/2,b3/2,b(nf)/2,b(nf+1)/2 s5=s2-i*s4; 由三角

24、系數(shù)得出指數(shù)形式的系數(shù)，s5= a1/2,a(nf+1)/2-i b(nf+1)/2s6=fliplr(s5); 對(duì)得到的系數(shù)矩陣反折后傳給s6，s6a(nf+1)/2-i b(nf+1)/2,a1/2n=nf*2*pi/t; n為最大諧波角頻率，等于所需展開(kāi)諧波次數(shù)乘以，k2=-n:2*pi/t:n; 形成n到n的變量序列，變量之間的間隔為，n,n,n,ns7=s6,s5(2:end); 形成n:n的的傅立葉指數(shù)對(duì)稱復(fù)系數(shù)，s7a(nf+1)/2-ib(nf+1)/2 a2/2- ib2/2, a1/2, a2/2-ib2/2a(nf+1)/2-i b(nf+1)/2subplot(2,1,1) 創(chuàng)建子圖1，位于第一行第一列x=fb_pulse 調(diào)用方波脈沖函數(shù)畫(huà)出方