fsignal.m

来自「这是一个对信号和噪声进行分析的程序」· M 代码 · 共 33 行

function signal=fsignal()
b=0.15;
a=15;
f=14.65;
h=pi/2;
s=zeros(101,1);
y1=zeros(101,1);
y=zeros(101,1);
n=1;
for t=0:0.01:1
if t>b
s(n)=20*sin(2*pi*f*t+h)*exp(-a*(t-b));
else
    s(n)=0;
end
y1(n)=2*sin(2*pi*t+pi/8);

n=n+1;
end
y2=s+y1;
t1=0:0.01:1;
figure(1)
plot(t1,s,t1,y2)
legend('冲击信号','加正弦信号后的信号');
figure(2)
y= awgn(y2,10,'measured'); % Add white Gaussian noise. 
plot(t1,y2,t1,y) % Plot both signals. 
legend('Original signal','Signal with AWGN');
figure(3)
plot(t1,y)
title('欲处理信号')
signal=y;