QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 11549|回复: 23
打印 上一主题 下一主题

非线性规划程序

[复制链接]
字体大小: 正常 放大
dx999        

1

主题

0

听众

19

积分

升级  14.74%

该用户从未签到

新人进步奖

跳转到指定楼层
1#
发表于 2004-12-9 20:04 |只看该作者 |倒序浏览
|招呼Ta 关注Ta

#include<iostream.h>3 I3 u8 r/ R/ ~# j: L# m% G #include<math.h> 3 t8 [6 w7 g# V1 gdouble fx(int n,double p[],double x) H. l" Z* H0 C) ~, D; _{6 n0 `; w) g- t+ @ double z=0,y=1;/ x, Z& V% l0 E8 Z' Y: E int i=0,j=0; 8 u3 l4 Q! c$ Y! s for(;i<n+1;i++)0 d- }8 F, j' s2 V: ] { 7 o) H l" ? t4 k! ? for(;j<i;j++) 7 F6 J1 f# ~% k2 G; l y*=x; 3 L! |8 r m# ~8 ]1 A3 z/ k$ J" [ z+=*(p+i)*y; 4 S5 u5 A: X" H3 N2 c d E/ z9 \ } + \; U+ d* [( l7 L7 p9 B: ^6 g return z;2 y a4 S' y) A" C }

