非线性规划程序

dx999

#include<iostream.h> 7 a3 E1 p4 C$ x6 G#include<math.h> ' B w9 w! b+ e6 ?5 N2 W1 c1 v- G4 Udouble fx(int n,double p[],double x) 2 x; f* z) }6 u/ W+ d{ double z=0,y=1; int i=0,j=0; + `* U. p2 }: t. C% @8 [! i for(;i<n+1;i++) ' v/ ` x ?& x { 2 H; A0 a1 |/ L; e for(;j<i;j++) y*=x; z+=*(p+i)*y; 7 }2 s9 D4 Z# }5 z5 c9 o8 g } 5 K2 F4 E* m5 E( o; D return z; }

int fi(int k) { if(k==0) return 1; ! S# x' a0 g0 I. b) M% E8 e else { 5 I3 z$ D# Y9 v% ^ if(k==1) return 1; 5 w; c* \, f' P else { int i=0,i1=1,i2=1,w=0; for(;i<k-1;i++) ! Z, E; b; G" N5 r( C& L { w=i1+i2; i1=i2; i2=w; 3 g: U/ n/ e! D, | } - F& [. F t0 k5 _ return w; 9 D. q( e$ K: S& \; u } } " h. E. B f/ ]; v# R}

, S5 R3 s' D: q* I% @

int r(double w,double a,double b) / [0 t( k( _+ i2 B& f9 g{ int i=0; for(;fi(i)<=(b-a)/w;i++) 3 a* ^% f0 r6 D4 r0 Y {} return i; } void main() { u' T' k/ `) e5 C# ` cout<<"输入函数的最高次幂"; : M' h7 I% Z% _* ^' q9 {) E. n& w int n; cin>>n; 3 S( k4 j4 B9 J# Z( K9 w double *p=new double[n+1]; 4 s" O# x5 A6 o [ double a,b,g=0,j=0,c=0; cout<<"依升降幂输入x项的系数"<<endl; 1 ]/ P1 L* ~. \1 b9 K5 R9 E for(int i=0;i<n+1;i++) 8 c3 q3 x6 { d+ t1 u/ u cin>>*(p+i); cout<<"输入误差值"<<endl; / u7 ^! `/ L; M% Q# l7 v" T double w; cin>>w; / w7 }0 u& R- D# u& Y' @ cout<<"输入初始区间"<<endl; cin>>a>>b; ' y3 ?* g, d5 n. }" q cout<<"输入辨别常数"<<endl; cin>>c; ^! Z/ e- e2 l' n; V2 y# k* M2 C n=r(w,a,b); i=0; " W# K* r, n% k2 @1 ^, J5 p4 E cout<<fi(n-i-2)<<endl; cout<<a<<" "<<b; 1 O. o& x( p* j' ~# [4 h* W4 b g=a+fi(n-2-i)*(b-a)/fi(n-i); j=a+fi(n-1-i)*(b-a)/fi(n-i); $ _+ N3 } ^& N) J5 H' `( _ cout<<g<<" "<<j; for(i=1;i<n-1;i++) 3 x. P, W5 O. H { % z B+ W, N: T- H! U# }: t& }- f , V5 B* d3 e6 ~9 X) F) a if(fx(n,p,g)>fx(n,p,j)) * m2 g3 d" T* H \6 }( [ { a=g; " s+ K9 ~) }$ T. j, x g=j; j=a+fi(n-1-i)*(b-a)/fi(n-i); . K- P. k$ @+ @# Y1 X, } } else { 0 l! O0 p$ Y9 z8 ^, F4 B6 v b=j; 4 O- l, ?" t9 P$ p j=g; ' t$ l h2 c7 l1 S g=a+fi(n-2-i)*(b-a)/fi(n-i); } , [7 l m* [/ z/ P, A+ Q & Y6 j/ x% I5 _# Q8 C$ b } + ^: v- \/ s Z& K$ Q" c! t j=g+c; ; C# }. |$ n1 [! k8 P! o if(fx(n,p,g)>fx(n,p,j)) a=g; else b=g; cout<<"函数值在 ("<<g<<","<<j<<")之间；"; , }+ Y R9 e) c! E2 m) I$ Q) ?0 S" J0 u}

duanzengbin

C语言吗？

Allos

这个程序可读性太差

zhanfeng

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

ilikenba

用的是什么算法？

harleych

没什么用

leonsun2002

什么算法？

zhyi

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

没有解时就找到某个区间

mengfanqi

可以慢慢研究!!!

garr

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