1. matlab帮助文档中通过计算方差的形式来计算OFDM一个符号的功率

1):>>doc 《QPSK and OFDM with MATLAB System Objects》

代码如下:为什么是利用方差呢?直接利用如下公式不行吗?对比结果发现,相差很少。

txSig_erg=sum(abs(txSig).^2);10*log10(txSig_erg/length(txSig));
dataIn = randi([0,1],frameSize);              % Generate binary data
qpskTx = qpskMod(dataIn);                     % Apply QPSK modulation
txSig = ofdmMod(qpskTx);                      % Apply OFDM modulation
powerDB = 10*log10(var(txSig));               % Calculate Tx signal power

注意: 功率和功率普密度不是同一个概念哦!!

2. 如何计算一个OFDM符号的功率谱密度

有两种计算方式,请参考如下matlab帮助文档

1):>>doc  periodogram函数,有详细的参考公式

2):>>doc  《Power Spectral Density Estimates Using FFT》

%% 对比两种计算功率谱密度的方法,一个采用FFT,一个采用periodogram函数
%% 以下是计算实信号的功率谱密度
rng default
Fs = 1000;
t = 0:1/Fs:1-1/Fs;
x = cos(2*pi*100*t) + randn(size(t));N = length(x);
xdft = fft(x);
xdft = xdft(1:N/2+1);
psdx = (1/(Fs*N)) * abs(xdft).^2;
psdx(2:end-1) = 2*psdx(2:end-1); % 此处为什么要乘以2呢???
freq = 0:Fs/length(x):Fs/2;plot(freq,10*log10(psdx))
grid on
title('Periodogram Using FFT')
xlabel('Frequency (Hz)')
ylabel('Power/Frequency (dB/Hz)')periodogram(x,rectwin(length(x)),length(x),Fs)
mxerr = max(psdx'-periodogram(x,rectwin(length(x)),length(x),Fs))%% 以下计算复数信号的功率谱密度
rng default
n = 0:999;
x = exp(1j*pi/4*n) + [1 1j]*randn(2,length(n))/sqrt(2);%% Use fft to obtain the periodogram. Because the input is complex-valued, obtain the
%% periodogram from  rad/sample. Plot the result.N = length(x);
xdft = fft(x);
psdx = (1/(2*pi*N)) * abs(xdft).^2;% 计算复信号时没有乘以2!!!这是为什么???
freq = 0:(2*pi)/N:2*pi-(2*pi)/N;
%%
plot(freq/pi,10*log10(psdx))
grid on
title('Periodogram Using FFT')
xlabel('Normalized Frequency (\times\pi rad/sample)')
ylabel('Power/Frequency (dB/rad/sample)')
%%
periodogram(x,rectwin(length(x)),length(x),'twosided')
mxerr = max(psdx'-periodogram(x,rectwin(length(x)),length(x),'twosided'))%% 归一化频率时的计算方式如下
Input with Normalized Frequency
Use fft to produce a periodogram for an input using normalized frequency. Create a signal consisting of a sine wave in N(0,1) additive noise. The sine wave has an angular frequency of  rad/sample. Use the default settings of the random number generator for reproducible results.
%%
rng default
n = 0:999;
x = cos(pi/4*n) + randn(size(n));
%%
Obtain the periodogram using fft. The signal is real-valued and has even length. Because the signal is real-valued, you only need power estimates for the positive or negative frequencies. In order to conserve the total power, multiply all frequencies that occur in both sets -- the positive and negative frequencies -- by a factor of 2. Zero frequency (DC) and the Nyquist frequency do not occur twice. Plot the result.
%% N = length(x);
xdft = fft(x);
xdft = xdft(1:N/2+1);
psdx = (1/(2*pi*N)) * abs(xdft).^2;
psdx(2:end-1) = 2*psdx(2:end-1);
freq = 0:(2*pi)/N:pi;
%% plot(freq/pi,10*log10(psdx))
grid on
title('Periodogram Using FFT')
xlabel('Normalized Frequency (\times\pi rad/sample)')
ylabel('Power/Frequency (dB/rad/sample)')%% Compute and plot the periodogram using periodogram. Show that the two results are identical.periodogram(x,rectwin(length(x)),length(x))
mxerr = max(psdx'-periodogram(x,rectwin(length(x)),length(x)))
mxerr = 1.4211e-14

二者计算的结果是一样的,第二种方法更加具体,是利用FFT来计算,相当于是第一种方法的具体实现过程

OFDM专题之如何计算OFDM一个符号的功率,功率谱密度相关推荐

  1. OFDM专题之输入的复信号从何而来?

    本博文承接上篇博文:OFDM专题之原理图理解 由于是专题系列,所以每篇博文都是有相互关系的,不是独立的,所以看这些博文要一篇一篇的看, 萌新也要多看几遍,不是一天多看几遍,而是很多天的很多次,一次有的 ...

  2. 【OFDM系列2】OFDM复信号符号信噪比EsN0、比特信噪比EbN0、SNR的含义及关系详解

    更基础的各类信噪比的关系请看博客 [干货2020.12.13更新]实信号.复信号Es.N0.符号信噪比EsN0.带内信噪比.比特信噪比EbN0.SNR的含义及关系详解 本文只讨论OFDM复信号的信噪比 ...

  3. matlab实现ofdm调制解调,基于ofdm的调制解调原理,ofdm调制解调fpga代码 - 全文

    OFDM介绍 OFDM ( Orthogona lFrequency D iv isionM ulTIplex ing) 即正交频分复用技术, 实际上 OFDM 是 MCM M ult-iC arri ...

  4. OFDM和F-OFDM的功率谱与峰均功率比仿真

    OFDM和F-OFDM的功率谱与峰均功率比仿真 王家尉1 (1.湖南大学通信工程系,湖南长沙 410082;) 摘  要: OFDM(Orthogonal Frequency Division Mul ...

  5. 云栖小镇不是“镇”,就像中关村不是“村”。小镇是一个符号,就像起建于50年前的硅谷的“谷”,和100年前爱迪生所在的门洛公园。...

    云栖小镇不是"镇",就像中关村不是"村".小镇是一个符号,就像起建于50年前的硅谷的"谷",和100年前爱迪生所在的门洛公园. 转载于:ht ...

  6. 【小题目】输入两个数字以及一个符号,输出这两个数字在这个符号下运算产生的结果

    import java.util.Scanner; public class SwitchCaseExer {public static void main(String[] args){// 获取两 ...

  7. php去掉最后一个,号,织梦当前位置去除最后一个''符号

    dede提供了面包屑的功能.但是最后面总会带一个'>'符号.甚是烦人.目标干掉它. 第一种解决方法:   1.  在系统参数里面修改这个符号. 第二种解决方法:   2. 使用截子符,runph ...

  8. 字数统计,一个英文字母算一个字,一个汉字算一个字,一个符号算一个字

    //============================================================================= // // 字数统计,一个英文字母算一个 ...

  9. 【日记本砸】21.06.11-20 复杂的式子和角标只是一个符号一个标记而已

    喜欢 就是我想把所有的安全感和喜欢都有所克制.恰如其分地传达给你, 不让你失落也不给你压力. 互联网+吹牛大赛 破事不值得浪费时间. 梦 生活费 我以为自己能结余,结果入不敷出- 买伞 是个小天才,想 ...

  10. C++基础编程题(27)输入一个数字,为其高,一个符号*,输出该符号组成的平行四边形形状

    C++基础编程题(27)输入一个数字,为其高,一个符号*,输出该符号组成的平行四边形形状

最新文章

  1. linux自学笔记——RAID级别特性以及软RAID的实现
  2. LNMP 环境遇到的权限问题
  3. java arraylist排序_最全Java集合笔记
  4. 2020-python小工能
  5. python算法应用(六)——搜索与排名2(PageRank算法及其拓展应用)
  6. 第三节课-损失函数和优化
  7. matlab超限像素平滑法_教程 | 不做手残党,玩转像素风LOGO的三大实用技巧!
  8. 学术资源不定期分享-【钱学森《工程控制论》英文原版】
  9. pyqt5写一个抽奖程序(再次美化)
  10. DSS与Schedulis的集成安装教程(单机)
  11. 二项分布期望与方差的证明
  12. jquery-seat-charts 使用-自定义座位号 及 重新加载数据
  13. mqtt 3.1 php代码,MQTT V3.1
  14. Android Studio创建app问题: Install repository and sync project等
  15. 敏之澳电商:入驻拼多多开店流程及费用多少?
  16. C语言 查找书籍(结构体)
  17. 如何破解百度文库的防复制功能
  18. Ubuntu系统下创建本地源
  19. gitlab复制project新建项目
  20. 中国城市码cityCode映射

热门文章

  1. 云原生|kubernetes|etcd集群详细介绍+安装部署+调优
  2. 好看的皮囊 · 也是大自然的杰作 · 全球高质量 · 美图 · 集中营 · 美女 · 2017-08-24期...
  3. IntelliJ IDEA 在使用manven后的纠结(每次修改代码都要重启tomcat才能看效果吗?),请各位大侠来看看问题
  4. 打印系统开发(9)——办公打印机使用教程
  5. Intel VT学习笔记(八)—— 开启EPT
  6. android+system+disk.vdi,VirtualBox压缩VDI格式硬盘/调整磁盘大小
  7. Visio 导出图片时字符间距错乱
  8. 有关企业合并时的10条IT措施
  9. 建筑企业收并购系列二:吸收合并政策影响
  10. 超微晶磁芯在开关电源中的应用