matlab fsolve函数原理,求助 关于matlab中fsolve函数的问题

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

function F=myfun(t)

F(1)=9.912.*t(1)-19.84.*t(2)+9.912.*t(3)+0.016.*t(14)-(-0.0867.*t(1).^4-0.2258.*t(2).^4-0.1675.*t(3).^4-0.0812.*t(4).^4-0.0385.*t(5).^4-0.0200.*t(6).^4-0.0112.*t(7).^4-0.0068.*t(8).^4-0.0043.*t(9).^4-0.0030.*t(10).^4-0.0022.*t(11).^4-0.0014.*t(12).^4-0.1617.*t(13).^4-0.4825.*t(14).^4-0.2367.*t(15).^4-0.0851.*t(16).^4-0.0396.*t(17).^4-0.0204.*t(18).^4-0.0114.*t(19).^4-0.0068.*t(20).^4-0.0043.*t(21).^4-0.0029.*t(22).^4-0.0020.*t(23).^4-0.0034.*t(24).^4).*0.0018;

F(2)=9.912.*t(2)-19.84.*t(3)+9.912.*t(4)+0.016.*t(15)-(-0.0483.*t(1).^4-0.1675.*t(2).^4-0.2528.*t(3).^4-0.1756.*t(4).^4-0.0843.*t(5).^4-0.0397.*t(6).^4-0.0204.*t(7).^4-0.0114.*t(8).^4-0.0068.*t(9).^4-0.0043.*t(10).^4-0.0030.*t(11).^4-0.0019.*t(12).^4-0.0779.*t(13).^4-0.2376.*t(14).^4-0.5052.*t(15).^4-0.2458.*t(16).^4-0.0887.*t(17).^4-0.0412.*t(18).^4-0.0211.*t(19).^4-0.0117.*t(20).^4-0.0069.*t(21).^4-0.0043.*t(22).^4-0.0029.*t(23).^4-0.0046.*t(24).^4).*0.0018;

F(3)=9.912.*t(3)-19.84.*t(4)+9.912.*t(5)+0.016.*t(16)-(-0.0246.*t(1).^4-0.0812.*t(2).^4-0.1756.*t(3).^4-0.2552.*t(4).^4-0.1764.*t(5).^4-0.0845.*t(6).^4-0.0397.*t(7).^4-0.0204.*t(8).^4-0.0114.*t(9).^4-0.0068.*t(10).^4-0.0043.*t(11).^4-0.0025.*t(12).^4-0.0434.*t(13).^4-0.0853.*t(14).^4-0.2458.*t(15).^4-0.5082.*t(16).^4-0.2469.*t(17).^4-0.0891.*t(18).^4-0.0413.*t(19).^4-0.0211.*t(20).^4-0.0117.*t(21).^4-0.0069.*t(22).^4-0.0043.*t(23).^4-0.0061.*t(24).^4).*0.0018;

F(4)=9.912.*t(4)-19.84.*t(5)+9.912.*t(6)+0.016.*t(17)-(-0.0136.*t(1).^4-0.0385.*t(2).^4-0.0843.*t(3).^4-0.1764.*t(4).^4-0.2554.*t(5).^4-0.1765.*t(6).^4-0.0845.*t(7).^4-0.0397.*t(8).^4-0.0204.*t(9).^4-0.0114.*t(10).^4-0.0068.*t(11).^4-0.0035.*t(12).^4-0.0261.*t(13).^4-0.0397.*t(14).^4-0.0887.*t(15).^4-0.2469.*t(16).^4-0.5086.*t(17).^4-0.2471.*t(18).^4-0.0892.*t(19).^4-0.0414.*t(20).^4-0.0211.*t(21).^4-0.0117.*t(22).^4-0.0068.*t(23).^4-0.0082.*t(24).^4).*0.0018;

F(5)=9.912.*t(5)-19.84.*t(6)+9.912.*t(7)+0.016.*t(18)-(-0.0082.*t(1).^4-0.0200.*t(2).^4-0.0397.*t(3).^4-0.0845.*t(4).^4-0.1765.*t(5).^4-0.2554.*t(6).^4-0.1765.*t(7).^4-0.0845.*t(8).^4-0.0397.*t(9).^4-0.0204.*t(10).^4-0.0112.*t(11).^4-0.0052.*t(12).^4-0.0168.*t(13).^4-0.0204.*t(14).^4-0.0412.*t(15).^4-0.0891.*t(16).^4-0.2471.*t(17).^4-0.5086.*t(18).^4-0.2471.*t(19).^4-0.0892.*t(20).^4-0.0414.*t(21).^4-0.0211.*t(22).^4-0.0113.*t(23).^4-0.0115.*t(24).^4).*0.0018;

F(6)=t(7)-t(6);

F(7)=t(8)-t(5);

F(8)=t(9)-t(4);

F(9)=t(10)-t(3);

F(10)=t(11)-t(2);

F(11)=t(12)-t(1);

F(12)=7.1.*t(14)-7.145.*t(13)+0.02.*t(1)+7.475-(-0.3339.*t(1).^4-0.1717.*t(2).^4-0.0827.*t(3).^4-0.0460.*t(4).^4-0.0277.*t(5).^4-0.0178.*t(6).^4-0.0122.*t(7).^4-0.0087.*t(8).^4-0.0064.*t(9).^4-0.0049.*t(10).^4-0.0037.*t(11).^4-0.0019.*t(12).^4-0.1654.*t(13).^4-0.1189.*t(14).^4-0.0671.*t(15).^4-0.0379.*t(16).^4-0.0228.*t(17).^4-0.0145.*t(18).^4-0.0098.*t(19).^4-0.0069.*t(20).^4-0.0050.*t(21).^4-0.0037.*t(22).^4-0.0026.*t(23).^4-0.0025.*t(24).^4).*0.0011;

