我有以下代码,它基本上绘制了图上的原始信号x
,并迭代地更新重建信号rec
。
plot(x); hold on
err = 100; tol = 0.1; err_vec = [];
while err > tol % iterations
% Low-pass filter xpg
REC = fft(rec);
REC(M+2:N-M) = 0;
rec = real(ifft(REC)); plot(rec, 'r*'); drawnow
% Restore the known samples %
rec(ks) = y(ks);
% Error
err = norm(rec - x)
err_vec = [err_vec err];
end
我喜欢在图上保留x
,并且只在每次重复时更新rec
,以便我可以看到rec
正逐渐接近x
。
但是,使用我当前的代码,虽然保留了x
,但每次迭代的rec
只是重叠在图上,这很烦人。我想仅从当前的迭代中显示rec
。
我应该如何更改代码呢?
答案 0 :(得分:1)