【正文】 )(2cos)(2)(1)sin()cos(...sin cos...sin cos sin cos)( ww w w w w w（２）指數形式傅立葉級數（b）求解，及合成信號波形所用程序: fｕnction [Ａ_sｙm,B_sym]＝ＣTFShｃhsym ％ 采用符號計算求一個周期內連續(xù)時間函數 f 得三角級數展開系數，再用這些 %展開系數合成連續(xù)時間函數ｆ、傅立葉級數 ％ 量值數是就都出輸入輸得數函?%數階得波諧 6=fN?% Nn數位確準得據數出輸?% 數系開展項 soc 波諧次、3,２,１是就次依素元后其，項流直是就素元 1 第?myｓ_A?% B_sym 第 2，３，４,、元素依次就是 1，2，次諧波 siｎ項展開系數 %taｏ=１tａo/T=0、2 sｙｍs t n k xT＝4。ａ=５。end134 5 4 1 Ox＝timｅ＿ｆun_x（t)。％獲取串數組Ａ所對應得 ASC2碼數值數組 B_sym（k+1)=ｄｏｕble(vpa（subs(Ｂs，ｎ，k）,Nｎ））。disp(ｃ）。disp(ｄ)。ｆ 0= ｃ(１)。% 2 次諧波f3=ｃ(4）、＊ｃos(2＊pi*3＊t／T)+d(4)、＊sin(2*ｐi*3＊t/T)。% ４次諧波f５=c（6）、*ｃos(2*ｐi*5*ｔ/T）+d(6）、*ｓｉn(2＊ｐi*5*t／T)。% 8次諧波f9 ＝ ｃ(1 ０)、＊ c ｏ ｓ(２ * ｐ i*9 ＊ ｔ /T)+d(10)、* ｓ in(2 ＊ ｐ i ＊ 9 ＊t/T）。％ 直流+基波+2 次諧波+3 次諧波f13＝f12＋f4+f5+f6；% 直流＋基波+2 次諧波+３次諧波+４次諧波+5 次諧波+6 次諧波f14＝f13+f７+f8+ｆ9+f１0；%0～10 次subｐlot(2,2，1)ｐlｏt(t，f0＋f1）,hoｌｄ ｏny＝tiｍe_ｆun＿ｅ(t)。ｐlot(t,y，’r：’）tiｔle(“1—3 次諧波+直流”）axis([—8,8,0、5,5］）sｕｂploｔ(２,2,3)ploｔ(t,f13），hoｌd ony＝ｔime_fun＿e（t)；plot(ｔ,y,’ｒ:’）ｔitle（“１—6 次諧波＋直流＇)axis(［8，8，－0、５,5］)ｓubpｌot（２,2，4)plｏt（t,ｆ14）,holｄ ony=timｅ_fun_e（t）。t=—3*Ｔ:0、0１:3*T。%連續(xù)時間函數—周期矩形脈沖 funcｔion x＝tiｍｅ_fun＿ｘ（t）% 該函數就是 CTFSｈcｈsym、ｍ得子函數。ｘ1=sｙm(”Heaｖisidｅ(ｔ)“)*h。%? 函數得輸入輸出都就是數值量 ％ Ｎｎ輸出數據得準確位數 % mys_A? 第１元素就是直流項,其后元素依次就是１,２,次諧波 coｓ項展開系數 ％ Ｂ_ｓyｍ 數系開展項 nis 波諧次、３,２,1 是就次依素元、,４，3,2 第?％TT=m*taｏ, 信號周期 ％ ? ｆN? 諧波得階數 ％m(m=T/tao)周期與脈沖寬度之比,如 m=4，8，16,1００等 ％tao脈寬：ｔａo＝Ｔ/ｍsｙｍs tnｙiｆ naｒgin＜3Nｆ=inｐut（’plｅａse Inｐut 所需展開得最高諧波次數:Nf＝’)。if nａrｇｉn〈5Ｎn=32。Ａs=2/Ｔ＊iｎｔ（y*cｏs（2*pi*n＊t/Ｔ）,ｔ,0,T/2)。Ａ＿sym(1)=double(vpa(Ａ０，Ｎｎ））。B＿sym(k+1)=ｄouble(ｖｐa(sｕbs(Bs，n，k),Ｎn))。%Ａ＿sym 得 1*k 陣擴展為 1＊（k+1)陣An=ｆｌｉplr（Aｎ)。%對 B_sym 陣左右對稱交換Ｂn（1,ｋ+1)=0。%對擴展后得 S３陣左右對稱交換回原位置FnR=Aｎ/2—i＊Ｂn/２；% 用三角函數展開系數 A、Ｂ值合成付里葉指數系數ＦnL=fｌｉplr（FnＲ）；N=Nf*2*pi/T。Fn＝[FnＬ，F(xiàn)nＲ（2:eｎｄ）]。% 調用連續(xù)時間函數周期矩形脈沖sｕｂplot(2,1，1)stｅm(k２,aｂs(Fn)）；％畫出周期矩形脈沖得頻譜(T=M*tao）ｔiｔle(＇連續(xù)時間函數周期三角波脈沖得雙邊幅度譜’)axｉs([－8０，80,0,０、12]）ｌine（[80，8０］,[0,０],＇colｏr“，’ｒ”）ｌｉｎe([０,０],[0，０、１２]，’colｏｒ’,＇r“）end ｆunctioｎ x=time＿fun_e(t)% 該函數就是ＣTFSｓhbpsｙm、ｍ得子函數。％ t組數間時是就?％ T 0=T/oat=yｔuｄ期周是就??T=５。ｘ=rectpuｌs（t,tａo)；％產生一個寬度 tａo=1 得矩形脈沖 sｕbplot(2,2,2）ｐlｏt(ｔ,x）ｈoｌｄ on x=rectｐuls(t—５，tao)。%產生一個寬度ｔao＝１得矩形脈，中心位置在 t=—５處 ploｔ(t，x)ｔitlｅ(”周期為 T＝5,脈寬 tａo=１得矩形脈沖＇)ａxis([10,10,0,2］)funｃtｉoｎ y＝time_fｕｎ＿s(t）syms t y＝1—abs(t)。ｘ=x1—sｙｍ(Heａvisidｅ(t2）’)。ｅzplｏt(t，y,[—10,1０])grid第四篇：信號與系統(tǒng)仿真實驗報告信號與系統(tǒng)仿真實 驗 報 告班級： 學號： 姓名： 學院：實驗一一、實驗者姓名：二、實驗時間：三、實驗地點：四、實驗題目：5(s2+5s+6)求三階系統(tǒng)H(s)=3的單位階躍響應，并繪制響應波形圖。六、試驗程序：num=[5 25 30]。step(num,den,10)。了解到了N階系統(tǒng)的單位階躍響應的計算方法，和系統(tǒng)的響應波形圖的函數應用和繪制方法。九、實驗改進想法：無。五、解題分析：離散卷積是數字信號處理中的一個基本運算，MTLAB提供的計算兩個離散序列卷積的函數是conv，其調用方式為 y=conv(x,h)。MATLAB函數conv的返回值y中只有卷積的結果，沒有y的取值范圍。許多離散LTI都可用如下的線性常系數的差分方程描述229。bx[kn]nnn=0n=0NN其中x[k]、y[k]分別系統(tǒng)的輸入和輸出。(an/a0)y[kn]+229。filter函數調用形式為 y=filter(b,a,x)。X表示輸入序列，y表示輸出序列。當序列的DTFT可寫成ejW的有理多項式時，可用MATLAB信號處理工具箱提供的freqz函數計算DTFT的抽樣值。若X（ejW）可表示為b0+b1ejW+...bMejWMB(ejW)X(e)==jWjWjWNA(e)a0+a1e+...+aNe則freqz的調用形式為 X=freqz(b,a,w)，其中的b和 a分別是表示前一個jW式子中分子多項式和分母多項式系數的向量，即a=[a0,a1,...,aN]，w為抽樣的頻率點，向量w的長度至少為2。在抽樣點w上的值。六、實驗程序:clc。close。a=[1 ]。subplot(2,1,1)。dstep(b,a,50)w=[0:1:512]*pi/512。freqz(b,a,w)七、實驗結果：沖擊響應圖及階躍響應圖:Impulse (sec)Step (sec)3035404550 100Magnitude(dB) Frequency(180。p rad/sample)八、實驗心得體會：通過實驗我們知道了使用Matlab來繪出出一個線性移不變系統(tǒng)的幅頻和相頻曲線。九、實驗改進想法：無。五、解題分析：在MATLAB信號處理工具箱中，MATLAB提供了4個內部函數用于計算DFT和IDFT，它們分別是：fft(x)，fft(x,N)，ifft(X)，ifft(X,N)。M是序列x的長度，即M=length(x)。若MN，則將原序列截短為N點序列，再計算其N點DFT；若Mifft(X)計算M點的IDFT。ifft(X,N)計算N點IDFT。clear。N=64。t=d*n。x=2*sin(4*pi*t)+5*cos(8*pi*t)。subplot(3,1,1)。title(‘source signal’)。plot(q,abs(y))。subplot(3,1,3)。title(‘phase’)。通過本次實驗，對幅值譜和相位譜有了更深的了解，并與課程《信號與系統(tǒng)》里的一些相關知識連接到了一起，使得學到的只是更加深刻、有意義。實驗四一、實驗者姓名：二、實驗時間：三、實驗地點：四、實驗題目：將信號x(t)=sin(240pt)做離散序列，比較原序列與經過FFT和IFFT變換后的序列，并做出說明。連續(xù)時間函數與離散時間函數在編程中的區(qū)別主要體現(xiàn)在如下兩個方面：第一，自變量的取值范圍不同，離散時間函數的自變量是整數，而連續(xù)時間函數的自變量為一定范圍內的實數；第二，繪圖所用的函數不同，連續(xù)函數圖形的繪制不止一個。用MATLAB表示一離散序列，可用兩個向量來表示。之后畫出序列波形。由于計算機內寸的限制，MATLAB無法表示一個無窮長的序列。六、實驗程序：t=0:1/255:1。y=real(ifft(fft(x)))。plot(t,x)。subplot(2,1,2)。七、實驗結果：八、實驗心得體會：通過對做信號的離散序列以及經FFT和IFFT的變換，了解了相關特性。九、實驗改進想法：無。已知帶通濾波器的系統(tǒng)函數為H(s)=五、解題分析：需要知道求頻率響應的方法，并繪制圖形。t=linspace(0,2*pi,1001)。U=[,1,]。a=[1,2,10001]。H=polyval(b,j*w)./polyval(a,j*w)。subplot(2,1,1),plot(w,abs(H)),grid。u21=abs(U(1)*H(1))*cos(99*t+angle(U(1)*H(1)))。u23=abs(U(3)*H(3))*cos(101*t+angle(U(3)*H(3)))。figure(2)。subplot(2,1,2),plot(t,u2)。九、實驗改進想法：無。解釋解釋分為定性的和定量的兩個內容﹐定性解釋是根據重力圖并與地質資料對比﹐初步查明重力異常性質和獲得有關異常源的信息?？拷|量重心﹐在地表投影處將觀測到最大異常。延伸異常相應于延伸的異常體﹐而等軸異常相應于等軸物體在地表的投影。在平面上出現(xiàn)幾個極值的復雜異常輪廓﹐表明存在幾個非常接近的激發(fā)體。一種常用的反演方法是選擇法﹐即選擇重力模型使計算的重力異常與觀測重力異常間的偏差小于要求的誤差。應用運用領域在區(qū)域地質調查﹑礦產普查和勘探的各個階段都可應用重力勘探﹐要根據具體的地質任務設計相應的野外工作方法。意義 重力勘探解決以下任務﹕研究地殼深部構造﹔研究區(qū)域地質構造﹐劃分成礦遠景區(qū)﹔掩蓋區(qū)的地質填圖﹐包括圈定斷裂﹑斷塊構造﹑侵入體等﹔廣泛用于普查與勘探可燃性礦床(石油﹑天然氣﹑煤)﹐查明區(qū)域構造﹐確定基底起伏﹐發(fā)現(xiàn)鹽丘﹑背斜等局部構造﹔普查與勘探金屬礦床(鐵﹑鉻﹑銅﹑多金屬及其他)﹐主要用于查明與成礦有關的構造和巖體﹐進行間接找礦﹔也常用于尋找大的﹑近地表的高密度礦體﹐并計算礦體的儲量﹔工程地質調查﹐如探測巖溶﹐追索斷裂破