24小时热门版块排行榜    

查看: 598  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 15/750 2026-09-01 12:23 by icm639
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +8 ziyangfang 2026-08-25 11/550 2026-09-01 09:51 by zzuzxg
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见