matlab模拟信号受噪声干扰分析,如何用matlab编写噪声调频干扰信号

也就是说，噪声调频信号具有如下的瞬时频率：

f(t)=f0+KFM*u(t)

% 摘自《信息时代的电子战》(29所)

% 1。产生高斯噪声电压

% 2。用6极点椭圆滤波器形成噪声带宽

% 3。将滤波器输出送到一压控震荡器，输出即为调频噪声干扰信号

clf;

clc;

close all;

clear all;

N = 4096;

tmax = 1/2;

t = 2*tmax*(-N/2:N/2-1)/N;

fs= N/2/tmax;

n = rand(1,N);

subplot(2,2,1);plot(t,n);grid;

xlabel('Time');ylabel('Noise Voltage');

title(['Gaussian Noise Voltage']);

y = erf(n);

[b,a] = ellip(6,。

5,。5,。005);

v = filter(b,a,n);

verf = filter(b,a,y);

f0 = 0;k = 50; j=sqrt(-1);

f= f0 + k*v;

f1 = f0 + k*verf;

x = cos(2*pi*f。

*t) + j*sin(2*pi*f。*t);

xerf = cos(2*pi*f1。*t) + j*sin(2*pi*f1。*t);

V = fft(v);

Pv = V。*conj(V);

Pv = fftshift(Pv);

X = fft(x);XE = fft(xerf);

Px = X。

*conj(X);

Pxe = XE。

*conj(XE);

Px = fftshift(Px);

Pxe = fftshift(Pxe);

Px = Px/max(Px);Pv = Pv/max(Pv);Pxe = Pxe/max(Pxe);

Pxdbm = 30 + 10*log10(Px+1e-6);

Pvdbm = 30 + 10*log10(Pv+1e-6);

Pxedbm = 30 + 10*log10(Pxe+1e-6);

fsam = fs*(-N/2:N/2-1)/N;

subplot(2,2,2);

plot(fsam,Pvdbm);grid;

subplot(2,2,3);

plot(fsam,Pxdbm);grid;

subplot(2,2,4);

plot(fsam,Pxedbm);grid;

。

全部