返回信息流我想请教一下,我在做循环平稳信号的高阶统计量的仿真时,用了循环体但是没有体现出循环。而且变量的幅值也是没有体现。clc
clear
close all
N = 128;
n = [0:N-1];
Nfft =128;
fs = 12000;
t = n/fs;
A=1;%定义输入信号幅值
fn=3000;
fa=12000;
fr=8000;
fc=10000;
fA=20;
P=1/20;
T=1/120;
y=randn(size(t));
e1=mean(y);%%%%均值方差
e1=0;
Ai=zeros(1,128*3);
e2=mean(Ai);
e2=P;
i=3;
for i=1:3
for j=1:length(t)
Ai(j+128*(i-1))=cos(2*pi*fA*t(j)+30)+1+y(j);
end
end
a=0.01*P;
mu=0;
sigma=a;
tao=normrnd(mu,sigma,1,128);
B=2;%定义衰减震荡函数幅值
w=zeros(1,128*3);
h=zeros(1,128*3);
i=3;
for i=1:3
for j=1:length(t)
w(j+128*(i-1))=t(j)-i*P-tao(j);
h=B*exp(-30*w).*sin(2*pi*fn*w);%定义衰减震荡函数
end
end
figure
plot(w,h);%输出t2-h图像
xlabel('w','fontsize',13);%x轴为t2,字号13
ylabel('h','fontsize',13);%y轴为h,字号13
title('h=exp(-w/3).*sin(2*pi*fn*w)');%显示标题
s=Ai.*h;
figure(2)
plot(w,s);%输出t2-h图像
xlabel('w','fontsize',13);%x轴为t2,字号13
ylabel('h','fontsize',13);%y轴为h,字号13
title('h=exp(-w/3).*sin(2*pi*fn*w)');%显示标题
这是一条镜像帖。来源:北邮人论坛 / matlab / #11457同步于 2014/12/29
Matlab机器人发帖
求助循环平稳信号的仿真??
liumi1125
2014/12/29镜像同步0 回复
订阅后,新回复会通过你的通知中心匿名送达。
0 条回复
暂无回复 · 你可以订阅本帖等待新回复。