策划的频谱

1325(30天)
aaa
aaa 2012年4月24日
回答: 弗朗茨Bouchereau 2021年7月29日
你好,
我有困难策划一个正弦波的频率谱。这段代码,我期望的主要频率成分集中在1 /(2 *π),但他们不是。有什么我错过在我的代码吗?
x =(0:0.1:2 *π);
y = sin (x);
z = fft (y);
z = fftshift (z);
N =长度(y);
f = (- N / 2: N / 2 - 1) / N;
情节(f、abs (z))

接受的答案

里克Rosson
里克Rosson 2012年4月24日
请尝试:
% %时间规范:
Fs = 100;%样本每秒
dt = 1 / f;%秒每样
StopTime = 1;%秒
t = (0: dt: StopTime-dt) ';
N =大小(t, 1);
% %正弦波:
Fc = 12;%赫兹
x = cos(2 *π* Fc * t);
% %傅里叶变换:
X = fftshift (fft (X));
% %频率规格:
dF = Fs / N;%赫兹
f = - f / 2: dF: Fs /第二代;%赫兹
% %绘制光谱:
图;
情节(f、abs (X) / N);
包含(的频率(赫兹));
标题(级响应的);
HTH。
瑞克
6个评论
安雅-盖特曼
安雅-盖特曼 2021年2月9日
卡兰,看起来像瑞克选择了100,他猜到了,或者看你的时间样本,发现其中0.01分开。无论你的样本时间步和1 /步骤。如果你的样品不均匀采样时间,您将需要使用一个函数像毁掉使其均匀采样进行分析。

登录置评。

更多的答案(7)

kittipob techawudtipat
kittipob techawudtipat 2016年4月16日
我的这个> > n = 1:2:5;> > x = -3.98 *(罪((π* n) / 2)) / (n *π)+ 32 *(罪((π* n) / 2)) / ((n。^ 3) *(π,^ 3)* (10 ^ 6));> > y = 1.996 + x;> >的情节(y(1:50))指数超过矩阵维度。
我怎样才能解决这个问题。图我以情节不发生谢谢兄弟

弗朗茨Bouchereau
弗朗茨Bouchereau 2021年7月29日
有多种方法可以计算和情节真实功率谱和功率谱密度在MATLAB(当我说‘真’我意味着输出功率谱值对应于实际功率值)。
1)如果你想计算功率谱,而无需指定许多参数和希望函数选择最佳参数,您可以使用 pspectrum 调用的函数没有输出的情节会给你计算功率谱。
2)如果你想计算功率谱和功率谱密度和想要完全控制窗口大小,窗口重叠,窗口类型,和FFT点数,你可以使用韦尔奇周期图 pwelch 函数。调用的函数没有输出的情节会给你计算功率谱。
3)如果你想想象功率谱,可以使用 信号分析仪的应用 。这个程序让你想象你的信号同时在时间,频率,和时频域。你可以放大信号在这些地区的兴趣和分析光谱放大区域。
4)如果你有你的信号分割成多个信号帧可以使用 频谱分析仪的范围
最后,这是一个受欢迎的MATLAB医生页面解释了FFT和真实功率谱之间的关系: 使用FFT功率谱密度估计

检波器博士
检波器博士 2012年4月24日
你的“f”看起来了。您的配置(使用fftshift):
N =长度(y);
dt = x - x (1) (2);%时间增量
《= 1 / (2 * dt);%奈奎斯特频率
df = 1 / (N * dt);%的频率增加
如果国防部(N, 2) = = 0% N是偶数
f =界线:df: Nyq-df;
其他的% N是奇数
f = [(1 * (df: df:界线))(0:df:界线)];
结束
1评论
检波器博士
检波器博士 2012年4月24日
忘记了if语句中的“= = 0”……哎呀!

登录置评。


reddy
reddy 2015年5月21日
你好每一个人
请你情节光谱变量的EXCEL文件。是哪一个?
谢谢。
请提供我代码也。
Phanindra

Deepanshu古普塔
Deepanshu古普塔 2017年11月23日
我试图阴谋的波频谱字符串摘沿着它的长度的一半。我没有看到任何图当我运行它。c = 343;%音速空气% L = 1000;% % w0字符串长度= 200;% % Bn位移距离= 0;%字符串静止,负部分是% x = 0 (0: L);L % %的距离
n = 2;%。%的模式
fn = (n * c) / 2 * L;
ω= 2 *π* fn;%ω%
omegaN = n *ω;% % w是位移
kn = w0 / c;%波不%
t = 1 / fn;%的时间%
一个= ((8 * w0) / n ^ 2 *π^ 2)*罪((n *π)/ 2);% + ve %入射波部分
W =总和((* cos (omegaN * t) + Bn * sin (omegaN * t)) * sin (kn * x));% n = 1:10%位移eq.数目
情节(W t);标题(铰链字符串)包含(位移)ylabel(“时间”)

Khathutshelo Mudau
Khathutshelo Mudau 2020年9月27日
Hy,有人能帮我确定或绘图错误声音信号的频谱

内尔敏哈姆丹
内尔敏哈姆丹 2020年11月24日
它要求0 0 1 0数据通过一个狭窄的通信通道传输。使用MATLAB程序,找到在接收机接收到的信号与带限信道的影响
第一个情节提交0 0 1 0数据在时间上飞机。然后找到并画出频谱。我不知道要做什么你能帮我请

社区寻宝

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

开始狩猎!