求用維納濾波方法處理雜訊的matlab程式必須

時間 2022-03-16 12:41:59

1樓:匿名使用者

假設乙個點目標在x,y平面上繞單位圓做圓周運動,由於外界干擾,其運動軌跡發生了偏移。其中,x方向的干擾為均值為0,方差為0.05的高斯雜訊;y方向干擾為均值為0,方差為0.

06的高斯雜訊。

1、試設計一fir維納濾波器,確定最佳傳遞函式,並用該濾波器處理觀測訊號,得到其最佳估計。(注:自行設定誤差判定閾值,根據閾值確定濾波器的階數或傳遞函式的長度)。

4、 附件的實驗報告中給出了解題思路,實現源程式、以及結果分析,分別繪製了x方向和y方向的期望訊號、雜訊訊號、觀測訊號、濾波後訊號、最小均方誤差訊號的曲線圖;

6、繪製了期望訊號、觀測訊號和濾波後點目標的運動軌跡。

2樓:匿名使用者

clear all

clck=8;

fs = 1000; %取樣率

n = 1000; %取樣點數

n = 0:n-1;

t = 0:1/fs:1-1/fs; %時間序列

signal=sin(2*pi*10*t);

noise=sin(2*pi*100*t); %前500點高斯分部白雜訊,後500點均勻分布白雜訊

xn= signal+k*noise; %構造的混合訊號

figure(1)

plot(xn); grid on;

title('原始訊號 ');

fs=fft(xn,512); %將訊號變換到頻域

afs=abs(fs); %訊號頻域圖的幅值

f=(0:255)*fs/512; %頻率取樣

figure(2)

plot(f,afs(1:256)); %濾波前的訊號頻域圖

grid on;

xlabel('頻率/赫茲');ylabel('幅度');title('訊號濾波前頻域圖');

figure(3);

rxx=xcorr(xn,xn); %得到混合訊號的自相關函式

m=100; %維納濾波器階數

for i=1:m %得到混合訊號的自相關矩陣

for j=1:m

rxx(i,j)=rxx(abs(j-i)+n);

endend

rxy=xcorr(xn,signal); %得到混合訊號和原訊號的互相關函式

for i=1:m

rxy(i)=rxy(i+n-1);

end %得到混合訊號和原訊號的互相關向量

h = inv(rxx)*rxy'; %得到所要涉及的wiener濾波器係數

signal_filter=filter(h,1, xn); %將輸入訊號通過維納濾波器

plot(signal_filter);grid on;

title('維納濾波後的訊號');

fs2=fft(signal_filter,512); %將訊號變換到頻域

afs2=abs(fs2); %訊號頻域圖的幅值

f=(0:255)*fs/512; %頻率取樣

figure(4)

plot(f,afs2(1:256)); %濾波前的訊號頻域圖

grid on;

xlabel('頻率/赫茲');ylabel('幅度');title('訊號濾波後頻域圖');

figure(5)

subplot(221)

plot(xn);grid on;

title('原始訊號 ');

subplot(222)

plot(f,afs(1:256)); %濾波前的訊號頻域圖

grid on;

xlabel('頻率/赫茲');ylabel('幅度');title('訊號濾波前頻域圖');

subplot(223)

plot(signal_filter);grid on;

title('維納濾波後的訊號');

subplot(224)

plot(f,afs2(1:256)); %濾波前的訊號頻域圖

xlabel('頻率/赫茲');ylabel('幅度');title('訊號濾波後頻域圖');

grid on;

c一道題目,求程式,用遞迴方法求解

include stdio.h define nlength 20 int itoa int nin,char sout int main itoa a,buf printf s n buf return 0 最多就只用了乙個printf.應該還是比較容易理解的吧,就不多做解釋了 include i...

用快速冪方法求組合數,怎麼用計算器求排列組合 求具體步驟

事實上,m e mod n可以直接計算,沒有必要先算m e。m e mod n叫做模取冪運算,根據簡單的數論知識,很容易設計乙個分治演算法。具體如下 設是整數b的二進位制表示 即b的二進位制有k 1位,b k 是最 高位 下列過程隨著c的值從0到b成倍增加,最終計算出a c mod n modula...

用什麼辦法就能在月內暴瘦2019斤求方法

那是不可能的吧,想要健康的 乙個月最多估計也就是10斤吧。不能只要求快速,欲速而不達,本身就是乙個漫長的過程 好好堅持就行。就要像 任涵去脂 一樣,堅持自己的運動和飲食 健康的去 重慶博視眼科 一頓就吃一碗飯 蔬菜,早上2根黃瓜,8點一根,10點一根。晚飯吃水果。適當跑跑步,堅持半個月見效,而且不傷...

高一求函式解析式什麼時候用什麼方法

一共有十一種 一 定義法,二 換元法,三 方程組法,四 特殊化法,五 待定係數法,六 函式性質法,七 反函式法,八 即時定義 法,九 建模法,十 影象法,十 一 軌跡法。但在高一階段只有六種 一 配湊法 把形如 f g x 內的 g x 當做整體,在解析式的右端整理成只含有 g x 的形式,再把 g...

補水保濕乳液使用方法?求介紹,補水乳液怎麼用

在給自己的臉上塗抹保濕乳液的時候,大家應該要時刻保證做到哪一塊區域比較缺水就先塗在哪一塊區域的原則。如果大家的臉上沒有塗抹眼霜的話,就先將保濕乳液塗抹在雙眼周圍的區域,接著便是臉部的兩側臉頰 唇部兩邊 唇部的上側與下側。福瑞達美業旗下生產的保濕乳液不同系列新增了不同的植物提取物,但主要成分都是透明質...