频率之间的转移价值和阴谋

                              lbd1 = 1053;
                             
                              f1 = 3 e8 / (lbd1 * 10 ^ 9)
                             
                                f1 = 2.8490 e + 14

f1和f2 =。

然而,采样频率

sample_rate = 1 e12汽油;

这是两个数量级小于f1和f2,当它应该是超过2倍大。

因此,频率在p别名

aliasfreq = 9.9715 e + 10;

情节的抵消(除以10 e8)之前

                              f1 - aliasfreq
                             
                                ans = 2.8480 e + 14

《旅,我有麻烦了代码和所有的偏移和操作,但我敢肯定这混叠是一个问题。

同时,因为这个问题,N是奇数f向量

f = (((N - 1) / 2): ((N - 1) / 2)) / N * sample_rate + f0;

这将比f在代码中略有不同。

5个评论
显示4年长的评论隐藏4年长的评论

Achraf AYEB 2023年5月3日

编辑:Achraf AYEB 2023年5月3日

清除所有

关闭所有

sample_rate = 1 e13; %采样率

T = 1 / sample_rate;

L = 1 e6;%信号长度

t = (0: L) * t; %时间向量

f0 = 3 e8/1053e-9;

%激光频率赫兹= >λ= C / F =成分

w = 2 *π* f0; %激光脉冲

调频= 3 e9; % GHz的调制频率

wmod = 2 *π*调频;%相位调制

N =长度(t);

lbd1 = 1053;

f1 = 3 e8 / (lbd1 * 10 ^ 9);w1 = 2 *π* f1;

lbd2 = 1053;

f2 = 3 e8 / (lbd2 * 10 ^ 9);w2 = 2 *π* f2;

m = 0; %调制深度

gaussian_order = 0;

g0 = 6.6; % ampli收益

l = 2;%材料长度

转变= 40 e-9; %时间转变

shift2 = 80 e-9; %时间转变

σ= 10 e-9; %高斯HWHM

