这样的曲线图如何绘制？

花如月

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

unknowno

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

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

花如月

原帖由 unknowno 于 2007-10-13 18:20 发表

                               
登录/注册后可看大图

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

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

VibrationMaster

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

花如月

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

unknowno

呵呵！谢谢，有思路了！

aleckxlm

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


关于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

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

今天无意中看到的帖子，下面借用花无月的那个程序，把中间的骨干曲线方程画出来了
function exaplot(e)
warning off all
for i=1:length(e)
    hold on;
    ezplot(@(s,A)myfun(s,A,e(i)),[0,3,0,7])
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 编辑 ]