这样的曲线图如何绘制？

unknowno

这样的图形如何绘制，尤其是象图中那样将一系列曲线画在一个图形中。还有那个θ-s图如何通过方程绘制！
方程为：
A/B=1/((1-s^2+3/4*epsilon*A^2)^2+(2*zeta*s)^2)^(1/2);
theta=arctan(2*zeta*s/(1-s^2+3/4*epsilon*A^2));
方程和示例图形如下, 右边的图形需要结合两个方程进行绘制！

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

花如月

直接绘制有问题么？还有啊，右边的图为什么方程里没有B,但是图里确有B=0.02呢？

[ 本帖最后由 花如月 于 2007-10-12 16:46 编辑 ]

unknowno

左边那个可以直接绘制，右边那个需要结合两个方程进行绘制！

还有就是我不知道如何将一系列的曲线画在同一个图形中！

eight

hold on
请仔细阅读本版置顶贴，里面就有介绍绘图的基本命令

[ 本帖最后由 ChaChing 于 2010-4-20 23:13 编辑 ]

unknowno

哦！我去看看！ 我查了一些matlab书，都没找到！

不好意思！ 看了还是不会！:@L
希望有人指点一下！

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

花如月

首先要相信自己。2个图情形差不多，我把左边的画了下。由于用得是隐函数绘图，因此程序标注比较麻烦（要去掉ezplot自带的label，并在特定的位置加text。还没想到实现的方法，知道的高手望不吝赐教呀）。你先自己在figure窗口里标注吧

1. >> e=[0 0.05 0.1 0.25 0.5 1];
   
2. >> exaplot(e)
   
3. >>

复制代码

结果见附图，exaplot.m的内容如下：

1. function exaplot(e)
   
2. warning off all
   
3. for i=1:length(e)
   
4.     hold on;
   
5.     ezplot(@(s,A)myfun(s,A,e(i)),[0,3,0,7])
   
6. end
   
7. function z = myfun(s,A,e)
   
8. z = (1/sqrt((1-s^2+3/4*0.04*A^2)^2+(2*e*s)^2))-A;

复制代码

[ 本帖最后由 ChaChing 于 2010-4-20 23:14 编辑 ]

unknowno

呵呵！我也是用ezplot，但没想到使用循环来画！
再弱弱的问一下：
@(s,A)这个代表什么意思？:@$
查帮助没查到！


另外你说的右边那个更好画我就不理解了，右边那个需要结合左边那个方程共同来画，怎么会更好画呢？

[ 本帖最后由 unknowno 于 2007-10-12 20:14 编辑 ]

无水1324

去掉label直接用
xlabel('')
ylabel('')
title('')
就可以把

花如月

应该可以的，去掉label问题不大。不过因为是隐函数，所以在特定位置加text不太好弄。

to 楼主：由于你事先没说2个图是相关的。我画第二图的时候发现波形不对，检查程序没有问题，再仔细看才发现是相关的。这样以来第二个图确实就比较麻烦了，需要先求方程（根据参数B求A然后代入）。方程是解出来了，不过解太长了，后边基本没法进行了。

3点建议供参考：
（1）是否可以根据A参数画theta的波形？
（2）考虑换个软件吧，我用matlab根据B求A^2,实际上是解个一元3次方程，虽然可以得解。
        但是化简后还是长得要命。
（3）再等等吧，我的水平也就只能帮你到这一步了，等待真正的高手的到来吧

至于你8楼的那个@的意思，在help里应该是可以搜索到的。matlab基础书上也应该有，或者百度一下匿名函数试试看

[ 本帖最后由 花如月 于 2007-10-12 21:16 编辑 ]

unknowno

谢谢你！:handshake
听了你的建议使我受益非浅，我同时也在考虑maple这个软件能不能解决！

呵呵！
想出解决的办法了，将右边方程中的A^2解出带入左边方程，就得到关于s和θ的方程，再使用6楼的方法就可以很容易画出了！
谢谢你！月月:lol

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

花如月

这个方法我昨天都想过了，似乎中间会有些麻烦。如果你画出来了，希望可以共享一下。:handshake

unknowno

还有一个问题就是，左边那个图形中表示峰值趋势的曲线我用什么方法可以把它画上去？

花如月

那个都不是问题，你先把右边的那个画法帖上来

而且在解A^2过程中会出现复数根，因为每次都可以解出3组。但是也有3个全是复数的情形，请问出现复系数下边怎么做呢？期待楼主的答案哦
>> temp(5,118:120,:)
ans(:,:,1) =
  17.8955 - 0.3057i  18.1525 - 0.3266i  18.4109 - 0.3460i
ans(:,:,2) =
   8.1974 +44.8221i   8.2733 +45.4772i   8.3536 +46.1347i
ans(:,:,3) =
   9.7258 -44.5164i   9.9062 -45.1506i  10.0835 -45.7887i
上边的3组解就都是复数，也就是你说的A^2:@D

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

unknowno

我使用的是matlab6.5，确实查不到@的用法，而且似乎不支持@这个函数，所以就做了如下的变动。


function exaplot(e)
syms x y  %x=s,y=theta
for i=1:length(e)
    hold on;
    ezplot(myfun(x,y,e(i)),[0.75,2,0,pi])
end
function z = myfun(x,y,e)
epsilon=0.2;xi=0.03;
z=2/3*3^(1/2)*((2*xi*x*cot(y)-1+x^2)/epsilon)^(1/2)/e-1/((cot(y)^2+1)^(1/2))/2/xi/x;



>> e=[0.02 0.06 0.1 0.15 0.2 0.25];
>> exaplot(e)

unknowno

你没看清楚我12楼的方法

你所解的是左边方程(A/B)中的A^2（这个方程确实很难处理）
我解的是右边方程(θ)中的A^2，然后再带入左边的方程(A/B)。（右边这个方程就好处理多了）
源程序我已经贴出！

我还有个问题：
能画出第二个图形的突破口，在于第二个方程相对第一个方程要简单许多。
问题是，如果第二个方程和第一个方程一样不能求出A^2的表达式，该怎么办？
换句话说，有没有直接根据两个隐函数来做出变量的曲线的方法。
（不知道叙述的清楚不清楚）

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