查看: 5961|回复: 3

求助！！数值分析程序（c语言）

字体大小: 正常放大

西西

1 主题	0 听众	20 积分

升级 15.79%

该用户从未签到

电梯直达

1^#

发表于 2004-5-10 08:46 |只看该作者 |正序浏览

|招呼Ta 关注Ta

<

>本人紧急求助，希望各位大虾们伸出你们友情之手，帮我一下。
<

>1。牛顿迭代法
<

>1。全主消元法
<

>3。改进平方根法
<

>4。牛顿向前，向后插值法
<

>5。加权最小二乘法
<

>6。龙贝格积分法
<

>拜托了。。。。

zan

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

使用道具举报

陈了

15 主题	3 听众	553 积分

升级 84.33%

TA的每日心情

	开心 2014-11-16 19:37

签到天数: 16 天

[LV.4]偶尔看看III

自我介绍: 数学，没有不可能

4^#

发表于 2010-3-22 11:13 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

哇哇哇哇哇哇。。。好复杂呀！看不懂

我也说一句

发表

使用道具举报

西西

1 主题	0 听众	20 积分

升级 15.79%

该用户从未签到

3^#

发表于 2004-5-14 10:54 |只看该作者 |招呼Ta 关注Ta

？

<

>谢谢你哦，呵呵

我也说一句

发表

使用道具举报

loveyuan

21 主题	2 听众	104 积分

该用户从未签到

2^#

发表于 2004-5-11 08:32 |只看该作者 |招呼Ta 关注Ta

<

>1。牛顿迭代法<

>
 #include "stdio.h"
 #include "math.h"
 int dnewt(x,eps,js)
 int js;
 double *x,eps;
 { extern void dnewtf();
int k,l;
double y[2],d,p,x0,x1;
l=js; k=1; x0=*x;
dnewtf(x0,y);
d=eps+1.0;
while ((d>=eps)&&(l!=0))
 { if (fabs(y[1])+1.0==1.0)
 { printf("err\n"); return(-1);}
 x1=x0-y[0]/y[1];
 dnewtf(x1,y);
 d=fabs(x1-x0); p=fabs(y[0]);
 if (p>d) d=p;
 x0=x1; l=l-1;
 }
*x=x1;
k=js-l;
return(k);
 }<

>全主消元法<

