首页 Matlab中求功率谱密度的代码(DOC)

Matlab中求功率谱密度的代码(DOC)

举报
开通vip

Matlab中求功率谱密度的代码(DOC)Matlab中求功率谱密度的代码(DOC) [code] clear; %随机产生长度大于1000的‘0’、‘1’信号序列,对其进行QPSK调制 %定义待仿真序列的维数 N global N N=2000; %定义产生‘1’的概率为 p global p p=0.5; %产生随机二进制序列 s_qpsk=randsrc(1,N,[1,0;p,1-p]); %画出生成的随机序列图 figure(1); stem(s_qpsk); axis([0 50 -0.5 1.5]); xlabel(...

Matlab中求功率谱密度的代码(DOC)
Matlab中求功率谱密度的代码(DOC) [code] clear; %随机产生长度大于1000的‘0’、‘1’信号序列,对其进行QPSK调制 %定义待仿真序列的维数 N global N N=2000; %定义产生‘1’的概率为 p global p p=0.5; %产生随机二进制序列 s_qpsk=randsrc(1,N,[1,0;p,1-p]); %画出生成的随机序列图 figure(1); stem(s_qpsk); axis([0 50 -0.5 1.5]); xlabel('维数N') ylabel('信号强度') title('0/1等概分布的信号') %********QPSK信号的数字调制******** [m_qpsk1,m_qpsk2]=qpsk_modulation(s_qpsk); figure(2); plot(m_qpsk1,m_qpsk2,'r*'); axis([-2 2 -2 2]); title('QPSK的信号空间图'); %********插值,相邻信号间插入7个零点******** insert_qpsk1=upsample(m_qpsk1,8); insert_qpsk2=upsample(m_qpsk2,8); %画出插值后的序列 figure(3); subplot(2,1,1); plot(insert_qpsk1(1:90),'ro'); axis([0 100 -1.5 1.5]); hold on; plot(insert_qpsk1(1:90)); xlabel('实部信号'); axis([0 100 -1.5 1.5]); title('QPSK插值后序列'); subplot(2,1,2); plot(insert_qpsk2(1:90),'yo'); axis([0 100 -1.5 1.5]); hold on; plot(insert_qpsk2(1:90)); xlabel('虚部信号'); axis([0 100 -1.5 1.5]); %********升余弦滤波器滤波******** out_qpsk1=rise_cos(insert_qpsk1,N,8*N); out_qpsk2=rise_cos(insert_qpsk2,N,8*N); %画出滤波后的信号 figure(5); subplot(2,1,1); n=1:100; plot(n,out_qpsk1(1:100),'.-r'); hold on; m=25:104; stem(m,insert_qpsk1(1:80),'o'); legend('滤波输出信号','输入信号'); title('通过平方根升余弦滤波器滤波得到QPSK实部输出信号(10个周期)'); subplot(2,1,2); plot(n,out_qpsk2(1:100),'.-r'); hold on; stem(m,insert_qpsk2(1:80),'y'); legend('滤波输出信号','输入信号'); title('通过平方根升余弦滤波器滤波得到QPSK虚部输出信号(10个周期)'); %********输出信号眼图******** %滤波后两路信号合并, 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 达成复数形式 eyediagram(out_qpsk1,5*8); title('QPSK实部眼图'); eyediagram(out_qpsk2,5*8); title('QPSK虚部眼图'); %********输出信号功率谱密度******** out_qpsk=out_qpsk1+i*out_qpsk2; R_I=xcorr(out_qpsk); power_qpsk=fft(R_I); figure(8); plot(10*log10(abs(power_qpsk(1:(length(power_qpsk)+1)/2)))-max(10*log10(abs(power_qpsk(1: (length(power_qpsk)+1)/2))))); grid on; xlabel('频率'); ylabel('dB'); title('QPSK功率谱密度'); [/code] 此程序来源于网络 作者:蜗_牛
本文档为【Matlab中求功率谱密度的代码(DOC)】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑, 图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
该文档来自用户分享,如有侵权行为请发邮件ishare@vip.sina.com联系网站客服,我们会及时删除。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
下载需要: 免费 已有0 人下载
最新资料
资料动态
专题动态
is_769014
暂无简介~
格式:doc
大小:15KB
软件:Word
页数:3
分类:
上传时间:2017-12-03
浏览量:87