p = exp (1/2 * ((t-shift) /(σ/ 2))。^ gaussian_order)。* exp(1我。* w1。* t) + exp (1/2 * ((t-shift2) /(σ/ 2)。^ gaussian_order)。* exp(1我。* w2。* t);%广场超高斯脉冲模式

x = p。* exp(1我* * cos (wmod * t));

y = fft (x); % fft supergaussian

y = y / max (abs (y)); %正常化fft

z = fftshift (y); %将傅里叶变换X通过将零频分量阵列的中心

f = (((N - 1) / 2): ((N - 1) / 2)) / N * sample_rate + f0;

%频率向量

n =长度(f);

% f向量的元素数量

%材料荧光

氟= g0 * exp (1/2 * ((3 e8. / f - 1053 - e - 9) / 0.7 e-9)。^ 2);

amp = sqrt (abs (z));

% z的规范

获得= amp。* exp (l *氟);

%的总增益

增益=增益/ max(获得);

yamp = fftshift (z。* exp (l *氟));%将傅里叶变换X通过将零频分量阵列的中心

p_freq = fft (p); % fft的平方脉冲

p_freqshifft = fftshift (p_freq); %将傅里叶变换X通过将零频分量阵列的中心

p_freqshifft = p_freqshifft / max (abs (p_freqshifft)); %正常化

%在ns

t = t * 10 e8;

%在GHz频率

f = f / 10 e8;

%显示图表2 * 2

次要情节(2,1,1)

%显示荧光

情节(f,氟/ max(氟))

%的荧光

传奇(氟)

持有

%冲击谱

情节(f (1: n)、abs (p_freqshifft (1: n)))

%的荧光

传奇(脉冲没有灯)

持有

标题(功率谱)

包含(“频数(GHz) ');

ylabel(“规范化振幅(U.A.)”)

传奇(氟)

%保存所有

次要情节(2,1,2)

情节(f (1: n),获得(1:n))

标题(功率谱)

包含(频率(GHz))

ylabel(“归一化强度(美)”)

保罗 2023年5月3日

编辑:保罗 2023年5月4日

在当前的代码,f1, f2 = 2.849 e14灯头。所以sample_rate需要大于2.849 e14灯头至少两倍。sample_rate 1 e12汽油和1 e13仍然太小。一旦解决了这个问题,还有别的事情需要考虑。

这里的代码,但我sample_rate大约5 x大于f1和f2

                                   % sample_rate = 1 e13; %采样率
                                  
                                   sample_rate = 1.4 e15;
                                  
                                   T = 1 / sample_rate;
                                  
                                   L = 1 e6;%信号长度
                                  
                                   t = (0: L) *;%的时间向量
                                  
                                   f0 = 3 e8/1053e-9;
                                  
                                   %激光频率赫兹= >λ= C / F =成分
                                  
                                   % w = 2π* * f0; %激光脉冲
                                  
                                   调频= 3 e9;% GHz的调制频率
                                  
                                   wmod = 2 *π*调频;%相位调制
                                  
                                   N =长度(t);
                                  
                                   lbd1 = 1053;
                                  
                                   f1 = 3 e8 / (lbd1 * 10 ^ 9);w1 = 2 *π* f1;
                                  
                                   lbd2 = 1053;
                                  
                                   f2 = 3 e8 / (lbd2 * 10 ^ 9);w2 = 2 *π* f2;
                                  
                                   m = 0;%调制深度
                                  
                                   gaussian_order = 0;
                                  
                                   g0 = 6.6;% ampli获得
                                  
                                   l = 2;%材料长度
                                  
                                   转变= 40 e-9;%时间转变
                                  
                                   shift2 = 80 e-9;%时间转变
                                  
                                   σ= 10 e-9;%高斯HWHM
                                  
                                   p = exp (1/2 * ((t-shift) /(σ/ 2))。^ gaussian_order)。* exp(1我。* w1。* t) + exp (1/2 * ((t-shift2) /(σ/ 2)。^ gaussian_order)。* exp(1我。* w2。* t);%广场超高斯脉冲模式
                                  
                                   x = p。* exp(1我* * cos (wmod * t));
                                  
                                   y = fft (x);% fft supergaussian
                                  
                                   y = y / max (abs (y));%正常化fft
                                  
                                   z = fftshift (y);%将傅里叶变换X通过将零频分量阵列的中心
                                  
                                   f = (((N - 1) / 2): ((N - 1) / 2)) / N * sample_rate + f0;
                                  
                                   %频率向量
                                  
                                   n =长度(f);
                                  
                                   % f向量的元素数量
                                  
                                   %材料荧光
                                  
                                   氟= g0 * exp (1/2 * ((3 e8. / f - 1053 - e - 9) / 0.7 e-9)。^ 2);
                                  
                                   amp = sqrt (abs (z));
                                  
                                   % z的规范
                                  
                                   获得= amp。* exp (l *氟);
                                  
                                   %的总增益
                                  
                                   增益=增益/ max(获得);
                                  
                                   % yamp = fftshift (z。* exp (l *氟));%将傅里叶变换X通过将零频分量阵列的中心
                                  
                                   p_freq = fft (p);%的fft平方脉冲
                                  
                                   p_freqshifft = fftshift (p_freq);%将傅里叶变换X通过将零频分量阵列的中心
                                  
                                   p_freqshifft = p_freqshifft / max (abs (p_freqshifft));%正常化

基于defnition p,我们应该期待它的DFT峰值为2.849 e14灯头赫兹,基于f1和f2的值。

freqp = (((N - 1) / 2): ((N - 1) / 2)) / N * sample_rate;

图

情节(freqp、abs (p_freqshifft)),包含(“赫兹”)

xlim ((2.8 - 2.9) * 1 e14灯头)

正是高峰应该是当密谋反对fp,这是正确的频率向量与p_freqshifft有关。

让我们保持一切赫兹

                                   %在ns
                                  
                                   t = t * 10 e8 %;
                                  
                                   %在GHz频率
                                  
                                   % f = f / 10 e8;

现在的情节

图

%显示图表2 * 2

次要情节(2,1,1)

%显示荧光

情节(f,氟/ max(氟))

%的荧光

传奇(“氟”)

持有

当前的情节举行

%冲击谱

情节(f (1: n)、abs (p_freqshifft (1: n)))

%的荧光

传奇(“脉冲没有灯”)

持有

目前公布的情节

标题(的功率谱)

包含(频率(GHz)”);

ylabel (“归一化幅度(U.A.)”)

传奇(“氟”,“p_freqshifft”,“翻译”,“没有”)

%保存所有

次要情节(2,1,2)

情节(f (1: n),获得(1:n))

标题(的功率谱)

包含(“频率(GHz)”)

ylabel (“归一化强度(美)”)

的峰值p_freqshifft似乎是在错误的位置。但这是因为因为它被谋害f,由f0转移(巧合的是,或不等于f1和f2)从freqp,峰值出现在

fpeak = f0 + f1

fpeak = 5.6980 e + 14

图

情节(f (1: n)、abs (p_freqshifft (1: n)))

参照线(fpeak“r”)

xlim ((5.695 - 5.705) * 1 e14灯头)

Achraf AYEB 2023年5月4日

谢谢你保罗,

现在我明白了,我必须为每一个轴曲线。祝你有美好的一天。

登录置评。

类别

信号处理信号处理工具箱数字和模拟滤波器数字滤波器的分析

找到更多的在数字滤波器的分析在帮助中心和文件交换

下载188bet金宝搏

释放

R2022b

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!

频率之间的转移价值和阴谋

0评论
显示1年长的评论藏1年长的评论

接受的答案

5个评论
显示4年长的评论隐藏4年长的评论

更多的答案(0)

另请参阅

类别

标签

下载188bet金宝搏

释放

社区寻宝

频率之间的转移价值和阴谋

0评论 显示1年长的评论藏1年长的评论

接受的答案

5个评论 显示4年长的评论隐藏4年长的评论

更多的答案(0)

另请参阅

类别

标签

下载188bet金宝搏

释放

社区寻宝

0评论
显示1年长的评论藏1年长的评论

5个评论
显示4年长的评论隐藏4年长的评论