| 查看: 281 | 回复: 0 | ||
[求助]
新手求助
|
|
我是新手，想通过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) |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
听说今天filecode变了
已经有34人回复
-
基金中了
已经有6人回复
-
filecode
已经有14人回复
-
固定端突然变了，今天
已经有7人回复
-
filecode
已经有11人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有6人回复
-
一个有机合成实验室都需要哪些设备？
已经有12人回复
-
求各位大神看下
已经有33人回复
-
咨询面上基金
已经有7人回复
50