【求助】关于复数开方的问题，兄弟姐妹帮各忙吧，我得毕设卡在这下不去了：（

passmatlab

你使用fft不对呀！还要使用还是函数fftshift。
只有两者的配合使用才能够得到你的结果。

molidong

我在www.matlab-world.com里早就回答你了。。。。

molidong

clear all
N=300  % no use
for t=0:1:299
ha(t+1)=exp(-0.05*t)*(1-exp(-0.05*t));   
end
subplot(2,1,1);
plot(0:1:299,ha);      
for t=0:1:299
Ha2(t+1)=fft(ha(t+1))*fft(ha(t+1));
end
for t=0:1:299
haa(t+1)=ifft(sqrt(Ha2(t+1))) ;
end
subplot(2,1,2);
plot(0:1:299,haa) ;                     
%e_mail:      molidong@yahoo.com.cn

molidong

没说的，叫我一声大哥我就知足了，不过我不明白，实际上你的程序（包括你原来的程序）始终是对常数求傅里叶变换及反变换，不知你这样做的意义何在。。

molidong

Ha2=fft(ha).*fft(ha) ％频域相乘％： 你再这里出现了复数，而不是实数，导致后面错了。原因我也不知道。。

molidong

噪声不过是随机数，你加上就是了，而且对于fft 和ifft 来说，不论是实数，还是复数，都等于自身，应该没什末问题。。fft(3+2j)=3+2j  
ifft(3+2j)=3+2j

molidong

噪声只是一个随机数，用rand(), randn()就能实现了，复数噪声我想你可以这样实现 K*rand()+j*w*rand()