QQ登录

只需要一步,快速开始

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

非线性规划程序

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

1

主题

0

听众

19

积分

升级  14.74%

该用户从未签到

新人进步奖

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

#include<iostream.h> 4 z4 Z$ k1 u( a3 \ P#include<math.h>+ n& y1 y. _: F( w double fx(int n,double p[],double x)+ k8 ?3 K& f4 o4 u { 8 }& t' [5 o" i& F/ ~. j' d double z=0,y=1; * ?1 x. e5 g, S$ P int i=0,j=0; - V# Y; w! h4 F( U3 M, ] for(;i<n+1;i++)' m* ^4 a( j1 ]- o$ ]* x { / @* l; ~8 E5 M" l9 A for(;j<i;j++) 6 a+ H7 t/ _$ x0 K# m+ `7 J$ ] y*=x; * \3 Y. ~, e: v7 O% O z+=*(p+i)*y;1 Q; V" T! q9 T5 A8 ^ } 9 @# N, ?4 Y/ P, Q: Z return z; * U* m. _2 O2 @- X; p4 Y5 k: W% g}

' p) N- S) s/ D& c3 i! K+ b

int fi(int k)3 r0 c1 J) A4 g t4 A/ C7 P" O+ B { 7 z% E1 w! Y2 r! G# V1 v if(k==0)8 ]+ w4 X+ T4 p' O2 z& {" k" J return 1;: f$ ^: ^7 E6 r, k* d else ) t# k* @% d& N { ! e! l) Y6 Q& J/ A if(k==1) ' n( P+ U0 |& X# I$ I' J return 1; + w4 c/ r6 d5 Q' H5 |7 I7 X else $ @3 ~ t Z. l4 w; W {! n) ]' k4 l& ]; w6 l int i=0,i1=1,i2=1,w=0;; U6 ~0 j* S8 x* P2 w6 E4 m( W for(;i<k-1;i++)/ B& l4 |4 A) ^, K {6 V' a2 Z# ]8 n' P. A w=i1+i2;0 S: W5 |' `0 A4 N; d i1=i2;1 R$ a; F3 H6 m% d" W i2=w;5 E: I) {3 u7 ?4 H- g2 ] } 3 R# ~/ @* f; N2 k( V8 S return w;+ S% R' ?, N+ |0 X" N3 @ }# v! i! b1 K. N c# [ J5 h } Y: G4 _' `4 M. ]$ `/ K}

5 `: C! K. Q( x4 A, `% N

