!!!高手请进!!!!!关于MATLAB的一个难题!!!!

jxw9987 · 发表于 2006-1-26 06:45:07

...x1.. ...x2... ...x3... .....x4 ... ....x5... 67852.30 45746.4 25994.00 56850.00 17486.40 6448.30 45556.80 28153.40 70125.60 30623.80 69289.30 41981.70 28730.40 111355.40 53119.20 73300.40 40946.90 27247.20 115514.80 63456.20 76617.80 42281.60 31680.00 128503.70 72518.00 78026.90 48169.70 46663.70 148086.30 82137.00 80663.30 46445.40 47652.30 162964.20 81917.20 79604.30 47147.30 47429.60 178951.20 88771.80 89890.60 43253.30 52092.40 219822.30 105751.80 90605.30 51151.50 58215.70 215536.00 114444.70 92534.60 59049.60 65231.30 211229.70 123419.30 83960.30 66947.50 73222.20 206933.30 138917.70 76054.00 60066.80 103450.20 202637.00 153416.00 94199.00 69304.20 115705.80 235220.00 130230.00 83954.00 95386.50 150335.50 228450.00 122907.00 80517.00 133594.70 178667.30 231543.00 137808.00 91311.50 148691.80 208054.70 240532.00 143005.00 Y=[351762.00 384302.00 383192.00 414541.00 421087.00 448728.00 ... 471299.00 496549.00 533040.00 587577.00 675053.00 ... 832609.00 1081916.00 1221742.00 1408238.00 1365448.00 1323221.00]';

=[ 55037 60249 60599 98093 111875 101280 113808 147206 ... 147602 199188 206813 236232 398587 540154 587452 573621 ... 590536]; 题大意是这样的:反映经济指标的 P 最大与 Y 最小两个同时满足的目标作为线性规划模型目标.这样我用DPS处理后是这样的两个方程: Y=665206.4-3.219x1-9.45x2+11.708x3-0.674x4+1.565x5 p=75041.3+1.38x1-5.580x2+6.815x3-0.344x4-0.121x5 我用MATLAB 就是做不出来,希望高手帮忙啊, 希望你把编写的程序写在下面!!!!!!!!!!!!!!!(要能够看到相关系数啊) 谢谢啦!!!!!!!!!!!!!!

[此贴子已经被作者于2006-1-25 23:00:11编辑过]

jxw9987 · 发表于 2006-1-26 06:46:20

<

>我是个菜鸟,希望高手帮帮忙啊!!!
<

>来北京我请吃烤鸭啊~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

young_ys · 发表于 2006-1-26 19:17:45

<

>其实只要问题解决就可以啦，没有必要找matlab来用的
<

>其他的软件，比如你说得dps也应该不错吧
<

>我看你的问题可以使用下lingo软件，或许轻松点，掌握matlab工具箱的应用还是比较麻烦的

数模雄鹰 · 发表于 2006-1-30 19:21:46

<

>下面是Matlab程序:
<

>X=[67852.30 45746.4 25994.00 56850.00 17486.40 6448.30 45556.80 28153.40 70125.60 30623.80 69289.30 41981.70 28730.40 111355.40 53119.20 73300.40 40946.90 27247.20 115514.80 63456.20 76617.80 42281.60 31680.00 128503.70 72518.00 78026.90 48169.70 46663.70 148086.30 82137.00 80663.30 46445.40 47652.30 162964.20 81917.20 79604.30 47147.30 47429.60 178951.20 88771.80 89890.60 43253.30 52092.40 219822.30 105751.80 90605.30 51151.50 58215.70 215536.00 114444.70 92534.60 59049.60 65231.30 211229.70 123419.30 83960.30 66947.50 73222.20 206933.30 138917.70 76054.00 60066.80 103450.20 202637.00 153416.00 94199.00 69304.20 115705.80 235220.00 130230.00 83954.00 95386.50 150335.50 228450.00 122907.00 80517.00 133594.70 178667.30 231543.00 137808.00 91311.50 148691.80 208054.70 240532.00 143005.00];
<