F(13)=7.1.*t(13)-14.258.*t(14)+7.1.*t(15)+0.016.*t(2)+12.558-(-0.1310.*t(1).^4-0.4825.*t(2).^4-0.2376.*t(3).^4-0.0853.*t(4).^4-0.0397.*t(5).^4-0.0204.*t(6).^4-0.0113.*t(7).^4-0.0068.*t(8).^4-0.0043.*t(9).^4-0.0029.*t(10).^4-0.0020.*t(11).^4-0.0012.*t(12).^4-0.1120.*t(13).^4-0.2268.*t(14).^4-0.1656.*t(15).^4-0.0798.*t(16).^4-0.0375.*t(17).^4-0.0193.*t(18).^4-0.0107.*t(19).^4-0.0063.*t(20).^4-0.0040.*t(21).^4-0.0026.*t(22).^4-0.0017.*t(23).^4-0.0025.*t(24).^4).*0.0018;

F(14)=7.1.*t(14)-14.258.*t(15)+7.1.*t(16)+0.016.*t(3)+12.558-(-0.0482.*t(1).^4-0.2367.*t(2).^4-0.5052.*t(3).^4-0.2458.*t(4).^4-0.0887.*t(5).^4-0.0412.*t(6).^4-0.0211.*t(7).^4-0.0117.*t(8).^4-0.0069.*t(9).^4-0.0043.*t(10).^4-0.0029.*t(11).^4-0.0016.*t(12).^4-0.0632.*t(13).^4-0.1656.*t(14).^4-0.2513.*t(15).^4-0.1747.*t(16).^4-0.0837.*t(17).^4-0.0393.*t(18).^4-0.0202.*t(19).^4-0.0111.*t(20).^4-0.0066.*t(21).^4-0.0041.*t(22).^4-0.0026.*t(23).^4-0.0035.*t(24).^4).*0.0018;

F(15)=7.1.*t(15)-14.258.*t(16)+7.1.*t(17)+0.016.*t(4)+12.558-(-0.0238.*t(1).^4-0.0851.*t(2).^4-0.2458.*t(3).^4-0.5082.*t(4).^4-0.2469.*t(5).^4-0.0891.*t(6).^4-0.0414.*t(7).^4-0.0211.*t(8).^4-0.0117.*t(9).^4-0.0069.*t(10).^4-0.0043.*t(11).^4-0.0022.*t(12).^4-0.0357.*t(13).^4-0.0798.*t(14).^4-0.1747.*t(15).^4-0.2547.*t(16).^4-0.1761.*t(17).^4-0.0843.*t(18).^4-0.0396.*t(19).^4-0.0203.*t(20).^4-0.0112.*t(21).^4-0.0066.*t(22).^4-0.0040.*t(23).^4-0.0047.*t(24).^4).*0.0018;

F(16)=7.1.*t(16)-14.258.*t(17)+7.1.*t(18)+0.016.*t(5)+12.558-(-0.0130.*t(1).^4-0.0396.*t(2).^4-0.0887.*t(3).^4-0.2469.*t(4).^4-0.5086.*t(5).^4-0.2471.*t(6).^4-0.0892.*t(7).^4-0.0414.*t(8).^4-0.0211.*t(9).^4-0.0117.*t(10).^4-0.0068.*t(11).^4-0.0032.*t(12).^4-0.0214.*t(13).^4-0.0375.*t(14).^4-0.0837.*t(15).^4-0.1761.*t(16).^4-0.2552.*t(17).^4-0.1763.*t(18).^4-0.0844.*t(19).^4-0.0396.*t(20).^4-0.0203.*t(21).^4-0.0111.*t(22).^4-0.0063.*t(23).^4-0.0065.*t(24).^4).*0.0018;

F(17)=7.1.*t(17)-14.258.*t(18)+7.1.*t(19)+0.016.*t(6)+12.558-(-0.0076.*t(1).^4-0.0204.*t(2).^4-0.0412.*t(3).^4-0.0891.*t(4).^4-0.2471.*t(5).^4-0.5086.*t(6).^4-0.2471.*t(7).^4-0.0892.*t(8).^4-0.0413.*t(9).^4-0.0211.*t(10).^4-0.0114.*t(11).^4-0.0048.*t(12).^4-0.0137.*t(13).^4-0.0193.*t(14).^4-0.0393.*t(15).^4-0.0843.*t(16).^4-0.1763.*t(17).^4-0.2553.*t(18).^4-0.1764.*t(19).^4-0.0844.*t(20).^4-0.0396.*t(21).^4-0.0202.*t(22).^4-0.0107.*t(23).^4-0.0092.*t(24).^4).*0.0018;

F(18)=t(13)-t(24);

F(19)=t(14)-t(23);

F(20)=t(15)-t(22);

F(21)=t(16)-t(21);

F(22)=t(17)-t(20);

F(23)=t(18)-t(19);

F(24)=(t(25)./777.14)+(t(25)./124)-t(1)+281;

F(25)=t(25)-10.*(299.*12-t(13)-t(14)-t(15)-t(24)-t(23)-t(16)-t(17)-t(18)-t(19)-t(20)-t(21)-t(22)-t(23)-t(24));

[t,fval]=fsolve('myfun',[283,283,283,283,283,283,283,283,283,283,283,283,283,283,283,283,283,283,283,283,283,283,283,283,50],optimset('Display','iter'))

这个程序为什么解出来的解是错的 能帮忙看看吗 求助！！