非线性规划程序

dx999

#include<iostream.h> #include<math.h> double fx(int n,double p[],double x) { double z=0,y=1; int i=0,j=0; 5 m9 |2 i4 E+ Y4 ? for(;i<n+1;i++) { * U& s K4 T& u5 o2 J# X+ L for(;j<i;j++) ! Y' K. g5 R5 x) W; s/ o' c& o6 N y*=x; 5 a' D" w- d! v' X" w z+=*(p+i)*y; } 2 H& h- b5 N s! F' b return z; ' V4 w4 \+ B, z: J* g}

9 c2 u6 M$ w2 Z6 p* R' P5 r; q8 A

int fi(int k) 8 r |/ l5 D R{ 6 [7 k |) C9 Z, g; \& l if(k==0) return 1; " u! ?- @6 d* B else { 9 q5 c8 |. S( d% R# y x, [ if(k==1) return 1; , z6 B; K6 ~9 V4 |! [$ { else ; s4 A* D% E/ [( [$ M { 6 t+ W# J$ R4 W9 e# a5 M& ^ int i=0,i1=1,i2=1,w=0; for(;i<k-1;i++) - o( Z4 `: K- P1 t! y, w { w=i1+i2; i1=i2; ( g6 k3 R( [% d2 G! H i2=w; . O1 ?& W" i# } } return w; } } }

1 {& M2 Q# _" y! K* L. c

int r(double w,double a,double b) { 2 Q% w. O) [. A6 w+ w9 x int i=0; ( |7 Q( p0 w8 X; J for(;fi(i)<=(b-a)/w;i++) {} return i; } $ W# J& y$ s% z0 ~5 H6 {. pvoid main() { & l' L- x7 y2 ?$ X' g( Z- l4 Y cout<<"输入函数的最高次幂"; ' ?' O- A: Q2 t3 @$ q int n; cin>>n; ' \0 p5 j& o5 L1 n double *p=new double[n+1]; " w" ]2 g* ]* o: t, C6 v5 ? double a,b,g=0,j=0,c=0; cout<<"依升降幂输入x项的系数"<<endl; for(int i=0;i<n+1;i++) & ~% q0 k$ l. m2 O0 |# N" K cin>>*(p+i); cout<<"输入误差值"<<endl; * y3 S: }7 E* F5 k9 b double w; ' t$ G9 ?4 e1 k cin>>w; cout<<"输入初始区间"<<endl; cin>>a>>b; , O6 Y) I: S* a! }: t" \5 m; }3 n cout<<"输入辨别常数"<<endl; cin>>c; n=r(w,a,b); i=0; 5 h" S, x( ` l8 X. R* z! o cout<<fi(n-i-2)<<endl; cout<<a<<" "<<b; / x2 a/ B: H5 t" `% O* n g=a+fi(n-2-i)*(b-a)/fi(n-i); j=a+fi(n-1-i)*(b-a)/fi(n-i); cout<<g<<" "<<j; for(i=1;i<n-1;i++) # C0 N0 G6 c$ `0 i { 4 a: h& ?6 Z) n# t5 V) J 1 D( ?- L4 \. F) i$ n$ j' l if(fx(n,p,g)>fx(n,p,j)) ' f7 d) }. y: D9 U. t: ?- n { # _ w* e- t( l* j- O: W a=g; ! u E' f$ N. r g=j; " k: g% l3 U! O1 f2 ` j=a+fi(n-1-i)*(b-a)/fi(n-i); } 6 U# F2 S; }/ l8 v" a else 5 p6 s; i% _! @. O a { b=j; j=g; / Y2 ?: F+ f7 F4 Z g=a+fi(n-2-i)*(b-a)/fi(n-i); } $ x; ?9 f# j1 K, e9 X9 U- \ } 3 d0 G8 v( ~# b# [" j$ @; n) j2 T5 d j=g+c; % y: p+ z3 E' ]2 d# h& f if(fx(n,p,g)>fx(n,p,j)) " p) F4 ^* I1 `6 C2 N2 \, R) A a=g; else b=g; cout<<"函数值在 ("<<g<<","<<j<<")之间；"; 7 W# Q+ P# N5 G+ _6 D1 q9 e ( ~9 }, B8 [% B& a! ?# G# }}

duanzengbin

C语言吗？

Allos

这个程序可读性太差

zhanfeng

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

ilikenba

用的是什么算法？

harleych

没什么用

leonsun2002

什么算法？

zhyi

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

没有解时就找到某个区间

mengfanqi

可以慢慢研究!!!

garr

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