非线性规划程序

dx999

#include<iostream.h> & i: @& j* F! v4 ]# }8 z2 a#include<math.h> ( l2 I& ~5 W7 J- w' _double fx(int n,double p[],double x) # H5 m7 g3 }' x1 _, r4 g. W O, o{ + T/ F9 P. ]; O8 W# S! z double z=0,y=1; int i=0,j=0; 0 o7 o; q2 |, R; N% G for(;i<n+1;i++) { 6 [0 L+ V- X( N. \: T* i. u for(;j<i;j++) y*=x; 4 A/ w/ |9 _+ L6 A4 G8 y) J6 |' b$ j z+=*(p+i)*y; } return z; * E& P0 I6 C( ]3 P( |- B}

int fi(int k) # h. S- c2 y* `1 R{ if(k==0) ! J r. C6 U# [0 m% K$ ~ return 1; + h' Z, M% d* a3 Z else M* b! \: j( b4 E' k7 m { if(k==1) ; N: U$ e% J7 t6 ? return 1; 8 r; e4 U9 ]4 n else & U" T: T# @/ i$ q0 Y1 H { 4 c$ i0 A1 L( q int i=0,i1=1,i2=1,w=0; for(;i<k-1;i++) $ Y6 d. X7 _. c4 \ { w=i1+i2; i1=i2; i2=w; - F3 g3 i0 D2 m$ k4 M } ' ?6 g1 T1 r2 p2 T return w; } } 8 S- P! C* A" Q4 L% ^5 D4 A}

3 a1 h% B* t9 h

int r(double w,double a,double b) " ~" _+ F; A# e' S, L% z; L: C0 m4 l{ 4 S* ~* @ [5 e int i=0; , l5 B F' R7 O# c: C: G0 g for(;fi(i)<=(b-a)/w;i++) 5 N+ D# w4 X9 I9 C: n5 o {} 3 m# W9 G6 W9 O G$ ] return i; 8 H+ k; s! d. ~0 e} . Y8 {6 U5 O; V- B6 }# pvoid main() { 2 C" I. u+ P9 C1 I cout<<"输入函数的最高次幂"; ' K* x6 o9 A9 g. ?4 | int n; cin>>n; double *p=new double[n+1]; : e; W- Y( k5 Y4 E7 ^ double a,b,g=0,j=0,c=0; cout<<"依升降幂输入x项的系数"<<endl; 0 @7 h1 b2 f! C/ i0 V3 T for(int i=0;i<n+1;i++) 3 G/ k. p. g( C- O: t! G9 s: u" ~ cin>>*(p+i); / E9 _3 ?/ `) A- D# E2 J1 @& W cout<<"输入误差值"<<endl; double w; cin>>w; . l0 b6 L$ G/ K6 { cout<<"输入初始区间"<<endl; cin>>a>>b; cout<<"输入辨别常数"<<endl; 6 Y- c$ _5 ]1 B/ p# E& u c cin>>c; " d9 {2 p, ^. k n=r(w,a,b); i=0; cout<<fi(n-i-2)<<endl; cout<<a<<" "<<b; , i& f( C9 b; b4 \. _) |# q g=a+fi(n-2-i)*(b-a)/fi(n-i); + ]. J, T5 {8 o j=a+fi(n-1-i)*(b-a)/fi(n-i); 1 F* {, n# X! H. |& k5 k cout<<g<<" "<<j; " p. c5 H$ p; K" l# b# Y for(i=1;i<n-1;i++) ) a/ i: Q- U1 S { # |9 w, l0 m2 H. }- F1 q5 |; h8 C) m if(fx(n,p,g)>fx(n,p,j)) / E% k; q! W; y9 \7 {5 u9 v4 b { $ {; y+ G' T" P9 T' Z* C" U) m$ Y a=g; g=j; ' Z: P/ H; D8 Y. I! y* ^- a j=a+fi(n-1-i)*(b-a)/fi(n-i); * Q0 q- a% ]9 n& H z } " c& i a0 b: m- }9 _ else + I4 r( p+ o w; Y# X { $ o% L3 B; ], x. X$ j. C' v b=j; 1 N# Y: h% M$ o$ ]6 p4 D( Q5 C# U j=g; # C2 e+ D, Y6 V4 M2 Y. P# p g=a+fi(n-2-i)*(b-a)/fi(n-i); 7 ]! \0 O( X9 U- }+ Y } ) h! o, {8 l3 N: t Q% w& ~ } j=g+c; : X; l/ ~9 [. Q d# E6 b if(fx(n,p,g)>fx(n,p,j)) a=g; 1 T7 D7 Q! U" ^! |, T) ~9 O else b=g; " a8 R3 A) J* o* [, ` cout<<"函数值在 ("<<g<<","<<j<<")之间；"; ' O9 \6 Z, g: g } - D5 S7 k1 H2 P7 j V" f

duanzengbin

C语言吗？

Allos

这个程序可读性太差

zhanfeng

如果加上注释会容易看一些！

ilikenba

用的是什么算法？

harleych

没什么用

leonsun2002

什么算法？

zhyi

有点像牛顿法求解非线性方程组

没有解时就找到某个区间

mengfanqi

可以慢慢研究!!!

garr

不是吧，用C语言求规划？汗。lingo不是能解吗？