完善资料让更多小伙伴认识你,还能领取20积分哦, 立即完善>
| 回复本帖可获得 1 分积分奖励! 每人限 1 次 | |
|
%plot_CCDF.m
%绘制图7.3中的CCDF曲线 clear;clc;clf Ns=2.^[6:10];%这里是使用的每一个元素相乘,也就是向量或者矩阵中对应元素相乘64,128,256,512,1024 b=2; M=2^b; Nblk=1e3; zdBs=[4:0.1:10]; N_zdBs=length(zdBs);% 61 CCDF_formula = inline('1-((1-exp(-z.^2/(2*s2))).^N)','N','s2','z');%内联函数 for n=1:length(Ns) N=Ns(n);% 64,128,256,512,1024 x=zeros(Nblk,N);%产生一个NBLK行 N列的全零矩阵 sqN=sqrt(N);%开平方根 for k=1:Nblk X=mapper(b,N);%size(X) 64 4 x(k,:)=ifft(X,N,1)*sqN;%%冒号是取出第k行的数据, CFx(k)=PAPR(x(k,:)); end s2=mean(mean(abs(x)))^2/(pi/2);%平均值 CCDF_theoretical=CCDF_formula(N,s2,10.^(zdBs/20));%理论值 for i=1:N_zdBs CCDF_simulated(i)=sum(CFx>zdBs(i))/Nblk;%仿真值 end semilogy(zdBs,CCDF_theoretical,'k-')%表示理论值Y坐标轴是对数坐标系 hold on,grid on semilogy(zdBs(1:3:end),CCDF_simulated(1:3:end),'k:*')%表示仿真值Y坐标是对数坐标系 end axis([zdBs([1 end]) 1e-2 1])%axis([xmin xmax ymin ymax]) 用来标注输出的图线的最大值最小值。 title ('OFDM system with N-point FFT') xlabel('PAPAR0[dB]') ylabel('CCDF=probability(PAPR>PAPR0)') legend('Theoretical','simulated')%添加图例的标注,显示理论值和仿真值 在运行后一直报错 颜色标注错误行数 错误提示是:赋值具有的非单一 rhs 维度多于非单一下标数 出错 CCDF (line 17) x(k,:)=ifft(X,N,1)*sqN 求各位大神帮忙看下好吧,谢谢各位了! |
|
相关推荐
3 个讨论
|
|
|
求各位大侠帮忙看看呀
|
|
|
|
|
|
|
|
回帖奖励 +1 分积分
顶一下,很厉害的样子
|
|
|
|
|
|
|
|
你正在撰写讨论
如果你是对讨论或其他讨论精选点评或询问,请使用“评论”功能。
我用matlab的coder 封装了一个c语言的dll, 但是在用labview调用时,会出现识别不到库的问题,有大神遇到过吗
2862 浏览 0 评论
5281 浏览 0 评论
在matlab中如何计算含有第一类修正的贝塞尔函数的积分算不出的问题?
9079 浏览 0 评论
怎么利用matlab得到95%,80%和70%的置信区间,并生成不同区间下的功率误差贝塔分布?
10523 浏览 0 评论
请问simulink的s-function模块如何添加多输入输出接口
14705 浏览 2 评论
/9
小黑屋| 手机版| Archiver| 电子发烧友 ( 湘ICP备2023018690号 )
GMT+8, 2025-12-4 05:50 , Processed in 0.711121 second(s), Total 65, Slave 51 queries .
Powered by 电子发烧友网
© 2015 bbs.elecfans.com
关注我们的微信
下载发烧友APP
电子发烧友观察
版权所有 © 湖南华秋数字科技有限公司
电子发烧友 (电路图) 湘公网安备 43011202000918 号 电信与信息服务业务经营许可证:合字B2-20210191

淘帖
8476