>#include "stdlib.h"
 #include "stdio.h"
 int acgas(ar,ai,n,br,bi)
 int n;
 double ar[],ai[],br[],bi[];
 { int *js,l,k,i,j,is,u,v;
double p,q,s,d;
js=malloc(n*sizeof(int));
for (k=0;k<=n-2;k++)
 { d=0.0;
 for (i=k;i<=n-1;i++)
 for (j=k;j<=n-1;j++)
 { u=i*n+j;
 p=ar*ar+ai*ai;
 if (p>d) {d=p;js[k]=j;is=i;}
 }
 if (d+1.0==1.0)
 { free(js); printf("err**fail\n");
 return(0);
 }
 if (is!=k)
 { for (j=k;j<=n-1;j++)
 { u=k*n+j; v=is*n+j;
 p=ar; ar=ar[v]; ar[v]=p;
 p=ai; ai=ai[v]; ai[v]=p;
 }
 p=br[k]; br[k]=br[is]; br[is]=p;
 p=bi[k]; bi[k]=bi[is]; bi[is]=p;
 }
 if (js[k]!=k)
 for (i=0;i<=n-1;i++)
 { u=i*n+k; v=i*n+js[k];
 p=ar; ar=ar[v]; ar[v]=p;
 p=ai; ai=ai[v]; ai[v]=p;
 }
 v=k*n+k;
 for (j=k+1;j<=n-1;j++)
 { u=k*n+j;
 p=ar*ar[v]; q=-ai*ai[v];
 s=(ar[v]-ai[v])*(ar+ai);
 ar=(p-q)/d; ai=(s-p-q)/d;
 }
 p=br[k]*ar[v]; q=-bi[k]*ai[v];
 s=(ar[v]-ai[v])*(br[k]+bi[k]);
 br[k]=(p-q)/d; bi[k]=(s-p-q)/d;
 for (i=k+1;i<=n-1;i++)
 { u=i*n+k;
 for (j=k+1;j<=n-1;j++)
 { v=k*n+j; l=i*n+j;
 p=ar*ar[v]; q=ai*ai[v];
 s=(ar+ai)*(ar[v]+ai[v]);
 ar[l]=ar[l]-p+q;
 ai[l]=ai[l]-s+p+q;
 }
 p=ar*br[k]; q=ai*bi[k];
 s=(ar+ai)*(br[k]+bi[k]);
 br=br-p+q; bi=bi-s+p+q;4 ^9 W* m \) c9 ]0 t) `( |; O: x8 X
 }( j6 C* P: N- V, H5 S Q- C' l
 } w2 c5 I% u' M+ _5 g5 a) q, L
u=(n-1)*n+n-1;; ?' v+ {3 d8 v
d=ar*ar+ai*ai;
' a0 R$ R; b3 o1 a/ S if (d+1.0==1.0)
; y! y; @$ _+ ~6 {( |2 t1 j: u) ?1 m { free(js); printf("err**fail\n");5 f6 F6 ^6 Q5 [: `
 return(0);0 S1 `0 j2 @7 w4 O& d9 J: ^
 }. o2 X7 S( B F- O% V. z" |% p
p=ar*br[n-1]; q=-ai*bi[n-1];9 A0 W1 ^2 u: x3 T, |9 F5 T; T8 |# i
s=(ar-ai)*(br[n-1]+bi[n-1]);
/ \ |$ V5 c- G' R* z br[n-1]=(p-q)/d; bi[n-1]=(s-p-q)/d;1 A6 M4 x0 y2 g" p7 f8 A0 `! v# M" z
for (i=n-2;i>=0;i--). p4 J# y. I; b$ r- J3 a9 |
for (j=i+1;j<=n-1;j++)' [" a! m9 W' K" V+ M
 { u=i*n+j;
* X; A/ Y! n I3 U- A+ g6 r p=ar*br[j]; q=ai*bi[j];
. f7 `$ U7 d( T( ^0 ], d s=(ar+ai)*(br[j]+bi[j]);
+ t+ J( ]- j3 W, ]" v/ V$ G6 x, ]0 E br=br-p+q;9 g i8 Y, ?! T
 bi=bi-s+p+q;
/ u3 ~' G2 b# B$ e* A1 f }, R4 K% N" s2 o. _
js[n-1]=n-1;
2 |) {/ Q! K% C1 y for (k=n-1;k>=0;k--)$ p1 E4 N6 R( {0 {& i( A( I
 if (js[k]!=k)% T6 T2 f. t# A. g* H
 { p=br[k]; br[k]=br[js[k]]; br[js[k]]=p;
. ~( g/ c. q7 F u; W* T. L p=bi[k]; bi[k]=bi[js[k]]; bi[js[k]]=p;8 X% c+ l; e6 w1 a5 L
 }
7 W2 M) S- W6 N2 m+ [6 Q free(js);8 R7 h# E! n4 F& P9 n/ U+ d) ^
return(1);' q8 R5 N& _6 `8 B+ M
 }<>平方根法<>#include "math.h"
 r. Y& v+ m" \' i #include "stdio.h"
2 m8 @: u9 z# a5 f% @' P int achol(a,n,m,d)
$ ]$ Q/ L% O( ]3 u4 m3 W int n,m;0 Y7 X7 V' a$ O4 g0 C8 F
 double a[],d[];
4 b2 D- k( N7 m3 m6 n7 Q; c { int i,j,k,u,v;
4 n0 `3 N# O2 @3 B" v( M5 a9 U if ((a[0]+1.0==1.0)||(a[0]<0.0)), }# Y" s2 @1 w, `9 D5 S
 { printf("fail\n"); return(-2);}8 \- W8 |, _& ?! C) ]0 c. Z6 D
