【正文】 。應(yīng)用條件應(yīng)用重力勘探的條件是﹕被探測(cè)的地質(zhì)體與圍巖的密度存在一定的差別﹔被探測(cè)的地質(zhì)體有足夠大的體積和有利的埋藏條件﹔干擾水平低。由于重力反演存在多解性﹐因此﹐必須依靠研究地區(qū)的地質(zhì)﹑鉆井﹑巖石密度和其他物探資料來(lái)減少反演的多解性。定量解釋是根據(jù)異常場(chǎng)求激發(fā)體的產(chǎn)狀要素建立重力模型。對(duì)稱異常曲線說(shuō)明質(zhì)量相對(duì)于通過(guò)極值點(diǎn)的垂直平面是對(duì)稱分布的﹔反之﹐非對(duì)稱曲線是由于質(zhì)量非對(duì)稱分布引起的。最大的水平梯度異常相應(yīng)于激發(fā)體的邊界。除某些構(gòu)造外﹐對(duì)一般地質(zhì)體重力異常的解釋可遵循以下的一些原則﹕極大的正異常說(shuō)明與圍巖比較存在剩馀質(zhì)量﹔反之﹐極小異常是由質(zhì)量虧損引起的。第五篇：信號(hào)與系統(tǒng)實(shí)驗(yàn)報(bào)告1長(zhǎng)江大學(xué)電工電子實(shí)驗(yàn)中心實(shí)驗(yàn)報(bào)告野外獲得的重力數(shù)據(jù)要作進(jìn)一步處理和解釋才能解決所提出的地質(zhì)任務(wù)﹐主要分3個(gè)階段﹕野外觀測(cè)數(shù)據(jù)的處理﹐并繪制各種重力異常圖﹔重力異常的分解(應(yīng)用平均法﹑場(chǎng)的變換﹑頻率濾波等方法)﹐即從疊加的異常中分出那些用來(lái)解決具體地質(zhì)問(wèn)題的異常﹔確定異常體的性質(zhì)﹑形狀﹑產(chǎn)狀及其他特征參數(shù)。七、實(shí)驗(yàn)結(jié)果：2101***222101234567八、實(shí)驗(yàn)心得體會(huì)：通過(guò)本次試驗(yàn)，了解了頻率響應(yīng)求法，加深了對(duì)輸出穩(wěn)態(tài)響應(yīng)的印象。subplot(2,1,1),plot(t,u1)。u2=u21+u22+23。u22=abs(U(2)*H(2))*cos(100*t+angle(U(2)*H(2)))。subplot(2,1,2),plot(w,angle(H)),grid。H=freqs(b,a,w)。u1=U*cos(w’*t+angle(U’)*ones(1,1001))。b=[2,0]。w=[99,100,101]。六、實(shí)驗(yàn)程序：clear。實(shí)驗(yàn)五一、實(shí)驗(yàn)者姓名：二、實(shí)驗(yàn)時(shí)間：三、實(shí)驗(yàn)地點(diǎn)：四、實(shí)驗(yàn)題目：2s，激勵(lì)信號(hào)22(s+1)+100x(t)=(1+cot)s*co1s0(t)0，求（1）帶通濾波器的頻率響應(yīng)；（2）輸出穩(wěn)態(tài)響應(yīng)并繪制圖形。通過(guò)計(jì)算機(jī)做出的信號(hào)波形圖，我們能夠很直白的看出原波形和經(jīng)過(guò)變換后的波形的差別。plot(t,y)。title(‘原波形’)。subplot(2,1,1)。x=sin(2*pi*120*t)。對(duì)于典型的離散時(shí)間信號(hào)，可用邏輯表達(dá)式來(lái)實(shí)現(xiàn)不同自變量時(shí)的取值。當(dāng)序列是從0開始時(shí)，可以只用一個(gè)向量x來(lái)表示序列。其中一個(gè)向量表示自變量的取值范圍，另一個(gè)向量表示序列的值。本實(shí)驗(yàn)中要求繪制離散時(shí)間信號(hào)圖，可以應(yīng)用MATLAB中的函數(shù)來(lái)實(shí)現(xiàn)。五、解題分析：此題需要對(duì)信號(hào)做離散序列，還要做FFT和IFFT變換，然后得到圖像進(jìn)行比較。九、實(shí)驗(yàn)改進(jìn)想法：無(wú)。七、實(shí)驗(yàn)結(jié)果：***0100806040200|F(k)|05101520Frequency253035***0100806040200|F(k)|05101520Frequency253035 4321|jW|0123405101520Frequency253035Step Response4000200002000 Amplitude40006000800010000120001400001234n(samples)5678八、實(shí)驗(yàn)心得體會(huì)：通過(guò)本次試驗(yàn)我知道了求取模擬信號(hào)在N等于一定值時(shí)的的DFT的幅值譜和相位譜的求法。plot(q,angle(y))。title(‘magnitude’)。subplot(3,1,2)。plot(t,x)。y=fft(x,N)。q=n*2*pi/N。n=0:63。close。若MN，則將原序列截短為N點(diǎn)序列，再計(jì)算其N點(diǎn)IDFT；若M六、實(shí)驗(yàn)程序：clc。M是序列X的長(zhǎng)度。fft(x,N)計(jì)算N點(diǎn)的DFT。fft(x)計(jì)算M點(diǎn)的DFT。實(shí)驗(yàn)三一、實(shí)驗(yàn)者姓名：二、實(shí)驗(yàn)時(shí)間：三、實(shí)驗(yàn)地點(diǎn)：四、實(shí)驗(yàn)題目：模擬信號(hào)x(t)=2sin(4pt)+5cos(8pt)，求N=64的DFT的幅值譜和相位譜。并知道了在《信號(hào)與系統(tǒng)》中得一些差分方程和各種響應(yīng)，譬如零輸入相應(yīng)、零狀態(tài)響應(yīng)、全響應(yīng)、自由響應(yīng)、強(qiáng)迫響應(yīng)、沖擊響應(yīng)、單位階躍響應(yīng)等等各種響應(yīng)在Matlab中的函數(shù)表達(dá)方式和他們的求法，以及系統(tǒng)的幅頻和相頻曲線的繪制都有了一定深刻的認(rèn)識(shí)。p rad/sample)(degrees) Frequency(180。figure(2)。dimpulse(b,a,20)subplot(2,1,2)。figure(1)。b=[1 01]。clear。注意一般情況下，函數(shù)freqz的返回值X是復(fù)數(shù)。返回值X就是DTFTb=[b0,b1,...,bM]。另外，可用MATLAB提供的abs、angle、real、imag等基本函數(shù)計(jì)算 DTFT的幅度、相位、實(shí)部、虛部。輸出序列的長(zhǎng)度和序列相同。其中 a=[a0,a1,...,aN]，b=[b0,b1,...,bM]，分別表示差分方程系數(shù)。(bn/b0)x[kn]n=1n=0NM利用MATLAB提供的filter函數(shù)，可方便地計(jì)算出上述差分方程的零狀態(tài)響應(yīng)。在已知差分方程的N個(gè)初始狀態(tài)y[k]，和輸入x[k]，就可由下式迭代計(jì)算出系統(tǒng)的輸出y[k]=229。ay[kn]=229。由離散序列卷積的性質(zhì)可知，當(dāng)序列x和h的起始點(diǎn)都為k=0時(shí)，y的取值范圍為k=0至length(x)+length(h)2。其中調(diào)用參數(shù)x,h為卷積運(yùn)算所需的兩個(gè)序列，返回值y是卷積結(jié)果。實(shí)驗(yàn)二一、實(shí)驗(yàn)者姓名：二、實(shí)驗(yàn)時(shí)間：三、實(shí)驗(yàn)地點(diǎn)：四、實(shí)驗(yàn)題目:一個(gè)因果線性移不變系統(tǒng)y(n)=(n2)+x(n)x(n2)，求：（1）H(z)；（2）沖激響應(yīng)h(n)；（3）單位階躍響應(yīng)u(n)；（4）H(ejw)，并繪出幅頻和相頻特性。為后面的實(shí)驗(yàn)打下基礎(chǔ)，并對(duì)信號(hào)仿真和《信號(hào)與系統(tǒng)》這門課程之間的聯(lián)系有所增加，對(duì)《信號(hào)與系統(tǒng)》這門課里的問(wèn)題也有了更加深入地了解。title(‘Step response’)七、實(shí)驗(yàn)結(jié)果：實(shí)驗(yàn)所得波形圖如下：Step (sec)678910八、實(shí)驗(yàn)心得體會(huì)：通過(guò)本次試驗(yàn)了解學(xué)會(huì)了一些新的函數(shù)的應(yīng)用。den=[1 6 10 8]。s+6s2+10s+8五、解題分析：要知道求單位階躍響應(yīng)需知道所用函數(shù)，以及產(chǎn)生波形圖所需要用到的函數(shù)。y=ｙ＊x。ｘ1=sym（＇Hｅavｉｓidｅ(t+2)’)。％產(chǎn)生一個(gè)寬度ｔao=1 得矩形脈，中心位置在ｔ＝５處 pｌot（t,ｘ)hoｌd ｏn x=rｅctpｕls(t＋5,ｔaｏ）。t=—２＊T:0、01：2*T； tao＝T/5。它由符號(hào)變量與表達(dá)式寫成。%subploｔ（3，3，3）%x＝tｉmｅ_fuｎ_ｅ(ｔ)。k2=—N:２＊pi／T:N。%B＿sｙｍ得 1＊k 陣擴(kuò)展為１*(k+１)陣Bn=fliｐlr（Bｎ)。%對(duì)擴(kuò)展后得 S１陣左右對(duì)稱交換回原位置Bn＝fｌiｐlr（B_ｓym)。end if nａrgout==0Ａn=flｉｐlr(Ａ_sym)；%對(duì) A_sym 陣左右對(duì)稱交換An（１，k+1）=A_sｙm(1)。for k=1：ＮfＡ_sym(k+１)＝double(vpa（subs(Ａs,n，k),Nn））。Bｓ＝2/T＊iｎt(ｙ*ｓin(2＊pi*n＊t/Ｔ）,t,０,Ｔ／2）。ｅnd y=tiｍe_ｆun_s(t)； A0＝2/Ｔ＊int（ｙ,t，0,T／2)。endT=ｉnpuｔ(’ｐleaｓe Iｎput 信號(hào)得周期 T＝”）。x＝x1－sym(’Heaｖisidｅ（t1）’）＊ｈ；已知周期為Ｔ=4 得三角波,在第一周期(2functｉon [A＿sym,B_sym]=CTFSｓhbpsyｍ(T，Nｆ)% 采用符號(hào)計(jì)算求［０，T]內(nèi)時(shí)間函數(shù)得三角級(jí)數(shù)展開系數(shù)。=1。ｅ1＝1/2+1/２、＊sｉgｎ(ｔ—０、５＋tａo／２)； e2＝１／2+1/２、*sign(t—0、5—tａo／2)； y=h、＊(e1—e2）。ploｔ(t，y,”r:’)tｉtlｅ(’1—１0 次諧波+直流“)axｉs（[－8，8,0、5，5］)hｏｌd off ｅnd fｕnｃtiｏn y=time_fun_e(t)% 該函數(shù)就是 CＴFShｃhsym、m 得子函它由符號(hào)函數(shù)與表達(dá)式寫成 a＝5； T=４； h＝１； taｏ=Ｔ/4。％調(diào)用連續(xù)時(shí)間函數(shù)－周期矩形脈沖pｌoｔ(t，y，”r：“）titｌe（”直流+基波’)axiｓ([８，８，0、５，5])sｕbploｔ（２,2,２）plot(t,f1２),hold oｎy=time_fun_e（t）。% 9 次諧波f １ 0=c（11)、*co ｓ(２ ＊ pi*10*t/T）+d(1 １)、*s ｉ n(2*pi ＊ 1 ０ ＊t/T)；％ 1０次諧波ｆ１1=f0+f１+f2；% 直流+基波＋2 次諧波f１2=f11+f3。% ５次諧波f6=c（7)、*cｏｓ（２*pi*6＊t/Ｔ）+d(７)、*sin（2＊pｉ*6*t/T)；％ 6 次諧波ｆ 7=c(8）、*cos(2*p ｉ ＊ ７ *t/T）+d(8）、＊ sin(2 ＊ ｐ ｉ ＊ 7 ＊t/T)；% 7 次諧波f8=c(9)、＊ｃos（2*ｐi＊8*ｔ/T)+d(9）、*sin（2*ｐi＊8＊t/T）。% 3次諧波f4=ｃ（５)、＊cos（2*ｐi*4＊t/T）+d（５)、＊siｎ(2*pi＊４*t/Ｔ)。%直流ｆ 1 ＝ c（2)、＊ co ｓ（2* ｐ i* １ * ｔ /T)＋ d(2）、＊ s ｉ n(2 ＊ pi* １ * ｔ /Ｔ）；% 基波ｆ 2= ｃ(３)、* ｃ o ｓ（2*pi ＊ 2 ＊ t/T)+d(3）、*sin（2 ＊ pi ＊ 2* ｔ ／T）。%輸出 d 為三角級(jí)數(shù)展開系數(shù)：第 2,3，4,、元素依次就是 1,２,次諧波ｓｉn 項(xiàng)展開系數(shù)ｔ=—3＊T:0、０１:3*T。%輸出 c 為三角級(jí)數(shù)展開系數(shù):第 1 元素就是直流項(xiàng)，其后元素依次就是 1,2,次諧波ｃos 項(xiàng)展開系數(shù) d=B＿ｓym。%獲取串?dāng)?shù)組Ｂ所對(duì)應(yīng)得 AＳC2 碼數(shù)值數(shù)組 end；if nａrgoｕｔ==0c=A＿syｍ。A0=iｎt(x,t，a,Ｔ+a）/Ｔ；％求出三角函數(shù)展開系數(shù)Ａ0 Ａｓ＝2/T*iｎｔ（x＊coｓ（2*pi*ｎ*t/Ｔ),t,a，T+ａ)；％求出三角函數(shù)展開系數(shù) As Bｓ=2/Ｔ*ｉnｔ(x*sin（2*pｉ＊n＊t/T)，t，a，T+ａ)；%求出三角函數(shù)展開系數(shù) Bs Ａ＿sym(1)=ｄoubｌe（vｐａ(A0,Ｎn）)；％獲取串?dāng)?shù)組 A0 所對(duì)應(yīng)得 ASC2碼數(shù)值數(shù)組 fｏr k＝１:Nf A＿sym(k＋1)＝ｄouble（ｖpａ(subｓ(Ａs，n，k）,Nn)）。ｉf ｎarｇｉnNf=10； enｄ if nａrｇｉn＜5Ｎｎ=３2。tao=T/4。=165。 229。242。（１）三角形式傅立葉級(jí)數(shù)dt t n t fTbdt t n t fTadt t fTat n b t n a at b t a t b t a t b t a a t fTT nTT nTTnnnnn n n n242。title（” “ 濾波器輸出波形’))實(shí)驗(yàn) 四周期信號(hào)得頻域分析１..實(shí)驗(yàn)?zāi)康胠 掌握周期信號(hào)傅立葉級(jí)數(shù)分解與合成得計(jì)算公式 l 掌握利用 MATＬAＢ實(shí)現(xiàn)周期信號(hào)傅立葉