MATLAB -信号的振幅和相位谱

98(30天)

显示旧的评论

兰迪·马什 2017年4月22日

0
链接

这个问题直接联系

//www.tatmou.com/matlabcentral/answers/336664-matlab-amplitude-and-phase-spectrum-of-a-signal

评论道: 大卫Goodmanson 2020年9月13日

剧情后信号的振幅和相位谱 s (t) 在MATLAB

信号是周期性的。

可能有人在这个例子显示如何绘制振幅和相位谱信号吗?

另外,你推荐什么免费的pdf书籍解决问题在MATLAB通信吗?

我在哪里可以找到这种类型的MATLAB解决问题?

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

登录置评。

在回答这个问题。

答案(4)

明星黾 2017年4月22日

1
链接

直接链接到这个答案

//www.tatmou.com/matlabcentral/answers/336664-matlab-amplitude-and-phase-spectrum-of-a-signal answer_264012

你需要的符号数学工具箱(除非你想花几个小时计算它自己)。

解决方案是开发线性方程的两个“屋顶”行(“房子”比喻太容易忽视)。这些都是典型的 “y = m * t + b” 方程。你有你需要的所有信息在图中估计和识别的参数方程。

然后将每一个方程(现在只有在 “t” ),乘以 exp (1 * w * t) 在哪里 ' w ' 弧度的频率,和集成产品对吗 “t” 在这个区域的定义。将从第一个积分 1 来 0 (正斜率线)和第二个 0 来 + 1 (负斜率直线)。然后评估后积分(现在只有在表达 ' w ' ),添加在一起,简化表达式。你应该得到一些条款取消,和其他人应该简化三角函数,使用欧拉恒等式。注意,您可以忽略竖线(你不需要方程),因为它们只是供参考。

这应该给你足够的信息来解决这个问题。如果你有问题,你做过什么,描述它的问题,留下一个清晰的线索在代码中我们可以明白你所做的一切。我们也许能够帮助。

我必须问我的教授如何处理当我第一次遇到这样的问题,所以我明白,他们起初并不明显。你没有其他人后您遇到的问题。

3评论
显示2年长的评论隐藏2年长的评论

明星黾 2017年4月22日

编辑:明星黾 2017年4月23日

我明白你想要什么。我不确定具体的地方给你参考,虽然 118bet金博宝app (链接)是一个很好的起点。我想复制粘贴图书的isbn你感兴趣的一个文本文件,然后做一个在线搜索来找到更多关于书籍。的 alibris (链接)网站是我第一次去哪里购买大多数书籍,尤其是技术书籍。(我没有金融或其他“alibris”的兴趣。)

有很多书可以在数字信号处理(或离散),Proakis是,奥本海姆被另一种流行的选择。(两位作者不同的书有不同的合作者)。我查阅了数字信号处理的书,找不到任何专门处理你发布的问题。

我的研究生课程讲义基本信号处理这些问题是作业。他们没有在课本里,处理硬件电路和滤波器的分析。

你发布的分析解决问题的办法是:

                                  信谊t w
                                 
                                  ltz = 2 + 1 * t;%左部(小于零)
                                 
                                  gtz = 2 - 1 * t;%对部分(大于零)
                                 
                                  fourier_kernel = exp(1我* w * t);%傅里叶变换的内核
                                 
                                  F_ltz = int (ltz * fourier_kernel t 1,0);%傅里叶变换的ltz
                                 
                                  F_gtz = int, gtz * fourier_kernel t 0, + 1);%傅里叶变换的“公司”
                                 
                                  F_signal =收集(F_ltz + F_gtz)
                                 
                                  F_signal =重写(F_signal,“要求”)
                                 
                                  F_signal =简化(F_signal,“步骤”,10)

我解决了几个步骤,这样你就可以看到这个想法发展,尽可能多的代码是如何工作的。(因此故意不高效的代码。)我将让你分析振幅和相位(考虑,解决方案完全是真实的,非常简单)。虽然频率成分扩展负来 +正 ,我将从情节 -10 *π 来 + 10 *π 。这将给你一个好主意你的信号的傅里叶变换是如何运作的。

