基于MATLAB的噪声仿真程序.docx_第1页
基于MATLAB的噪声仿真程序.docx_第2页
基于MATLAB的噪声仿真程序.docx_第3页
基于MATLAB的噪声仿真程序.docx_第4页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

clear all;close all;clc;%清除变量 T=100e-6;%采样时间 fs=300e6;%采样频率N=T*fs;%采样点数 detlf=20e6;%滤波器截止频率f1=100e6;%调制信号中心频率 m=0.5;%调制度 kfm=5e6;%调频斜率 kpm=5;%调相斜率 M=100;%积累次数p=fft(fir1(N-1,detlf/fs*2);%滤波器频谱s=0; for i=1:100 xn=ifft(fft(random(Normal,0,1,1,N).*p);%高斯白噪声通过滤波器 j=abs(fft(xn); end s=s/M; j=s; figure(1) t=0:1/fs:T-1/fs; plot(t*1e6,xn); xlabel(us); title(视频噪声时域波形); figure(2) plot(f*1e-6,20*log10(j.2/max(j.2);%视频噪声功率谱 axis(-1 22 -8 0); xlabel(MHZ); title(视频噪声功率谱)n=1:N; zn=(1+m*cos(2*pi*xn).*cos(2*pi*f1/fs*n);%噪声调幅干扰表达式 figure(3) plot(t*1e6,zn); title(噪声调幅干扰时域波形); xlabel(us); s=0; for i=1:100 zn=(1+m*cos(2*pi*xn).*cos(2*pi*f1/fs*n); j=abs(fft(zn); s=s+j; end s=s/M; j=s; figure(4) plot(f*1e-6,20*log10(j.2/max(j.2);%噪声调幅干扰功率谱title(噪声调幅干扰功率谱); xlabel(MHZ); axis(90 110 -200 0);sum(1)=0; for i=1:N-1sum(i+1)=xn(i)+sum(i); end wn=cos(2*pi*f1*t+2*pi*kfm*xn);%噪声调频干扰表达式figure(5) plot(t*1e6,wn); title(噪声调频干扰时域波形); xlabel(us); s=0; for i=1:100 xn=ifft(fft(random(Normal,0,1,1,N).*p); sum(1)=0; for i=1:N-1; sum(i+1)=xn(i)+sum(i); wn=cos(2*pi*f1*t+2*pi*kfm*xn); j=abs(fft(wn); s=s+j; s=s/M; j=s; figure(6) plot(f*1e-6,20*log10(j.2/max(j.2);%
人人文库> 全部分类> 应用文书 > 研究报告

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论