>Y=[351762.00 384302.00 383192.00 414541.00 421087.00 448728.00 ... 471299.00 496549.00 533040.00 587577.00 675053.00 ... 832609.00 1081916.00 1221742.00 1408238.00 1365448.00 1323221.00]'; [m,n]=size(X); t=ones(m,1); X=[t,X];%第一列置1 P=[ 55037 60249 60599 98093 111875 101280 113808 147206 ... 147602 199188 206813 236232 398587 540154 587452 573621 ... 590536]';
%求Y与X的线性回归 [b1,bint1,r1,rint1,stats1]=regress(Y,X,0.05); fprintf('回归方程为:Y=%7.1f+%7.1f*x1+%7.1f*x2+%7.1f*x3+%7.1f*x4+%7.1f*x5\n',b1(1),b1(2),b1(3),b1(4),b1(5),b1(6)); 
 %求P与X的线性回归 [b2,bint2,r2,rint2,stats2]=regress(P,X,0.05); fprintf('回归方程为

=%7.1f+%7.1f*x1+%7.1f*x2+%7.1f*x3+%7.1f*x4+%7.1f*x5\n',b2(1),b2(2),b2(3),b2(4),b2(5),b2(6));

结果为:
回归方程为:Y=521271.2+ -0.8*x1+ -9.6*x2+ 11.8*x3+ -1.1*x4+ 1.9*x5 回归方程为

=134209.5+ 0.0*x1+ -5.5*x2+ 6.7*x3+ -0.1*x4+ -0.2*x5
该结果与楼主给出的结果有差异,可能是楼主数值输入有误.

[此贴子已经被作者于2006-1-30 11:24:43编辑过]

jxw9987 · 发表于 2006-2-2 11:17:49

<

>>> X=[67852.3 45746.4 25994 56850 17486.4 64448.3 45556.8 28153.4 70125.6 30623.8 69289.3 41981.7 28730.4 111355.4 53119.2 73300.4 40946.9 27247.2 115514.8 63456.2 76617.8 42281.6 31680 128503.7 72518 78026.9 48169.7 46663.7 148086.3 82137 80663.3 46445.4 47652.3 162964.2 81917.2 79604.3 47147.3 47429.6 178951.2 88771.8 89890.6 43253.3 52092.4 219822.3 105751.8 90605.3 51151.5 58215.7 215526 114444.7 92534.6 59049.6 65231.3 211229.7 123419.3 83906.3 66947.5 73222.2 206933.3 138917.7 76054 60066.8 103450.2 202637 153416 94199 69304.2 115705.8 235220 130230 83954 95386.5 150335.5 228450 122907 80517 133594.7 178667.3 231543 137808 91311.5 148691.8 208054.7 240532 143005]; >> Y=[351762 384302 383192 414541 421087 448728 471299 496549 533040 587577 675053 832609 1081916 1221742 1048238 1365448 1323221]'; >> [m,n]=size(X); >> t=ones(m,1); >> X=[t,X]; >> P=[55037 60249 60599 68093 111875 101280 113808 147206 147602 199188 206813 ... 236232 398587 540154 587452 573621 595036]; >> [b1,bint1,r1,rint1,stats1]=regress(Y,X,0.05); >> fprintf('回归方程为:Y=%7.1f+%7.1f*x1+%7.1f*x2+%7.1f*x3+%7.1f*x4+%7.1f*x5\n',b1(1),b1(2),b1(3),b1(4),b1(5),b1(6)); 回归方程为:Y=231824.3+ 2.6*x1+ -5.1*x2+ 8.4*x3+ -2.8*x4+ 4.5*x5 >> [b2,bint2,r2,rint2,stats2]=regress(P,X,0.05); ??? Error using ==> or Matrix dimensions must agree.
<

>Error in ==> regress at 68 wasnan = (isnan(y) | any(isnan(X),2));
<

>>> fprintf('回归方程为

=%7.1f+%7.1f*x1+%7.1f*x2+%7.1f*x3+%7.1f*x4+%7.1f*x5\n',b2(1),b2(2),b2(3),b2(4),b2(5),b2(6)); 回归方程为

=231824.3+ 2.6*x1+ -5.1*x2+ 8.4*x3+ -2.8*x4+ 4.5*x5 >> 
我按你说的编程后不行啊,和你的结果也不一样啊
我用的是MATLAB7.1
请高手指教啊

jxw9987 · 发表于 2006-2-2 11:18:39

<

>在程序上好象还有语法上的错误啊

		自动登录	找回密码
密码			注-册-帐-号