| 查看: 868 | 回复: 2 | ||
lourong88铜虫 (初入文坛)
|
[求助]
球大神帮忙写个matlab程序 已有1人参与
|
|
网上下的程序，出来的图和文献完全不一样，球大神指导 T=1/50; % 基带信号宽度，也就是频率 fc=2/T; % 载波频率 ml=2; % 调制信号类型的一个标志位 nb=100; % 传输的比特数 delta_T=T/200; % 采样间隔 fs=1/delta_T; % 采样频率 SNR=0; % 信噪比 t=0:delta_T:nb*T-delta_T; % 限定t的取值范围 N=length(t); % 采样数 %------------------------- 调制部分 % 基带信号的产生 data=randn(1,nb)>0.5; datanrz=data.*2-1; data1=zeros(1,nb/delta_T); for q=1:nb data1((q-1)/delta_T+1:q/delta_T)=datanrz(q); end % 串并转换，将奇偶位数据分开 idata=datanrz(1:ml  nb-1)); qdata=datanrz(2:ml:nb); % QPSK信号的调制 ich=zeros(1,nb/delta_T/2); for i=1:nb/2 ich((i-1)/delta_T+1:i/delta_T)=idata(i); end for ii=1:N/2 a(ii)=sqrt(2/T)*cos(2*pi*fc*t(ii)); end idata1=ich.*a; qch=zeros(1,nb/2/delta_T); for j1=1:nb/2 qch((j1-1)/delta_T+1:j1/delta_T)=qdata(j1); end for jj=1:N/2 b(jj)=sqrt(2/T)*sin(2*pi*fc*t(jj)); end qdata1=qch.*b; s=idata1+qdata1; ss=abs(fft(s)); h=spectrum.welch; figure(1) subplot(3,2,1); plot(data1),title('基带信号'); axis([0 20000 -2 2]); subplot(3,2,2); psd(h,data1,'fs',fs),title('基带信号功率谱密度'); subplot(3,2,3); plot(s),title('调制信号'); axis([0 500 -3 3]); subplot(3,2,4); psd(h,ss,'fs',fs),title('调制信号功率谱密度');  IMG_1152.JPG  IMG_1153.JPG  IMG_1155.JPG  IMG_1161.JPG |
回复此楼
jxlpyjx
木虫 (正式写手)
小兵向前
- 应助: 0 (幼儿园)
- 金币: 5861.5
- 散金: 66
- 帖子: 926
- 在线: 53.5小时
- 虫号: 3186536
- 注册: 2014-05-07
- 性别: GG
- 专业: 人工智能与知识工程
2楼2016-07-06 00:10:27
sunyuanxin
禁虫 (著名写手)
|
本帖内容被屏蔽 |
3楼2016-07-27 17:30:54
10