lise2001 发表于 2008-4-26 23:58

求助如何根据表达式画幅频曲线

我这里有个表达式，就是根据这个表达式画出图形。

而根据相关资料，我找到了参考图。

但是我编程后仿真图与参考图相差很大。

所以想请哪位高手帮忙指点下，哪里出了问题。

我是用.m文件编程的。
公式如下：http://photo.store.qq.com/http_imgload.cgi?/rurl2=e54267b68cfb571d7628cbe959ed8421a844ef2cb926aa91176ae76986f6384075208545ef07e93ad4a9d623f0daa638e0b9d9088233fe01aebd8f522cc9ab3c86c861d141a1f07aa2cf0f30142b9a116ce09e7f

而参考图形如下：http://photo.store.qq.com/http_imgload.cgi?/rurl2=7c7befa4be9cba07dc894aeeed22caefd28767fd25d69bc5eb8d4a726d9ee356f7efa94e9da3a565f3b1266b09383164bce390410bf7668cb719a32cb894b9da818f31e3912b76abc5a5f3362c5a667db6cc12c6

对表达式进行说明：
该模型描述了典型的电力线信道的频率响应，概括了500KHz-20MHz的电力线的信道特征。


我的程序如下：echo on
t0=2.*10.^7;
fs=0.2;
ts=1/fs;
df=0.3;
f=;
m=0.64.*exp(-1.56.*(10.^(-7)).*f).*exp(-2.66.*j.*pi.*f)+0.38.*exp(-1.73.*(10.^(-7)).*f).*exp(-2.97.*j.*pi.*f)-0.15.*exp(-1.9.*(10.^(-7)).*f).*exp(-3.26.*j.*pi.*f)+0.05.*exp(-2.09.*(10.^(-7)).*f).*exp(-3.57.*j.*pi.*f);
plot(f,m(1:length(f)))
xlabel('Frequency')
title('The Discrete MagnitudeSpectrum')


但是这个程序运行后的仿真图与参考图大不一样。。。


所以想请哪位高手帮忙指点下哪里出了问题？

[ 本帖最后由 eight 于 2008-4-28 18:20 编辑 ]

ch_j1985 发表于 2008-4-27 06:35

回复 楼主 的帖子

初到请看置顶帖：
http://forum.vibunion.com/forum/thread-48179-1-1.html
另外，把图片放在附件里面就可以显示啦

lise2001 发表于 2008-4-28 18:01

好的，谢谢楼上的提醒~~~~我会好好去弄下的哦~~

ch_j1985 发表于 2008-4-28 22:28

原帖由 lise2001 于 2008-4-28 18:01 发表 http://www.chinavib.com/forum/images/common/back.gif
好的，谢谢楼上的提醒~~~~我会好好去弄下的哦~~

我看了一下你的程序，画幅频特性时应该不带j的，你最好把图片和公式帖出来，不然代码不好修改:@)

lise2001 发表于 2008-4-29 00:15

回复 4楼 的帖子

你好，很感谢你肯帮我
只是它说发帖时间与现在已经超过30分钟了
所以不能重新编辑
我可能要重新发下帖子。。。。

[ 本帖最后由 eight 于 2008-4-29 09:52 编辑 ]
页: [1]
查看完整版本: 求助如何根据表达式画幅频曲线