马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
function xrms = Signal_p(x)
t=0:0.001:1;
%x=square(2*pi*1000*t);
n=randn(size(t));%白噪声
f=n;
figure(1);
subplot(2,1,1);
plot(f,'r');%给白噪声作图
ylabel('幅值(v)');
xlabel('时间(s)');
title('原始信号');
y=fft(n,1000);%对信号进行傅里叶变换,参加DFT采样点的个数为1000
subplot(2,1,2);
m=abs(y);
f1=(0:length(y)/2-1)*10000/length(y);%这句什么意思呢?
plot(f1,m(1:length(y)/2));%为什么长度又和上一句不一样了
ylabel('幅值的模');
xlabel('时间(s)');
title('原始信号的傅里叶变换');
%用周期图法估计功率谱密度
p=y.*conj(y)/1000;%计算功率谱密度
ff=10000*(0:499)/1000;%计算变换后不同点对应的频率值,这句的计算公式是怎样的呢
figure(2);
plot(ff,p(1:500));
ylabel('幅值');
xlabel('频率hz');
title('周期谱密度');
将信号由时域转化到频域的时候,是怎样控制频率的长度与信号匹配问题的呢,还有频率的取值是根据什么取得呢?求相关资料推荐!
|