24小时热门版块排行榜    

查看: 584  |  回复: 4
当前主题已经存档。

xwz118

金虫 (正式写手)

[交流] [求助]求助用plotyy语句合成下面的例图

请教各位大侠,一般我们用matlab的plotyy语句把下面图合成一个图,具体要求如下:
   它们X轴相同,范围0.6:0.0001:7
         左边Y轴代表牵引力,单位是N
         右边Y轴代表的是总效率,0-1
两个分程序如下:
x=0.6:0.0001:7;
G=235000;
k=222186.2;
f=0.1;
fkp=(0.999*k*G+0.1094*f*k*G)./(0.1094*k+G*x)-G*f;
plot(x,fkp,'-b*')
hold on
title('行走系统牵引力理论特性曲线')
xlabel('潜孔钻机行走速度(km/h)')
ylabel('牵引力(N)')
grid on



v=0.6:0.0001:7;
delta_p=[8000000 15000000 25000000 35000000];
Cs=0.8e-9;
Csm=2.0e-9;
u=0.051;
n=2200;
Cv=0.2e+6;
Cvm=0.6e+6;
Cf=0.01;
Cfm=0.015;
Tc=5;
Vmax=45e-6;
nc=0.9;
rd=0.3223;
for i=1:1:length(delta_p)
for j=1:1:length(v)
    npv(i,j)=1-(1.07882*Cs)*(60*delta_p(i))/((u*n)*v(j));
    npt(i,j)=1/(1+(Cv*(u*n)/(60*delta_p(i)+Cf+(2*3.141592*Tc)/(delta_p(i)*Vmax))*(1.07882/v(j))));
    np(i,j)=npv(i,j)*npt(i,j);
    nmv(i,j)=1/(1+(Csm*delta_p(i))/(u*(v(j)/(0.377*rd*0.95))));
    nmt(i,j)=1-Cfm-(Cvm*(u*(v(j)/(0.377*rd*0.95))/delta_p(i)));
    nm(i,j)=nmv(i,j)*nmt(i,j);
    n1(i,j)=np(i,j)*nm(i,j)*nc;
end
end
figure(1)
plot(v,n1(1,,'-b')
hold on
plot(v,n1(2,,':r')
plot(v,n1(3,,'-.b')
plot(v,n1(4,,'--k')
title('泵—马达系统的理论效率分析')
xlabel('潜孔钻机行走速度(km/h)')
ylabel('总效率')
legend('8MPa','15MPa','25MPa','35MPa')
grid on
hold off

如何用plotyy语句把它们连起来,得出例题样子的图形

[ Last edited by xwz118 on 2009-3-31 at 22:03 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dnp

荣誉版主 (知名作家)

小木虫浪子

优秀版主

What would Jesus do?
2楼2009-03-30 16:38:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xwz118

金虫 (正式写手)

这个我知道
怎么写呢?
3楼2009-03-30 16:46:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abingchem

木虫 (著名写手)

★ ★
kuhailangyu(金币+1,VIP+0):xiexie 3-30 21:01
xwz118(金币+1,VIP+0):谢谢 3-31 17:13
help plotyy

PLOTYY Graphs with y tick labels on the left and right.
    PLOTYY(X1,Y1,X2,Y2) plots Y1 versus X1 with y-axis labeling
    on the left and plots Y2 versus X2 with y-axis labeling on
    the right.

    PLOTYY(X1,Y1,X2,Y2,FUN) uses the plotting function specified by
    the string FUN instead of PLOT to produce each plot.  FUN
    should be a plotting function like 'plot','semilogx','semilogy',
    'loglog','stem', etc. that accepts the syntax H = FUN(X,Y).

    PLOTYY(X1,Y1,X2,Y2,FUN1,FUN2) uses FUN1(X1,Y1) to plot the data for
    the left axes and FUN2(X1,Y1) to plot the data for the right axes.

    [AX,H1,H2] = PLOTYY(...) returns the handles of the two axes created in
    AX and the handles of the graphics objects from each plot in H1
    and H2. AX(1) is the left axes and AX(2) is the right axes.  

    See also PLOT.
冰是从最寒冷的那天开始融化的
4楼2009-03-30 18:58:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunxiao

荣誉版主 (著名写手)

小木虫钓鱼协会会长


xwz118(金币+1,VIP+0):谢谢 3-31 17:14
给你个例子吧!
x=0:0.1*pi:2*pi
y=sin(x);
z=exp(x);
plotyy(x,y,x,z,'plot','semilogy'))
西南交大隧道方向招硕、博士
5楼2009-03-30 22:43:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xwz118 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 yPv6ndDtiRxQ 2026-08-11 4/200 2026-08-11 21:15 by thgx1vqMqrLw
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 yPv6ndDtiRxQ 2026-08-10 7/350 2026-08-11 20:39 by thgx1vqMqrLw
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 yPv6ndDtiRxQ 2026-08-10 8/400 2026-08-11 19:56 by thgx1vqMqrLw
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[基金申请] FileCode能看出啥? +7 要乐观耀哥 2026-08-10 17/850 2026-08-11 10:50 by 要乐观耀哥
[基金申请] filecode +9 documentary 2026-08-10 9/450 2026-08-11 10:08 by wk7465
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
信息提示
请填处理意见