OFDM专题之如何计算OFDM一个符号的功率,功率谱密度
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一个符号的功率,功率谱密度相关推荐
- OFDM专题之输入的复信号从何而来?
本博文承接上篇博文:OFDM专题之原理图理解 由于是专题系列,所以每篇博文都是有相互关系的,不是独立的,所以看这些博文要一篇一篇的看, 萌新也要多看几遍,不是一天多看几遍,而是很多天的很多次,一次有的 ...
- 【OFDM系列2】OFDM复信号符号信噪比EsN0、比特信噪比EbN0、SNR的含义及关系详解
更基础的各类信噪比的关系请看博客 [干货2020.12.13更新]实信号.复信号Es.N0.符号信噪比EsN0.带内信噪比.比特信噪比EbN0.SNR的含义及关系详解 本文只讨论OFDM复信号的信噪比 ...
- matlab实现ofdm调制解调,基于ofdm的调制解调原理,ofdm调制解调fpga代码 - 全文
OFDM介绍 OFDM ( Orthogona lFrequency D iv isionM ulTIplex ing) 即正交频分复用技术, 实际上 OFDM 是 MCM M ult-iC arri ...
- OFDM和F-OFDM的功率谱与峰均功率比仿真
OFDM和F-OFDM的功率谱与峰均功率比仿真 王家尉1 (1.湖南大学通信工程系,湖南长沙 410082;) 摘 要: OFDM(Orthogonal Frequency Division Mul ...
- 云栖小镇不是“镇”,就像中关村不是“村”。小镇是一个符号,就像起建于50年前的硅谷的“谷”,和100年前爱迪生所在的门洛公园。...
云栖小镇不是"镇",就像中关村不是"村".小镇是一个符号,就像起建于50年前的硅谷的"谷",和100年前爱迪生所在的门洛公园. 转载于:ht ...
- 【小题目】输入两个数字以及一个符号,输出这两个数字在这个符号下运算产生的结果
import java.util.Scanner; public class SwitchCaseExer {public static void main(String[] args){// 获取两 ...
- php去掉最后一个,号,织梦当前位置去除最后一个''符号
dede提供了面包屑的功能.但是最后面总会带一个'>'符号.甚是烦人.目标干掉它. 第一种解决方法: 1. 在系统参数里面修改这个符号. 第二种解决方法: 2. 使用截子符,runph ...
- 字数统计,一个英文字母算一个字,一个汉字算一个字,一个符号算一个字
//============================================================================= // // 字数统计,一个英文字母算一个 ...
- 【日记本砸】21.06.11-20 复杂的式子和角标只是一个符号一个标记而已
喜欢 就是我想把所有的安全感和喜欢都有所克制.恰如其分地传达给你, 不让你失落也不给你压力. 互联网+吹牛大赛 破事不值得浪费时间. 梦 生活费 我以为自己能结余,结果入不敷出- 买伞 是个小天才,想 ...
- C++基础编程题(27)输入一个数字,为其高,一个符号*,输出该符号组成的平行四边形形状
C++基础编程题(27)输入一个数字,为其高,一个符号*,输出该符号组成的平行四边形形状
最新文章
- linux自学笔记——RAID级别特性以及软RAID的实现
- LNMP 环境遇到的权限问题
- java arraylist排序_最全Java集合笔记
- 2020-python小工能
- python算法应用(六)——搜索与排名2(PageRank算法及其拓展应用)
- 第三节课-损失函数和优化
- matlab超限像素平滑法_教程 | 不做手残党,玩转像素风LOGO的三大实用技巧!
- 学术资源不定期分享-【钱学森《工程控制论》英文原版】
- pyqt5写一个抽奖程序(再次美化)
- DSS与Schedulis的集成安装教程(单机)
- 二项分布期望与方差的证明
- jquery-seat-charts 使用-自定义座位号 及 重新加载数据
- mqtt 3.1 php代码,MQTT V3.1
- Android Studio创建app问题: Install repository and sync project等
- 敏之澳电商:入驻拼多多开店流程及费用多少?
- C语言 查找书籍(结构体)
- 如何破解百度文库的防复制功能
- Ubuntu系统下创建本地源
- gitlab复制project新建项目
- 中国城市码cityCode映射
热门文章
- 云原生|kubernetes|etcd集群详细介绍+安装部署+调优
- 好看的皮囊 · 也是大自然的杰作 · 全球高质量 · 美图 · 集中营 · 美女 · 2017-08-24期...
- IntelliJ IDEA 在使用manven后的纠结(每次修改代码都要重启tomcat才能看效果吗?),请各位大侠来看看问题
- 打印系统开发(9)——办公打印机使用教程
- Intel VT学习笔记(八)—— 开启EPT
- android+system+disk.vdi,VirtualBox压缩VDI格式硬盘/调整磁盘大小
- Visio 导出图片时字符间距错乱
- 有关企业合并时的10条IT措施
- 建筑企业收并购系列二:吸收合并政策影响
- 超微晶磁芯在开关电源中的应用