/ i6 }/ u) V( m

int fi(int k) 6 B6 a3 f- }: o{4 s+ {( y2 {# b+ H if(k==0) 5 D ]$ f @' ?# h* M6 G8 U return 1;: m4 r: C0 s1 p3 O7 n else : t, H8 e8 x$ }5 q7 t% @ { / I# [% B# Q8 t. F# n0 v8 Z7 V% R if(k==1) & n( h7 s9 R) Y1 z$ c return 1;% `# p, S8 Y* U) q else 0 b; X. U5 u2 r; l j/ P { # o0 r9 I, w6 x; v: S( i" o3 ` int i=0,i1=1,i2=1,w=0;4 H7 e6 B% ^8 `6 @8 D- g for(;i<k-1;i++) j6 a8 E2 ^3 h. F9 ^4 B {" P4 `) F% c/ b8 s3 B K w=i1+i2; L2 G$ e2 [0 _% C e i1=i2; l. o# i& ~9 U/ _& S) S* X% P5 n i2=w; + O, b1 R, C/ a } 8 ^$ S$ r4 e0 P5 n0 ^ return w;/ l- I3 }2 e {* @; k9 O }: y% n4 a/ e# e* r }/ \1 k2 g0 h& O2 | }

) o+ X3 n+ |( b5 i5 m9 A

int r(double w,double a,double b) N/ P" m* B8 {) g{* S, D. A- P8 @% y9 z! ~ int i=0;: T' }7 T% I. I for(;fi(i)<=(b-a)/w;i++)9 E* r* w9 ?8 z/ ^! J2 C& Q3 N% [ {}. [7 G( ^: R$ o9 ^7 Q! b5 D. T$ R return i; 6 z7 m+ o( _' K* h}: x: v; u' {2 {7 U7 x) ]; i void main() . S2 ?5 O; {! S! t t R4 b{: v; G) K8 C" l/ K4 d- ?' C cout<<"输入函数的最高次幂";& v" S1 X) S5 L/ T, O int n;7 [2 B) n$ v" t( K cin>>n; ?; Y7 s! s7 L- ] double *p=new double[n+1];# C2 }* k6 T+ A5 M( d; ^9 I8 R- u double a,b,g=0,j=0,c=0; & r( y+ F% P4 n7 ^7 r1 M9 c2 Z. |/ j cout<<"依升降幂输入x项的系数"<<endl;* Y- ?+ `$ H* w$ i/ M' w for(int i=0;i<n+1;i++) ) |+ ~- p: s" Q1 b( F6 Q6 u& s cin>>*(p+i);0 s5 r# o/ x* P5 I cout<<"输入误差值"<<endl;: F' S7 I7 Z* n2 l" f7 B6 o double w; . x9 c! n% i- D1 n cin>>w;5 ]" y0 H& ?& A9 H6 s( ` cout<<"输入初始区间"<<endl;1 f/ ^8 Q( X& w, t cin>>a>>b; Z: X' T7 H7 h- c cout<<"输入辨别常数"<<endl; ; ?5 O! \8 a6 G) F. \( X3 i cin>>c;9 C( L( g# d. i5 c. y, W n=r(w,a,b);" A4 t1 H# I8 S- p* H2 G. Z/ X i=0; 1 c$ s7 U3 W, [( P6 U. o. R cout<<fi(n-i-2)<<endl; 2 ]9 m2 y9 m/ ~* W/ L* s1 D0 ` cout<<a<<" "<<b;7 F4 {, |3 ~, R o$ c g=a+fi(n-2-i)*(b-a)/fi(n-i); * y) `" z. \. ]* q$ S. R6 k j=a+fi(n-1-i)*(b-a)/fi(n-i);) _/ e9 `! O/ _ cout<<g<<" "<<j; 7 S& e- \: S# E2 c3 D N Z6 D for(i=1;i<n-1;i++) 5 t* Z' O9 ~5 w/ [ M% B7 Q( Q1 r. _ { 5 C0 A8 v' r( Y, M* q9 t/ _ ( s4 V4 n$ _/ b . w4 x2 v2 D* T1 g4 \( k A if(fx(n,p,g)>fx(n,p,j)). {0 c1 y) P R9 Z: h {7 [' P( i* D+ {8 b# b) r. ~ a=g;: J- E6 ]) s2 C3 a2 [, y g=j; 3 W8 @; s8 s; C9 ` j=a+fi(n-1-i)*(b-a)/fi(n-i); u5 u) d) o. E }! u& H# t7 E# x' z! x2 n7 t; [# D else , F. h) L) q% k8 Q' j3 t { 0 S1 _0 l$ V6 T! d: ?! y. Z b=j; m- R4 ?! U M8 c j=g; ; b8 F5 e! M9 x5 O g=a+fi(n-2-i)*(b-a)/fi(n-i); 7 Z+ S# C" c# a M- u; i } ( e" K" t# w% B( x 8 J7 ` L- \5 ? } ) ~2 P* i; _* a( h$ o* q j=g+c; 1 P! w! F- c0 R( _) A6 ?/ d if(fx(n,p,g)>fx(n,p,j)) & a, M: q- C4 R/ z a=g;4 x" E e7 P1 A# j% f2 t else b=g; 7 \0 G) f Z7 X1 y! u6 t* \; n2 W cout<<"函数值在 ("<<g<<","<<j<<")之间;"; + h. j$ X- m0 t8 j0 q ) `+ d- |- ~0 ?' K }. E: g# n; u4 v2 ?* m+ [( l( v

zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持2 反对反对0 微信微信

0

主题

0

听众

18

积分

升级  13.68%

该用户从未签到

新人进步奖

回复

使用道具 举报

Allos        

13

主题

2

听众

103

积分

升级  1.5%

该用户从未签到

回复

使用道具 举报

zhanfeng        

8

主题

2

听众

152

积分

升级  26%

该用户从未签到

回复

使用道具 举报

ilikenba 实名认证       

1万

主题

49

听众

2万

积分

  • TA的每日心情
    奋斗
    2024-6-23 05:14
  • 签到天数: 1043 天

    [LV.10]以坛为家III

    社区QQ达人 新人进步奖 优秀斑竹奖 发帖功臣

    群组万里江山

    群组sas讨论小组

    群组长盛证券理财有限公司

    群组C 语言讨论组

    群组Matlab讨论组

    回复

    使用道具 举报

    harleych        

    0

    主题

    0

    听众

    19

    积分

    升级  14.74%

    该用户从未签到

    新人进步奖

    回复

    使用道具 举报

    0

    主题

    0

    听众

    17

    积分

    升级  12.63%

    该用户从未签到

    新人进步奖

    回复

    使用道具 举报

    zhyi        

    4

    主题

    2

    听众

    25

    积分

    升级  21.05%

    该用户从未签到

    新人进步奖

    回复

    使用道具 举报

    mengfanqi        

    2

    主题

    2

    听众

    150

    积分

    升级  25%

    该用户从未签到

    回复

    使用道具 举报

    garr        

    0

    主题

    0

    听众

    19

    积分

    升级  14.74%

    该用户从未签到

    新人进步奖

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-4-16 01:35 , Processed in 1.504568 second(s), 101 queries .

    回顶部