int r(double w,double a,double b) - X" q; u" l v{: [9 r. @8 g8 { int i=0;5 y, r: I ? V1 E for(;fi(i)<=(b-a)/w;i++) : E/ b, N4 V* |& Y4 F {}( Y2 |6 q' F' k/ O" c; M* \( U8 e! r return i;- ?8 [ O% T, l2 Q0 n I. B4 ` }. J" X" D: x; Q void main() 0 J, A1 s) ?" U/ I" I) |{( `" R8 i8 v& m/ \ cout<<"输入函数的最高次幂";! x E( D. L/ L! P9 J int n; 9 j1 }, z' U) H6 M8 w1 T cin>>n;8 c; ?& e% L/ o) i+ S double *p=new double[n+1]; 1 l1 s/ ?7 w X% F C% G, F double a,b,g=0,j=0,c=0;: Q. u; q9 U2 z! C/ G cout<<"依升降幂输入x项的系数"<<endl;: y# T6 _+ Q1 P- S for(int i=0;i<n+1;i++) # t' U7 C5 D' N0 K1 `7 O* c cin>>*(p+i); " M' S% N# }" _& g; E cout<<"输入误差值"<<endl; 7 Y4 X( c. Y6 r% ]2 }- h: B double w; 1 s, l4 O6 X- `( A$ i; G cin>>w; 2 j$ @( n( F0 A4 d$ K9 a cout<<"输入初始区间"<<endl; - c7 w; V& y. b! G cin>>a>>b;& j+ {; l+ m. Q5 a7 T' {' n cout<<"输入辨别常数"<<endl; ; C- \" I) V8 j3 A# P; H6 ~$ v cin>>c;6 R, n) H. d1 M2 ~6 l n=r(w,a,b);; P4 [5 }: t# {+ E* z i=0; p1 l6 R' K' g6 s; }% H9 E cout<<fi(n-i-2)<<endl; & ]8 q) ~+ R" X6 x1 h& V' x cout<<a<<" "<<b;5 N' G" r# \7 Y2 n A g=a+fi(n-2-i)*(b-a)/fi(n-i);" [9 A5 ^% M6 ^- M4 _# f j=a+fi(n-1-i)*(b-a)/fi(n-i); " ^$ X9 U" e8 \ cout<<g<<" "<<j;) H! G' D( f; H: x$ p- d) x- o; c for(i=1;i<n-1;i++) + |# `5 X2 B8 R; F- g* P& } { 6 `6 M% {& t9 D7 Z8 ` & Q8 Y7 X5 ?3 T, U7 l; a / k4 o# m, `* V" I R3 Q if(fx(n,p,g)>fx(n,p,j)) , Y# n6 ^' {# k8 {7 k+ c {4 Y9 c- D8 i0 b Y a=g;9 k! P) |, f2 k2 ?- v% C g=j;* j5 |- X: M% m j=a+fi(n-1-i)*(b-a)/fi(n-i);4 u' L& A( B& D3 Q# M }0 x- _ _: c# o$ j# Q else . i6 }' G; y$ s1 x! S2 h+ E {# Y* P" `' T7 g% ^: M b=j;7 X3 |1 l5 ?" r4 y( t j=g; ! g0 L+ m1 r+ A* C! D9 ` [5 k g=a+fi(n-2-i)*(b-a)/fi(n-i);) Z" s, C. E, `7 R) Y } ! w: l$ U* B7 g ' [# Q5 v5 h9 }4 h } 9 S, j2 h% x. a y2 D2 x1 r j=g+c; # s+ m# s" C3 ]: X( q% O if(fx(n,p,g)>fx(n,p,j))2 T! h8 o t5 Z H a=g; 0 t. x* j( I% B% J1 ?/ Y) S else b=g;+ D/ n% A7 x- ]3 q# E0 K0 I( w cout<<"函数值在 ("<<g<<","<<j<<")之间;"; - P, b$ {, S0 x% y: F & }) D; M( B0 \1 W! _9 H ` W} & j/ v J9 S) U+ w' |! `( K

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

0

主题

12

听众

112

积分

升级  6%

  • TA的每日心情

    2015-5-29 09:03
  • 签到天数: 8 天

    [LV.3]偶尔看看II

    邮箱绑定达人 社区QQ达人

    群组2015美赛备战交流群组

    回复

    使用道具 举报

    9

    主题

    5

    听众

    399

    积分

    升级  33%

  • TA的每日心情
    开心
    2017-3-19 15:39
  • 签到天数: 112 天

    [LV.6]常住居民II

    社区QQ达人

    群组学术交流A

    群组MCM优秀论文解析专题

    群组2013认证赛A题讨论群组

    回复

    使用道具 举报

    lsyzzhl        

    0

    主题

    4

    听众

    42

    积分

    升级  38.95%

    该用户从未签到

    回复

    使用道具 举报

    0

    主题

    4

    听众

    9

    积分

    升级  4.21%

    该用户从未签到

    沙发,板凳都没偶的份了。
    0 L- Y& ?) @2 D- B9 k) c; W/ L: a% g4 u5 i7 J: {& z+ k, P
    % X2 \2 Q; x; {  @# b9 \  Q
    " K: V3 Y, I& v' D$ J% N
    * I7 ~  R' Z) x. O& _5 S$ K$ i
    ) D) r( k6 `5 d  V9 ^3 p. N

    1 k; o; {! b4 A7 o& Y& c3 ~8 N9 U; ]- j$ A% U

    + @. n% X( |4 G  z7 D. p
    $ W) x; F& L  V  H" a5 n6 t% e& Q# t
    : q" v* K  _, Y& |2 ^
    . Y" e/ x' }( m4 s9 b" L2 A9 P: {4 ]( X) d) m: ^

    * F. \/ q) r. V0 A/ B. Q8 a& t 淘宝专享地址:zhfqw.taobao.com/
    0 d/ |8 I2 H. F+ p5 q% ^1 R" }9 ]
    2 e' [# e; T9 a6 l+ B如何将宝贝卖向全世界           索马里论坛**598专卖
    # ^' S. ^8 ^6 S3 V, w* `5 Z# T! {
    ! ^8 q$ [" g3 ^/ f4 r7 t) x
    如何将宝贝卖向全世界           索马里论坛**598专卖
    回复

    使用道具 举报

    0

    主题

    4

    听众

    744

    积分

    升级  36%

  • TA的每日心情
    无聊
    2012-2-23 09:12
  • 签到天数: 108 天

    [LV.6]常住居民II

    回复

    使用道具 举报

    21

    主题

    7

    听众

    3435

    积分

    升级  47.83%

  • TA的每日心情

    2014-5-25 20:58
  • 签到天数: 20 天

    [LV.4]偶尔看看III

    新人进步奖 优秀斑竹奖

    群组Matlab讨论组

    群组小草的客厅

    群组数学趣味、游戏、IQ等

    群组C 语言讨论组

    群组我行我数

    回复

    使用道具 举报

    515genius        

    0

    主题

    3

    听众

    12

    积分

    升级  7.37%

    该用户从未签到

    回复

    使用道具 举报

    0

    主题

    3

    听众

    19

    积分

    升级  14.74%

    该用户从未签到

    回复

    使用道具 举报

    999mmg        

    0

    主题

    3

    听众

    22

    积分

    升级  17.89%

    该用户从未签到

    新人进步奖

    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-7-27 19:36 , Processed in 2.106091 second(s), 102 queries .

    回顶部