求教怎么仿真sigmadelta调制器功率谱密度

yjj_123

目前仿真代码就是参考网上那种fft求功率谱密度的方法，但是结果总是和论文里不一样

wkp1992101

X轴坐标换成和论文一样的就可以

yjj_123

wkp1992101 发表于 2021-8-18 09:22
X轴坐标换成和论文一样的就可以

纵坐标要不要也改一下，差的也有点多感觉

yjj_123

wkp1992101 发表于 2021-8-18 09:22
X轴坐标换成和论文一样的就可以

兄弟，你有代码可以给我参考一下吗

spartan313

你X 轴整体除以10^7 看一下

yjj_123

geo24 发表于 2021-8-18 11:01
你X 轴整体除以10^7 看一下

只改横坐标不太行啊

spartan313

yjj_123 发表于 2021-8-18 11:11
只改横坐标不太行啊

你y轴也得Normalize的

wxl_123

将XY坐标轴修改成和论文里一样试一下

yjj_123

geo24 发表于 2021-8-18 14:49
你y轴也得Normalize的

兄弟，我改了不行啊，低频的地方应该到-200才对

yjj_123

wxl_123 发表于 2021-8-18 15:53
将XY坐标轴修改成和论文里一样试一下

clc;
D=importdata('D:\Xilinx\project\sn_mash11_v1\DATA3.txt');%载入数据
M1=D-42.6;%average=222.5
x=M1(1:2^21);
%x=M(:,1);%数据第一列数据赋予到V
Fs = 50e6;

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 = 02*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)')
参考的这个代码