怎么用matlab画幅频曲线，急啊，求大家帮个忙！
谁有编好的程序啊，小弟谢谢了！！[ 本帖最后由 joker 于 2007-8-10 12:13 编辑 ] ezplot也可以把，你知道了幅频方程，那其中就只有幅、频两个未知数。
回复 #1 joker 的帖子
不知道你画出了没有？我曾用plot直接画过，就是不光滑．用maple画效果差不多，不理想．希望互相交流、学习。:@)回复 #3 mechanic05 的帖子
不光滑哦？为什么不光滑，会不会是你取得步长太大我怎么没遇到过呢 前辈好！好久不见，暑期过的可好？
不知道前辈用什摸方法画的，编写ＦＯＲＴＲＡＮ或者ＭＡＴＬＡＢ程序吗？那样得反解出显式表达式，但是当幅频关系为高次时，是很困难甚至是不可能的．我画出的不光滑，是用ＭＡＰＬＥ和ＭＡＴＬＡＢ（直接用ＰＬＯＴ语句）都这样，我也正为这个郁闷．付上我的ＭＡＴＬＡＢ图，请前辈指点．下面是ＭＡＰＬＥ语句：
restart;omega:=1;with(plots);
> mu:=0.3;
> mu:=0.6;
> G:=0.3;
> G:=1;
> G:=2;
> G:=10;
> Gamma:=5*G/(8*omega);
> Gamma:=5*G/(8*omega);
> Gamma:=5*G/(8*omega);
> Gamma:=5*G/(8*omega);
> alpha:=4;
> alpha:=0.5;
> alpha:=0.2;
> beta:=1;
> Gamma:=(3*alpha-3*alpha+alpha)/(8*omega);
> Q:=beta/(2*omega);
> #Q:=beta/(4*omega);
omega := 1
Warning, the name changecoords has been redefined
[animate, animate3d, animatecurve, arrow, changecoords, complexplot,
complexplot3d, conformal, conformal3d, contourplot,
contourplot3d, coordplot, coordplot3d, cylinderplot,
densityplot, display, display3d, fieldplot, fieldplot3d,
gradplot, gradplot3d, implicitplot, implicitplot3d, inequal,
listcontplot, listcontplot3d, listdensityplot, listplot,
listplot3d, loglogplot, logplot, matrixplot, odeplot, pareto,
pointplot, pointplot3d, polarplot, polygonplot, polygonplot3d,
polyhedra_supported, polyhedraplot, replot, rootlocus,
semilogplot, setoptions, setoptions3d, spacecurve,
sparsematrixplot, sphereplot, surfdata, textplot, textplot3d,
tubeplot]
mu := .3
mu := .6
G := .3
G := 1
G := 2
G := 10
Gamma := .1875000000
Gamma := 5/8
Gamma := 5/4
Gamma := 25/4
alpha := 4
alpha := .5
alpha := .2
beta := 1
Gamma := 1.337500000
Q := 1/2
> Y:=(mu*a+mu*(a^2)-Gamma*a^5)^2+(sigma*a-Gamma*a^3)^2-Q^2;
2 5 2
Y := (.3 a + .6 a- .1875000000 a )
3 2
+ (sigma a - 1.337500000 a )- 1/4
> Y:=(mu*a+mu*(a^2)-Gamma*a^5)^2+(sigma*a-Gamma*a^3)^2-Q^2;
> Y:=(mu*a+mu*(a^2)-Gamma*a^5)^2+(sigma*a-Gamma*a^3)^2-Q^2;
> Y:=(mu*a+mu*(a^2)-Gamma*a^5)^2+(sigma*a-Gamma*a^3)^2-Q^2;
> implicitplot((mu*a+mu*(a^2)-Gamma*a^5)^2+(sigma*a-Gamma*a^3)^2-Q^2=0,sigma=0..10,a=0..5);
2 5 2 3 2
Y := (.3 a + .6 a- 5/8 a )+ (sigma a - 1.337500000 a )
- 1/4
2 5 2 3 2
Y := (.3 a + .6 a- 5/4 a )+ (sigma a - 1.337500000 a )
- 1/4
2 5 2 3 2
Y := (.3 a + .6 a- 25/4 a )+ (sigma a - 1.337500000 a )
- 1/4
implicitplot((mu*a+mu*(a^2)-Gamma*a^5)^2+(sigma*a-Gamma*a^3)^2-Q^2=0,sigma=0..10,a=0..5);
implicitplot((mu*a+mu*(a^2)-Gamma*a^5)^2+(sigma*a-Gamma*a^3)^2-Q^2=0,sigma=0..10,a=0..5);
implicitplot((mu*a+mu*(a^2)-Gamma*a^5)^2+(sigma*a-Gamma*a^3)^2-Q^2=0,sigma=0..10,a=0..5);
附加三幅图．谢谢！ＭＡＴＬＡＢ的ＰＬＯＴ语句更简单，和ＭＡＰＬＥ差不多．错在那里？请指点．谢谢！
回复 #4 咕噜噜 的帖子
不好意思:loveliness: 没加上．回复 #4 咕噜噜 的帖子
可能前辈很忙，麻烦抽空指点一、二，不胜感激！:@) 可以用MATLAB啊。把 频率值间隔取小点就行了。
比如 0.1 0.2 ...
如果还不行就 0.1 0.15 0.2 ...... 直到光滑为止。
回复 #11 pengweicai 的帖子
谢谢！不过我不明白您是用什么方法？是ＭＡＴＬＡＢ编程吗？能就我的问题给示范一下吗？非常感谢！！:@) 再问一句，直接用ＰＬＯＴ语句为什么画出的图很不光滑？请指教。
[ 本帖最后由 xinyuxf 于 2007-9-24 15:47 编辑 ] 无水前辈好！暑假过的可好？
我刚搜索了一下，好象前辈遇到过相同的问题，不知道现在作的怎样？用的什么方法？能给指教一下吗？我用ＰＬＯＴ语句直接画幅频响应图，有什么问题吗？谢谢！
[ 本帖最后由 xinyuxf 于 2007-9-24 15:49 编辑 ]
回复 #2 无水1324 的帖子
前辈好！不知道关于幅频响应曲线的画法研究的怎样了？我的问题是推导出了幅频响应方程（隐式的）．按理说非线性的运动幅频响应曲线应该存在跳跃和多值性．我该怎摸做呢？望指教．谢谢！ 你这个能不能把具体问题贴上来，我是说你的公式？我不是很懂maple，所以看的不是很明白，给出问题的话我可以帮你编写matlab看看在matlab里面不知道幅频响应曲线也可以画图，通常使用传递函数 呀，抱歉！公式给错了．那个是为了看懂人家的程序下载的公式．我的公式应该是见附件．
回复 #16 咕噜噜 的帖子
谢谢前辈！我做的是关于梁的大振幅振动问题．得到的幅频响应方程见上面附件．麻烦给看看．请问，解这个隐式方程并画图，为什摸直接用ＭＡＴＬＡＢ语句画出的图很不光滑并且不出现跳跃和多值？非常感谢！:loveliness: 我觉得用ezplot就行了，我也在做这样的频响曲线阿，就是这样画出来的啊！
页:
[1]
2