24小时热门版块排行榜    

查看: 207  |  回复: 0

qiaobo

金虫 (小有名气)

[求助] 新手求助

我是新手,想通过Matlab出张图,可总是计算不出来。大家帮我看看这个程序有什么问题吗?怎么才能运算出来我要的图呢?
clc
data=[23.09992        0        2.63E-03  %T P V
23.1382        10        -1.18E-03
23.38315        20        -4.84E-03
23.62804        30        -8.54E-03
23.77342        40        -0.0119876
23.84993        50        -1.52E-02
23.94173        60        -1.84E-02
24.00293        70        -2.15E-02
24.01057        80        -2.44E-02
28.26291        0        5.52E-03
28.33147        10        1.58E-03
28.81895        20        -2.24E-03
29.10829        30        -5.93E-03
29.33666        40        -9.56E-03
29.36711        50        -1.29E-02
29.44321        60        -0.0163378
29.55737        70        -1.94E-02
29.51932        80        -2.25E-02
33.41611        0        7.72E-03
33.56775        10        3.79E-03
33.99977        20        5.02E-05
34.32558        30        -3.80E-03
34.43921        40        -7.35E-03
34.49223        50        -0.0108851
34.56798        60        -1.41E-02
34.65128        70        -1.75E-02
34.74216        80        -2.05E-02
38.10681        0        1.02E-02
38.35591        10        6.26E-03
38.84641        20        0.0024524
39.21603        30        -1.37E-03
39.46488        40        -5.04E-03
39.44981        50        -8.66E-03
39.51012        60        -1.21E-02
39.6006        70        -1.54E-02
39.6383        80        -1.85E-02
43.18383        0        1.35E-02
43.2139        10        9.33E-03
43.74744        20        5.25E-03
43.95779        30        1.57E-03
44.21315        40        -2.24E-03
44.33331        50        -5.81E-03
44.40841        60        -9.21E-03
44.44595        70        -1.26E-02
44.491        80        -1.58E-02
48.16497        0        1.65E-02
48.35963        10        1.20E-02
48.62913        20        7.87E-03
48.77882        30        3.71E-03
48.83869        40        -1.64E-04
48.98836        50        -4.10E-03
48.9734        60        -7.41E-03
49.00333        70        -1.07E-02
49.09313        80        -1.40E-02
52.68051        0        0.019018
52.85966        10        1.46E-02
53.46416        20        1.03E-02
53.8745        30        6.30E-03
54.23253        40        2.38E-03
54.38169        50        -1.37E-03
54.4861        60        -5.11E-03
54.59795        70        -8.42E-03
54.65015        80        -1.16E-02
58.25566        0        2.20E-02
58.4045        10        0.0173564
58.82123        20        1.29E-02
58.97004        30        8.66E-03
59.10396        40        4.55E-03
59.223        50        8.98E-04
59.31971        60        -2.76E-03
59.34203        70        -6.29E-03
59.42386        80        -9.65E-03
63.47479        0        0.0248829
63.60104        10        2.02E-02
64.06893        20        0.0157432
64.40309        30        1.14E-02
64.65555        40        7.28E-03
64.87827        50        3.36E-03
64.79662        60        -3.32E-04
64.78176        70        -3.89E-03
64.82632        80        -7.45E-03
68.31429        0        0.0278794
68.49234        10        2.31E-02
68.98933        20        1.86E-02
69.33794        30        1.42E-02
69.61237        40        9.92E-03
69.74586        50        5.97E-03
69.7978        60        2.10E-03
69.82745        70        -1.37E-03
69.90163        80        -4.89E-03
73.40152        0        3.09E-02
73.55723        10        2.60E-02
74.11326        20        2.13E-02
74.49137        30        1.68E-02
74.77309        40        1.25E-02
74.87688        50        8.46E-03
75.01033        60        4.51E-03
75.05481        70        8.67E-04
75.06222        80        -2.75E-03
78.20571        0        3.41E-02
78.28728        10        2.90E-02
78.86563        20        2.42E-02
79.32535        30        1.96E-02
79.55522        40        1.53E-02
79.7851        50        1.11E-02
79.93341        60        7.19E-03
80.00015        70        3.32E-03
80.1188        80        -2.32E-04
];
pp=1.279;%1.032177183;


V=1/pp+data(:,3);T=data(:,1);P=data(:,2);

x=optimresults.x;
c0=0.0894;b3=0.94685;b4=0.000569488;
b1=x(1);b2=x(2);%;

V2=(1-c0.*log(1+P./(b1.*exp(b2.*T)))).*(b3+b4*T);
[V V2];


scatter3(T,P,V,'filled')
hold on
x1fit = min(T):1:max(T);
x2fit = min(P):1:max(P);
[X1FIT,X2FIT] = meshgrid(x1fit,x2fit);
YFIT =(1-c0.*log(1+X2FIT./(b1.*exp(b2.*X1FIT)))).*(b3+b4*X1FIT);
mesh(X1FIT,X2FIT,YFIT)
xlabel('Temp')
ylabel('press')
zlabel('volu')
%surf(xx,yy,zz)
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qiaobo 的主题更新
信息提示
请填处理意见