Matlab 필터 시뮬레이션 오류

Question

나는 위너 필터 시뮬레이션을 matlab에 만들었지 만, 결과가 정확하지 않아서 코드에서 실수를 한 것처럼 보입니다. 누군가 코드를 훑어보고 오류를 지적하면 고맙겠습니다.

k = 1:100; 
sigma = 0.1; 
s1 = randn(1,100); 
n1 = randn(1,100) * sigma; 
Pnum = [1 0.1]; 
Pdenum = [1 0.9]; 
Gnum = [1 0.9]; 
Gdenum = [1 0.1]; 
n = filter(Pnum, Pdenum, n1); 
s = filter(Gnum, Gdenum, s1); 
x = n + s; 
rxx = xcorr(x); 
rxs = xcorr(x,s); 
toeplitz_rxx = toeplitz(rxx); 
wopt = inv(toeplitz_rxx) * rxs; 

wopt = inv(toeplitz_rxx) * transpose(rxs); 
s_hat = filter(wopt,1,x); 
error = x - s_hat; %this is where the wrong result is being calculated 
plot(k,error, 'r',k,s_hat, 'b', k, x, 'g'); 

문제는 산출 오차가 정확하게 s_hat의 미러 이미지이다 : 여기

코드이다. 분명히이 말이 맞을 수 없습니다.

미리 감사드립니다.

Answer 1

문제 없습니다.

나는 당신이 무엇을하려고하는지에 대해 아이디어가 완벽하게 입니다. 그러나 당신이하려고하는 것이 신호 대 잡음비를 얻고, 일부 선형 1 차 시스템을 통해 그 효과를 모방하려한다면, P과 G 다항식을 수정해야합니다. 그렇다면 완전히 틀린 것입니다. help filter을 입력하고 A 및 B 정의를 참조하십시오. 토플 리츠에 병이 조절은 여전히있다

k = 1:100; 
sigma = 0.1; 
s1 = randn(1,100); 
n1 = randn(1,100) * sigma; 
Pnum = [0 0.1]; 
Pdenum = [1 -0.9]; 
Gnum = [0 0.9]; 
Gdenum = [1 -0.1]; 
n = filter(Pnum, Pdenum, n1); 
s = filter(Gnum, Gdenum, s1); 
x = n + s; 
rxx = xcorr(x); 
rxs = xcorr(x,s); 
toeplitz_rxx = toeplitz(rxx); 

wopt = inv(toeplitz_rxx) * transpose(rxs); 
s_hat = filter(wopt,1,x); 
error = x - s_hat; %this is where the wrong result is being calculated 
plot(k,error, 'r',k,s_hat, 'b', k, x, 'g'); 

:

이 올바른 코드입니다. 나는 그런 식으로 두지 만, 당신은 여전히 ​​다루기에 재미가 있습니다. 이것은 쉽게 해결할 수 있습니다.

그들은 사실입니다. StackOverflow는 숙제를 해결하는 것이 아니라 빛을 찾는 데 도움을주기위한 것입니다.