花如月 发表于 2007-10-13 11:22

回复 #19 unknowno 的帖子

（1）6.5里确实没有匿名函数
（2）似乎ezplot也没法达到你所说的功能

unknowno 发表于 2007-10-13 18:20

1、左边那个图形中表示峰值趋势的曲线我用什么方法可以把它画上去？
这个不是问题的问题该怎么解决，实在无从下手，给个指点就可以了！

2、我想把19楼的问题在其他数学软件版块再发一遍，不知道合适不合适？

花如月 发表于 2007-10-15 16:48

原帖由 unknowno 于 2007-10-13 18:20 发表 http://www.chinavib.com/forum/images/common/back.gif
1、左边那个图形中表示峰值趋势的曲线我用什么方法可以把它画上去？
这个不是问题的问题该怎么解决，实在无从下手，给个指点就可以了！


峰值的趋势，可以考虑绘图的时候用函数句柄。然后通过函数句柄访问绘图数据，求出每条曲线最大值处的xy坐标。然后用plot进行绘制就可以了，如果不平滑可以考虑插值。能想到的也就这个思路，应该还是具有可行性的，自己多动手！

VibrationMaster 发表于 2007-10-15 17:49

为什么不能通过求导找极值方式将虚线解出来? 这样也显得文邹邹一点.

花如月 发表于 2007-10-15 17:52

回复 #23 VibrationMaster 的帖子

我没用过求隐函数的极值的函数，应该也是可行的。方法一定有不少，只要肯去发现

unknowno 发表于 2007-10-15 20:39

回复 #22 花如月 的帖子

呵呵！谢谢，有思路了！

aleckxlm 发表于 2007-11-30 01:13

感谢楼上各位大侠
最近在学动力学, 遇到了相同的问题, 看了贴子很有收藏
谢谢


关于ezplot：
If your function has additional parameters, for example k in myfun:
function z = myfun(x,y,k)
z = x.^k - y.^k - 1;
then you can use an anonymous function to specify that parameter:
ezplot(@(x,y)myfun(x,y,2))

[ 本帖最后由 ChaChing 于 2010-1-17 23:05 编辑 ]

laozhou 发表于 2010-4-20 15:56

回复 18楼 花如月 的帖子

其实不用那么麻烦的，在非线性问题中它叫骨干曲线方程，就是等效固有频率曲线！它是关于振幅的函数，直接画出来就可以了

今天无意中看到的帖子，下面借用花无月的那个程序，把中间的骨干曲线方程画出来了
function exaplot(e)
warning off all
for i=1:length(e)
    hold on;
    ezplot(@(s,A)myfun(s,A,e(i)),)
end
A=0:0.01:7;
s=sqrt(1+3/4*0.04*A.^2);%骨干曲线方程
plot(s,A,'r--');
function z = myfun(s,A,e)
z = (1/sqrt((1-s^2+3/4*0.04*A^2)^2+(2*e*s)^2))-A;

[ 本帖最后由 ChaChing 于 2010-4-20 21:52 编辑 ]
页: 1 [2]
查看完整版本: 这样的曲线图如何绘制？