|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
function y=xinzaobi(s)
clc
close all
D=0.6;
a=1;
b=1;
fs=10; %参样频率,可调参数
f0=0.01;%信号频率,可调参数
Ts=1/fs;%参样时间
h=1/fs; %时间步长
t=0:Ts:4095*Ts;
s=0.3*sin(2*pi*f0*t); %信号
x0=s+sqrt(2*D)*randn(size(t)); %噪声,强度为D
%四阶龙格库塔法对双稳态输出信号求解
x1=sr(a,b,h,x0);
%输出信号求傅立叶变换
y=fft(x1,4096);
pyy=y.*conj(y)/4096,
ff=fs*(0:2048)/4096;
S=abs(pyy).^2;
N=sum(abs(pyy).^2)-abs(pyy).^2;
SNR=10*log10(S/N);
plot(ff,SNR)
这个程序输出的只是么一个频率下的SNR曲线,并且是一条直线,我怎么用matlab去画出SNR与噪声强度D之间的关系曲线呢,还望高手教教我啊 |
|