| 查看: 853 | 回复: 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 |
回复此楼» 猜你喜欢
寻找调剂
已经有8人回复
-
289求调剂
已经有16人回复
-
312求调剂
已经有8人回复
-
263求调剂
已经有10人回复
-
一志愿C9材料与化工专业总分300求调剂
已经有6人回复
-
296求调剂
已经有5人回复
-
0703化学/290求调剂/本科经历丰富/工科也可
已经有5人回复
-
284求调剂
已经有4人回复
-
化学工程085602 305分求调剂
已经有12人回复
-
考研调剂
已经有4人回复
sunyuanxin
禁虫 (著名写手)
|
本帖内容被屏蔽 |
3楼2016-07-27 17:30:54
jxlpyjx
木虫 (正式写手)
小兵向前
- 应助: 0 (幼儿园)
- 金币: 5861.5
- 散金: 66
- 帖子: 926
- 在线: 53.5小时
- 虫号: 3186536
- 注册: 2014-05-07
- 性别: GG
- 专业: 人工智能与知识工程
2楼2016-07-06 00:10:27
