傅里叶大小和相位谱图形使用fft,错误:数组索引必须是正整数或逻辑值。

3视图(30天)

显示旧的评论

路加福音化 2022年10月21日

0
链接

这个问题直接联系

//www.tatmou.com/matlabcentral/answers/1832618-graphing-fourier-magnitude-and-phase-spectra-using-fft-error-array-indices-must-be-positive-intege

评论道: 保罗 2022年10月22日

我写一个脚本得到一个锯齿波函数的频率谱一段2π

我面临的问题是我不能找到一种方法来在我的代码中使用频域。我试着fft()和傅里叶()但是已经会见了错误。任何提示赞赏!

目标:

策划magitude和阶段的Cn和Dn频率(见附加图片我的意思)

                          fft =[];安= [];cnplot =[];问= [];jnt = [];
                         
                          T = 2π*;
                         
                          = 1;
                         
                          funct1 = @ t / t (t);
                         
                          功能= @ (t)的(大小(funct1));
                         
                          w0 = 2 *π* (1 / T);
                         
                          a0 = (1 / (T)) *积分(funct1, 0, T);
                         
                          为n = 1:1:20变化百分比范围从3,100
                         
                          为t = 1:1:360 * 5
                         
                          anfunct = @ (t) funct1 (t) * cos (n * t * w0);
                         
                          bnfunct = @ (t) funct1 (t) * sin (n * t * w0);
                         
                          一个= (1 / (T / 2)) *积分(anfunct, 0, T);
                         
                          bn = (1 / (T / 2)) *积分(bnfunct, 0, T);
                         
                          cn =(一个+ bn * * bn) ^ (1/2);
                         
                          θ=每股((bn) /);
                         
                          dn = (5 * cn) *θ;
                         
                          英国《金融时报》(t) =一个* cos (n * w0 *函数(t)) + bn *罪(n * w0 *函数(t));
                         
                          x (t) = - cn * cos (n * w0 *函数(t) +θ);
                         
                          d (t) = dn * exp (1 * n * w0 * t);
                         
                          k (w) = fft (x (t));%的错误:数组索引必须是正整数或逻辑值。
                         
                          结束
                         
                          fft = (fft;英国《金融时报》);
                         
                          问=[问;x];
                         
                          jnt = (jnt; d);
                         
                          结束
                         
                          F_A = (fft) + a0总和;
                         
                          C_A =(问)+ a0总和;
                         
                          D_A = (jnt) + a0总和;
                         
                          函数tt = 1:1: t);
                         
                          次要情节(1,1)
                         
                          F_A情节(tt)
                         
                          标题(“信号X (t)源自trigonomtric形式”)
                         
                          包含(“t”)
                         
                          ylabel (“f (t)”)
                         
                          次要情节(1、2)
                         
                          C_A情节(tt)
                         
                          标题(“信号X (t)源自紧凑trigonomtric形式”)
                         
                          包含(“t”)
                         
                          ylabel (“f (t)”)
                         
                          ww = (1:1: w);
                         
                          nn = (1:1: n);
                         
                          %以下任何代码不工作
                         
                          次要情节(3,1,3);
                         
                          情节(w, k (w))
                         
                          标题(“振幅谱”)
                         
                          包含(“w”)
                         
                          ylabel (“k”(w))

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

登录置评。

在回答这个问题。

答案(1)

保罗 2022年10月22日

0
链接

直接链接到这个答案

//www.tatmou.com/matlabcentral/answers/1832618-graphing-fourier-magnitude-and-phase-spectra-using-fft-error-array-indices-must-be-positive-intege answer_1080988

你好卢克,

线的问题,我们看到x (t)是非整数数值,不能用来索引数组fft。此时在代码中,w是没有定义的,因此也会成为一个问题。

它看起来像fft旨在成为一个二维数组;你确定你想指数与单个号码吗?

最后,考虑重命名变量fft以免影子的内置函数 fft 。

                              fft =[];安= [];cnplot =[];问= [];jnt = [];
                             
                              T = 2π*;
                             
                              = 1;
                             
                              funct1 = @ t / t (t);
                             
                              功能= @ (t)的(大小(funct1));
                             
                              w0 = 2 *π* (1 / T);
                             
                              a0 = (1 / (T)) *积分(funct1, 0, T);
                             
                              为n = 1:1:20变化百分比范围从3,100
                             
                              为t = 1:1:360 * 5
                             
                              anfunct = @ (t) funct1 (t) * cos (n * t * w0);
                             
                              bnfunct = @ (t) funct1 (t) * sin (n * t * w0);
                             
                              一个= (1 / (T / 2)) *积分(anfunct, 0, T);
                             
                              bn = (1 / (T / 2)) *积分(bnfunct, 0, T);
                             
                              cn =(一个+ bn * * bn) ^ (1/2);
                             
                              θ=每股((bn) /);
                             
                              dn = (5 * cn) *θ;
                             
                              英国《金融时报》(t) =一个* cos (n * w0 *函数(t)) + bn *罪(n * w0 *函数(t));
                             
                              x (t) = - cn * cos (n * w0 *函数(t) +θ);
                             
                              d (t) = dn * exp (1 * n * w0 * t);
                             
                              x (t)
                             
                              k (w) = fft (x (t));%的错误:数组索引必须是正整数或逻辑值。
                             
                              结束
                             
                              fft = (fft;英国《金融时报》);
                             
                              问=[问;x];
                             
                              jnt = (jnt; d);
                             
                              结束
                             
                                ans = -0.0056

                                 数组索引必须是正整数或逻辑值。

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

保罗 2022年10月22日

我没有仔细看,如果,在代码实际上是想计算,所以恐怕我不能评论。

登录置评。

在回答这个问题。

类别

MATLAB 数学傅里叶分析和过滤

找到更多的在傅里叶分析和过滤在帮助中心和文件交换

下载188bet金宝搏

MATLAB

释放

R2022a

社区寻宝

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

开始狩猎!