建模案例—飞行管理问题要点.ppt

data: alpha=0.000000,5.391190,32.230953,5.091816,20.963361,2.234507, 5.391190,0.000000,4.8040024,6.813460,5.807866,3.815925, 32.230953,4.804024,0.000000,4.364672,22.833654,2.125539, 5.091816,6.613460,4.363673,0.000000,4.537692,2.989819, 20.963361,5.807866,22.833654,4.537692,0.000000,2.309841, 2.234507,3.815925,2.125539,2.989819,2.309841,0.000000; beta=0.000000 109.263642 -128.250000 24.179830 173.065051 13.474934 109.263642 0.000000 -88.871096 -42.243563 -92.304847 9.000000 -128.250000 -88.87096 0.000000 12.476311 -58.786243 0.310809 24.179830 -42.243563 12.476311 0.000000 5.969234 -3.525606 174.065051 -92.304846 -58.786244 5.969234 0.000000 1.914383 14.474934 9.000000 0.310809 -3.525606 1.913383 0.000000 ; enddata end 用MATLAB计算编程如下 function [f,g]=plane(x) alph=[0.000000,5.391190,32.230953,5.091816,20.963361,2.234507, 5.391190,0.000000,4.8040024,6.813460,5.807866,3.815925, 32.230953,4.804024,0.000000,4.364672,22.833654,2.125539, 5.091816,6.613460,4.363673,0.000000,4.537692,2.989819, 20.963361,5.807866,22.833654,4.537692,0.000000,2.309841, 2.234507,3.815925,2.125539,2.989819,2.309841,0.000000]; bet=[0.000000 109.263642 -128.250000 24.179830 173.065051 13.474934 109.263642 0.000000 -88.871096 -42.243563 -92.304847 9.000000 -128.250000 -88.87096 0.000000 12.476311 -58.786243 0.310809 24.179830 -42.243563 12.476311 0.000000 5.969234 -3.525606 174.065051 -92.304846 -58.786244 5.969234 0.000000 1.914383 14.474934 9.000000 0.310809 -3.525606 1.913383 0.000000]; f=abs(x(1))+abs(x(2))+abs(x(3))+abs(x(4))+abs(x(5))+abs(x(6)); g(1)= alpha(1,2) -abs(beta(1,2)+0.5*x(1)+0.5*x(2)); g(2)= alpha(1,3) -abs(beta(1,3)+0.5*x(1)+0.5*x(3)); g(3)= alpha(1,4) -abs(beta(1,4)+0.5*x(1)+0.5*x(4)); g(4)= alpha(1,5) -abs(beta(1,5)+0.5*x(1)+0.5*x(5)); g(5)= alpha(1,6) -abs(beta(1,6)+0.5*x(1)+0.5*x(6)); g(6)= alpha(2,3) -abs(beta(2,3)+0.5*x(2)+0.5*x(3)); g(7)= alpha(2,4) -abs(beta(2,4)+0.5*x(2)+0.5*x(4)); g(8)= alpha(2,5) -abs(beta(2,5)+0.5*x(2)+0.5*x(5)); g(9)= alpha(2,6) -abs(beta(2,6)+