zstp_sig.m

这是关于“噪声调频信号”的仿真模型

% 摘自《信息时代的电子战》（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;