编辑 - - - - - -

阴谋的代码是:

                                  图(1)
                                 
                                  fplot (F_signal [-10 * pi + 10 *π])
                                 
                                  网格
                                 
                                  标题([“美元”乳胶(潜艇(F_signal,' w ',“ω”))“美元”),“翻译”,“乳胶”)

“振幅”是绝对你的信号的傅里叶变换的价值。定义为“阶段” “量化(图像放大(F_signal),真正的(F_signal))” 。

登录置评。

Afshin Aghayan 2017年7月24日

1
链接

直接链接到这个答案

//www.tatmou.com/matlabcentral/answers/336664-matlab-amplitude-and-phase-spectrum-of-a-signal answer_275290

编辑:Afshin Aghayan 2017年8月2日

看下面的Matlab函数,它可以计算相位谱以及振幅谱与完美的准确性:

//www.tatmou.com/matlabcentral/fileexchange/63965-amplitude-and-phase-spectra-of-a-signal--fourier-transform-

这个程序计算输入信号的振幅和相位谱与可接受的精度特别是相位谱的计算。代码是否为计算三个主要工作振幅和相位谱。首先,它扩展了输入信号到正无穷;因为计算傅里叶变换(FT) (fft在Matlab函数),我们认为我们与无限的信号周期波长,代码创建一个super_signal通过将原始信号本身,直到旁边super_signal的长度大约是1000000的样品,为什么我选择1000000样品吗?实际上,它只是基于尝试和错误! !对大多数信号,我试过,晚餐与1000000样品有最好的输出信号。

其次,计算fft在Matlab可以选择不同的决议,Mathwork文档和帮助使用NFFT = 2 ^ nextpow2(长度(信号金宝搏官方网站)),它肯定是不够的,需要高精度输出。在这里,我选择最适合的解决NFFT = 100000信号。

第三,英国《金融时报》的代码过滤结果通过阈值,这是非常重要的一步!计算相位谱,其结果非常吵闹,因为浮动四舍五入误差,它会导致在计算“反正切”甚至小四舍五入误差产生显著噪音相位谱的结果,对于抑制这种噪声可以定义一个阈值。这意味着如果特定频率幅度小于预定义的阈值(您必须定义它)把零代替它。

这三个步骤有助于显著提高振幅和相位谱的结果。

如果你在研究中使用这个程序,请引用以下文章:

Afshin Aghayan, Priyank贾斯瓦尔,哈米德Reza Siahkoohi (2016)。“地震去噪使用冗余提升计划”。GEOPHYSICS, 81(3), V249-V260. https://doi.org/10.1190/geo2015-0601.1

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

登录置评。

图像分析 2017年4月22日

0
链接

直接链接到这个答案

//www.tatmou.com/matlabcentral/answers/336664-matlab-amplitude-and-phase-spectrum-of-a-signal answer_264020

试着使信号的fft。这是一个开始:

                             s = [1 . 01: 2, 1.99: -0.01: 1)
                            
                             t = linspace(1, - 1,长度(s));
                            
                             情节(t, s,“线宽”2);
                            
                             网格在;
                            
                             ylim ([0, 3]);
                            
                             s_spectrum = fft (s);

4评论
显示3年长的评论隐藏3年长的评论

大卫Goodmanson 2020年9月13日

你好兰迪,

s (t)的情节有点误导。由于信号是周期,你可以暂时考虑其中的几个“房屋”加入。第一个房子的终点是1。第二个房子的起点是1。这意味着垂直落差从1到0(和垂直上升,从0到1)实际上并没有发生。或者,你可以说他们发生但彼此抵消这房子连在一起的。所以你最终得到一系列屋顶线,来回从1到2和2比1。找到谱你将只需要一个房子,但起始和结束点是1,没有降至0。

登录置评。