a[0]=sqrt(a[0]);* Q0 w' Y, P4 L. |+ J# E& P3 [
for (j=1; j<=n-1; j++) a[j]=a[j]/a[0];
1 ~7 k9 R$ M: p5 h2 E" F- _/ q! d9 G for (i=1; i<=n-1; i++)8 S" U: s5 s; `: m9 u9 ]
 { u=i*n+i;; _! X2 r( h7 P8 H
 for (j=1; j<=i; j++)
. s% l# b- M' J9 l1 v+ i$ K! o; l { v=(j-1)*n+i;, ~1 V9 K, L R
 a=a-a[v]*a[v];: F1 s! Q( |- @- X! E
 }
; v& `7 a( J- X6 Q$ q if ((a+1.0==1.0)||(a<0.0))
+ W8 O1 Z0 J. k6 H- ~( y { printf("fail\n"); return(-2);}
1 P2 K! n( j% G, z( v5 F) \ a=sqrt(a);" C; @8 g$ A4 R- e$ @
 if (i!=(n-1)) ?' M/ G! p7 `) L& W
 { for (j=i+1; j<=n-1; j++)6 A% ~$ ~+ s4 E! E! ~- _+ h
 { v=i*n+j;3 w" g8 r0 W, H6 N
 for (k=1; k<=i; k++)
& V1 }7 Q ~, @& R$ T a[v]=a[v]-a[(k-1)*n+i]*a[(k-1)*n+j];: |: C$ z2 F2 E8 U5 e) H+ j; U
 a[v]=a[v]/a;5 y8 H8 W& ?2 x5 d
 }
0 c) d- `$ r/ f }
1 q* v: ?" ]- q# X }7 x2 [5 R" @2 F- B3 j3 q" t2 s
for (j=0; j<=m-1; j++), U5 W- @5 y& d( V, W5 [; j
 { d[j]=d[j]/a[0];
9 T0 b8 N4 X7 s# q" U- U% Q for (i=1; i<=n-1; i++)0 p' l3 {0 W6 d- [, R! n: v
 { u=i*n+i; v=i*m+j;
& k! Z+ J/ {( ]4 Y for (k=1; k<=i; k++)
0 V* a* S6 \. C* G$ M* E d[v]=d[v]-a[(k-1)*n+i]*d[(k-1)*m+j];' C E& V; I9 ^* D- I, E/ e2 E, u
 d[v]=d[v]/a;
' `& U( f4 V# a- B* S0 ^ }, Z$ n: U: B0 z2 s9 o: \
 }* K% _/ c& m2 z0 U" R
for (j=0; j<=m-1; j++), o* {# R) U! I) [
 { u=(n-1)*m+j;
* v) V; y( M0 E/ ~% S+ O; S6 @ d=d/a[n*n-1];3 V# d! a: f; h g
 for (k=n-1; k>=1; k--)
- x$ F6 ^3 d3 r { u=(k-1)*m+j;
6 ]+ P* E5 y6 V! T, g5 P9 C- u' ]8 l5 ] for (i=k; i<=n-1; i++)% u" y5 p3 |+ F$ B
 { v=(k-1)*n+i;
, Q8 r* Y5 o* a d=d-a[v]*d[i*m+j];
$ _4 D7 S$ [/ j" H& F { }! v' W2 j. v# B7 \, B( M
 v=(k-1)*n+k-1;" ^+ O$ w. K4 E c) p- f+ z
 d=d/a[v];4 y1 ?: J3 L9 }: ?
 }4 i6 b1 h! `8 G! N6 _
 }7 m) e* o+ J+ `
return(2);
6 v0 z* `- ^/ W; B; } }<>牛顿向前，向后插值法<>double eelgr(x0,h,n,y,t)% r: b6 b3 c ~! N
 int n;8 [' n- L& R, u3 a/ v: }& A
 double x0,h,t,y[];
 k) o- N4 J( `& d; K1 u9 v { int i,j,k,m;
! K0 K8 b: d U' W double z,s,xi,xj;
4 G* U/ D. N: g; h4 {0 Q, b float p,q;, N% t7 h* t0 F5 Y4 u2 l- `+ K
z=0.0;
* y& ?( n% n; f }+ o if (n<1) return(z);
# ]$ F+ N+ W7 b1 S, n9 w if (n==1) { z=y[0]; return(z);}
- Q6 H1 z- w8 \- o/ X if (n==2)
! E/ s Y1 d' m) [+ I: ~) Y { z=(y[1]*(t-x0)-y[0]*(t-x0-h))/h;! e9 E, K' n# s6 M
 return(z);
* q' |, |0 A# h! D ]/ h }
# l; t& X4 L' R2 L- C if (t>x0)1 Y3 [ u9 _4 I3 e( s* U5 k0 A( P
 { p=(t-x0)/h; i=(int)p; q=(float)i;0 q/ E3 x r- W3 L% ?! N
 if (p>q) i=i+1;- K( \% @0 F& k- g
 }( Z( K) N( ^! |! j7 b+ {
else i=0;
, B7 z- Y9 I9 f/ u. q k=i-4;
' Y) ?4 J& R. L if (k<0) k=0;
# m! L0 v6 l2 f5 H% Q/ [! c m=i+3;( w6 |7 L/ w8 x- a9 Q
if (m>n-1) m=n-1;( z; S) ?# ~2 \2 ]+ [
for (i=k;i<=m;i++)
* i7 C% Q- _/ z: H; e( S+ S { s=1.0; xi=x0+i*h;9 ^ s2 R, u L f7 i
 for (j=k; j<=m; j++)+ C6 t9 {; G5 u- U3 l
 if (j!=i)* Y- S' y5 u( U7 C8 O# u
 { xj=x0+j*h;
3 a3 R8 d, \) m' q5 w5 n% h: e s. b s=s*(t-xj)/(xi-xj);5 R# r6 N- D0 n0 f
 }" Z e# ]8 E A# y5 K2 ?5 U& |6 y
 z=z+s*y;
5 [9 Z. Q5 d) w9 o }
. X: {- P8 u2 E- k. w6 d/ ? return(z);# Q0 J: \: @- W4 \7 ]( ?, r
 }; b: p+ m2 |0 o, U
向前，向后是一样的思想！！<>加权最小二乘法<>#include "math.h"6 k$ ?, y1 a' ?- y
 void hpir1(x,y,n,a,m,dt)
( c0 P7 P: L' B: j A int n,m;
: `6 z" x$ ~! t double x[],y[],a[],dt[];0 T" V8 E) E9 c3 l( c* _
 { int i,j,k;
2 F1 _$ a$ a. `+ [ double z,p,c,g,q,d1,d2,s[20],t[20],b[20];
6 ^7 o6 b, {# y* k5 {6 V' O: A for (i=0; i<=m-1; i++) a=0.0;
' h: R6 t6 d: q% U1 g if (m>n) m=n;8 y9 E- Q! q/ O
if (m>20) m=20;7 R4 d, o {9 A! W8 i
z=0.0;
' ?( k4 A3 s) F& Z; b6 L for (i=0; i<=n-1; i++) z=z+x/(1.0*n);# b- _# B( U7 ~! z+ F5 t% I5 G
b[0]=1.0; d1=1.0*n; p=0.0; c=0.0;1 @) c1 M& A) @- J) l
for (i=0; i<=n-1; i++)
. b* b4 U; y# x { p=p+(x-z); c=c+y;}
! e1 `. p% q( B" z2 f. q) m8 @5 s c=c/d1; p=p/d1;
- m$ ?6 i/ r( S# y1 M& M0 F) m a[0]=c*b[0];5 _" O5 h4 B/ ]6 U8 o' Q) ~' q
if (m>1)
: E! G# V; y. ]8 c7 i { t[1]=1.0; t[0]=-p;
. I% X. f1 G# o d2=0.0; c=0.0; g=0.0;
3 B4 \* V3 s: g! s! }. ]3 H for (i=0; i<=n-1; i++)2 e, K- l8 A0 g# |4 J- S3 P z( v
 { q=x-z-p; d2=d2+q*q;
' Q1 q% E6 h8 @) t o" E, o" S c=c+y*q;
* o* s& c0 H+ {0 H0 j g=g+(x-z)*q*q;
. X+ K7 w& W F% R: H }
7 [' i1 g) I- e8 b" e, R c=c/d2; p=g/d2; q=d2/d1;! Z5 `* V$ a3 e4 D( N
 d1=d2;3 Q% j1 x+ Q( `' }( x, H( [7 B
 a[1]=c*t[1]; a[0]=c*t[0]+a[0];
$ o9 M. G1 ], `0 }) D! T }
5 n% l0 @; p0 _ for (j=2; j<=m-1; j++)
& k% m& i: d% @$ _4 f0 |; ~' R { s[j]=t[j-1];, V, E; a, h2 t" Y
 s[j-1]=-p*t[j-1]+t[j-2];
0 g& r1 p. B; F# e2 q if (j>=3)4 U9 Z% F$ j1 Q/ f
 for (k=j-2; k>=1; k--)
/ J: j# {0 f) z& P' u s[k]=-p*t[k]+t[k-1]-q*b[k];
1 b `1 c$ ?' d8 I6 u- K s[0]=-p*t[0]-q*b[0];
' l3 ]/ r! N5 i+ F# E' I d2=0.0; c=0.0; g=0.0;% y: W8 D& T, |; T0 H6 B: H
 for (i=0; i<=n-1; i++)
* N* d+ c2 r4 z$ i1 n) a* z { q=s[j];/ C* \0 S, n: A6 C$ ~( a
 for (k=j-1; k>=0; k--)- i! i% `' K# e* S; i6 ~+ A. l
 q=q*(x-z)+s[k];
6 Z% P3 F) a! s5 v2 f4 D d2=d2+q*q; c=c+y*q;0 ?* {. `# i: V" \9 [
 g=g+(x-z)*q*q;" z0 p( z8 D9 {3 N0 v
 }
$ m: c# x4 `( L9 j0 ~ c=c/d2; p=g/d2; q=d2/d1;
( o1 r' u" z' i8 N6 G d1=d2;
- ?% j4 _+ K* _4 \8 ^ a[j]=c*s[j]; t[j]=s[j];
- v! \0 u5 r0 ~1 ?& n2 z for (k=j-1; k>=0; k--)# }5 M8 i' @% o: A7 S: c
 { a[k]=c*s[k]+a[k];6 A6 Q& d, I- j5 P
 b[k]=t[k]; t[k]=s[k];7 u+ O. K) J2 l
 }
6 ]4 y/ x( Q7 R6 Q+ u/ N, B. |2 f }
0 A# J2 n$ o' C. C" y1 x dt[0]=0.0; dt[1]=0.0; dt[2]=0.0;
0 _7 n# B6 `. s for (i=0; i<=n-1; i++)
% |$ o' ^- V0 ?) l I2 ?0 I { q=a[m-1];
1 V- E0 I# m1 x. ^4 t- e1 N6 \) m; Y for (k=m-2; k>=0; k--)
. @* K2 x8 @* S q=a[k]+q*(x-z);/ C' ^( W% u# y7 t9 U; Q% s
 p=q-y;
% K t: I' r% v if (fabs(p)>dt[2]) dt[2]=fabs(p);& _0 R# n- j! b' V
 dt[0]=dt[0]+p*p;3 Y" z! |, _) m$ b9 {# X
 dt[1]=dt[1]+fabs(p);
/ Y, p* }7 A1 {- `8 b: n3 f3 O* P }' _; d* C8 @, E h3 Z6 O8 N
return;
7 w( Y6 o/ a% ~" ?' c }<>龙贝格积分法<>#include "math.h"+ X; i% |+ S) ?# G# F5 p, M' g5 d
 double fromb(a,b,eps)
" X4 `" y% }3 D- ?- m6 {; z double a,b,eps;5 v0 Q- b% ] b B9 P
 { extern double frombf();2 J- x' u7 k/ ~9 ~0 A
int m,n,i,k;: D7 u' L# w6 h# W |" o4 K( c9 ^
double y[10],h,ep,p,x,s,q;
; A; S+ e6 S. ^$ g h=b-a;
, H/ m0 [' x- u y[0]=h*(frombf(a)+frombf(b))/2.0;9 k4 a) H0 d* M3 Z8 ?/ M! E
m=1; n=1; ep=eps+1.0;2 Y4 y6 ~' t2 f! Z# Q# W% |
while ((ep>=eps)&&(m<=9)). \8 j! z+ E' `5 z
 { p=0.0;; i4 l x v; d
 for (i=0;i<=n-1;i++)
8 y# f! k7 z- q' ]* w { x=a+(i+0.5)*h;
; s2 ^1 P% a: R l& \; E. M p=p+frombf(x);
1 h0 J: m" W* A8 G% x! w }
5 u( R2 {* W( S( z0 ^ p=(y[0]+h*p)/2.0;# w Z3 ]6 ~& w2 z9 X, C
 s=1.0;
# E1 e- o! M6 @3 l) s/ {$ N for (k=1;k<=m;k++)2 M& e% |1 W6 g0 W# E/ q
 { s=4.0*s;: g# z; h# w# T `( F' H/ F% J$ d
 q=(s*p-y[k-1])/(s-1.0);5 G' \& L L% V5 w
 y[k-1]=p; p=q;4 j, y+ X+ I6 N( G# ^; e
 }
& z/ I3 i" }* i6 A j7 E ep=fabs(q-y[m-1]);
; [! k+ x8 T6 o+ v, t& G5 O& h- o: D m=m+1; y[m-1]=q; n=n+n; h=h/2.0;2 a% D) L# N$ R) W, F
 }
8 O" E/ W& m4 j0 L7 v( Z; j3 ]2 \ return(q);
$ b. N9 a! v+ D }<>呵呵希望对你有用！！

我也说一句

发表

May God grant me the wisdom to discover the right and strength to make it endure 因为陌生，所以勇敢，因为距离，所以美丽。

使用道具举报

返回列表

帐号		密码		只需要一步，快速开始		注册地址	找回密码

求助！！数值分析程序（c语言）

？

浏览过的版块

新人进步奖

元老勋章

QQ

电话咨询

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

求助！！ 数值分析程序（c语言）

？

浏览过的版块

新人进步奖

元老勋章

QQ

电话咨询

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

求助！！数值分析程序（c语言）