版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡介
1、仿真一:證明所有樣本函數(shù)的變換相等。變換的平方的均值與自相關(guān)函數(shù)的%驗(yàn)證樣本函數(shù) clc;clear; x=zeros(1000,1000); a1=zeros(1,1999); y=;%求 1000 個(gè)樣本函數(shù)for k=1:1000變換平方和均值的頻譜等于自相關(guān)函數(shù)的頻譜變換x(k,:)=randn(1, 1000);y(k,:)=fft(x(k,:);end%求 1000 個(gè)樣本函數(shù)for k=1:1000if k=1y1=abs(y(k,:).2;else變換的平方和y1=y1+abs(y(k,:).2;endend y=y1;y=y/1000;%即 1000 個(gè)樣本函數(shù)k=fix(l
2、ength(y)/2); fs=1000;f=fs*(0:k-1)/length(y);變換平方和的均值%畫出 1000 個(gè)樣本函數(shù)figure(1) subplot(2,1,1)plot(f,abs(y(1:k)/length(y)變換平方和均值的頻譜title(樣本函數(shù)變換平方和均值的頻譜)xlabel(Hz)%求 1000 個(gè)樣本函數(shù)各自的自相關(guān)函數(shù)FFT for k=1:1000a,b=xcorr(x(k,:),unbiased); a1=a1+a;yr=fft(a,length(y);%對(duì)自相關(guān)函數(shù)做if k=1c=yr;elsec=c+yr;endend yr=c/1000;變換%
3、求 1000 個(gè)樣本自相關(guān)函數(shù)FFT 的平均值,作為平穩(wěn)序列的自相關(guān)函數(shù)的FFTa1=a1/1000;%求 1000 個(gè)樣本自相關(guān)函數(shù)的平均值,作為平穩(wěn)序列的自相關(guān)函數(shù)%畫出平穩(wěn)序列自相關(guān)函數(shù)的頻譜 kr=fix(length(yr)/2);f=fs*(0:kr-1)/length(yr); subplot(2,1,2) plot(f,abs(yr(1:kr)title(自相關(guān)函數(shù)的頻譜)xlabel(Hz)%畫出自相關(guān)函數(shù)圖形 figure(2)plot(b,a1)title(自相關(guān)函數(shù))%驗(yàn)證樣本函數(shù)變換平方和均值的頻譜等于自相關(guān)函數(shù)的頻譜figure(3) k=1000;kr=lengt
4、h(y);error=abs(y(1:k)/length(y)-abs(yr(1:kr); plot(error)title(error)%輸出最大誤差和平均誤差 display(最大絕對(duì)誤差max_error:) max_error=max(abs(error)display(平均誤差mean_error:) mean_error=mean(error)display(誤差的方差var_error:)var_error=var(error)輸出結(jié)果:最大絕對(duì)誤差max_error: max_error =0.1289平均誤差mean_error:mean_error =-0.0015誤差的方
5、差var_error:var_error=0.0020仿真二:證明所有樣本函數(shù)的 z 變換的平方的均值與自相關(guān)函數(shù)的z 變換相等。%根據(jù)z 變換的定義計(jì)算 z 變換,即對(duì)x(n)z(-n)求和 clc;clear;%初始化z1=zeros(1,1000);%存隨機(jī)序列z 變換 z2=zeros(1,1000);%存自相關(guān)函數(shù)z 變換 a1=zeros(1,1999);y=zeros(1,1000);%隨機(jī)序列z 變換平方的yy=zeros(1,1000);%求 1000 個(gè)樣本函數(shù)Z 變換for k=1:10 x=randn(1,1000);%生成隨機(jī)序列 n=1;a,b=xcorr(x,un
6、biased);%求自相關(guān)函數(shù) a1=a1+a;for z=linspace(1.1,2.1,1000) for i=1:1000z1(n)=z1(n)+x(i)*z(-i);%對(duì)隨機(jī)序列z 變換z2(n)=z2(n)+a(i)*z(-i);%對(duì)自相關(guān)函數(shù)z 變換end n=n+1;end y=y+z1.2; yy=yy+z2;endy=y/10;%求 10 個(gè)樣本 z 變換平方的平均值figure(1) subplot(2,1,1) plot(y)title(樣本z 變換平方的均值)yy=yy/10;%求 10 個(gè)樣本自相關(guān)函數(shù) z 變換的平均值,作為平穩(wěn)序列自相關(guān)函數(shù)的z 變換subplo
7、t(2,1,2) plot(abs(yy)title(樣本自相關(guān)函數(shù)的z 變換)figure(2) plot(b,a1/10)title(樣本的自相關(guān)函數(shù))輸出結(jié)果:輸出結(jié)果誤差比較大,以下是選取誤差較小的圖仿真三:證明cross-spectra 等于各自變換的乘積%驗(yàn)證cross-spectra 等于各自clc;clear; x=zeros(10,1000); y=zeros(10,1000); xx=zeros(1,1000); yy=zeros(1,1000); for k=1:10 x(k,:)=randn(1,1000);y(k,:)=randn(1,1000); x+x(k,:);
8、yy=yy+y(k,:);end%x=randn(1,1000);%y=randn(1,1000);/10;y=yy/10; X=fft(x);Y=fft(y); XY=X.*Y;a,b=xcorr(x,y,unbiased);變換的乘積r=fft(a,length(XY);%對(duì)自相關(guān)函數(shù)做變換figure(1) k=fix(length(XY)/2); fs=1000;f=fs*(0:k-1)/length(XY); subplot(2,1,1) plot(f,abs(XY(1:k)/length(XY)title(各自變換的乘積)kr=fix(length(r)/2); f=fs*(0:k
9、r-1)/length(r); subplot(2,1,2) plot(f,abs(r(1:kr) title(cross-spectra)%驗(yàn)證cross-spectra 等于各自 figure(2) k=length(XY);kr=length(r);error=abs(XY(1:k)/length(XY)-abs(r(1:kr); plot(error)title(error)變換的乘積%輸出最大誤差和平均誤差 display(最大絕對(duì)誤差max_error:) max_error=max(abs(error)display(平均誤差mean_error:) mean_error=mea
10、n(error)display(誤差的方差var_error:) var_error=var(error)輸出結(jié)果:最大絕對(duì)誤差max_error: max_error =0.5811平均誤差mean_error: mean_error =-0.1137誤差的方差var_error: var_error =0.0143仿真四:白噪聲、平滑濾波%白噪聲、平滑濾波(中值濾波) clc;clear;N=10000;x=randn(1,N);%產(chǎn)生均值為 0,方差為 1 的figure(1)隨機(jī)序列,也就是白噪聲序列plot(x);%白噪聲title(白噪聲)n=N/10;a,b=xcorr(x,n,
11、biased); %求自相關(guān)函數(shù) figure(2)plot(b,a)%畫出自相關(guān)函數(shù) title(自相關(guān)函數(shù))yr1=fft(a);%求自相關(guān)函數(shù) y=abs(yr1);%求功率譜 figure(3)plot(y);title(功率譜密度) axis(0 2*n -5 5)變換%中值濾波 yr2=y;d=6;k=length(yr2)-d; for i=1:kxr=yr2(i:i+d); yr2(i)=median(xr);end figure(4)plot(yr2(1:length(yr2)-d)%畫出中值濾波后的功率譜密度axis(0 2*n -5 5)title(中值濾波后的功率譜密度)display(中值濾波前的均值和方
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 2023-2024學(xué)年湖南省長沙市高一上學(xué)期入學(xué)考試英語試題(解析版)
- 成長作文之成長的故事作文教案
- 2024年電能計(jì)量配套產(chǎn)品合作協(xié)議書
- 2024年汽車?yán)鋮s風(fēng)扇項(xiàng)目建議書
- 2024年船用舾裝件項(xiàng)目合作計(jì)劃書
- 2024年防靜電海綿項(xiàng)目建議書
- 2024年配電控制設(shè)備項(xiàng)目發(fā)展計(jì)劃
- 2024年放射性污染防治項(xiàng)目合作計(jì)劃書
- 2024年生產(chǎn)L型氨基酸的新酶種項(xiàng)目合作計(jì)劃書
- 高端醫(yī)療區(qū)渣土運(yùn)輸服務(wù)
- 2024年房產(chǎn)贈(zèng)與合同范本(31篇)
- 2024年中國移動(dòng)校園招聘高頻難、易錯(cuò)點(diǎn)500題模擬試題附帶答案詳解
- 《雕梁畫棟》課件-2024-2025學(xué)年人美版(2024)初中美術(shù)七年級(jí)上冊(cè)
- 物理透鏡 課件 2024-2025學(xué)年蘇科版八年級(jí)上冊(cè)物理
- 2024秋四年級(jí)英語上冊(cè) Unit 5 Dinners ready第2課時(shí)(Let's learn Lets play)教案 人教PEP
- 【智慧醫(yī)療】醫(yī)療健康產(chǎn)業(yè)園概念策劃方案(XQ)
- 2024年計(jì)算機(jī)二級(jí)WPS考試題庫380題(含答案)
- 2024中國移動(dòng)公司招聘高頻500題難、易錯(cuò)點(diǎn)模擬試題附帶答案詳解
- 2.1 網(wǎng)絡(luò)改變世界 課件-2024-2025學(xué)年統(tǒng)編版道德與法治八年級(jí)上冊(cè)-1
- 青島市教育局智慧教育中心平臺(tái)建設(shè)方案
- 22G101三維彩色立體圖集
評(píng)論
0/150
提交評(píng)論