查看: 5964|回复: 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;
5 R% n& F3 m+ S X }! H7 ?6 s( l% C% }
 }& }* Z) V3 x4 _& f
u=(n-1)*n+n-1;6 q- }! U- {. Z/ X5 T2 ]/ f
d=ar*ar+ai*ai;
4 x9 N$ v7 s" j- B ^ if (d+1.0==1.0)) A% r8 N# E0 y8 }# \% z$ C5 c
 { free(js); printf("err**fail\n");
( Z, C& Y4 |8 n2 `7 c {% U0 u6 W return(0);
 B& I) B; a' C7 [: } }
7 g) u- `! C C$ Q p=ar*br[n-1]; q=-ai*bi[n-1];& ?( [ A* H, `9 f9 O9 Y
s=(ar-ai)*(br[n-1]+bi[n-1]);/ ]" e; C1 o% ?+ C
br[n-1]=(p-q)/d; bi[n-1]=(s-p-q)/d;& ^" N8 s: B1 _
for (i=n-2;i>=0;i--)9 z Q% u# m3 V/ Q
for (j=i+1;j<=n-1;j++)$ P8 Q1 R; |4 ?9 F' q9 L
 { u=i*n+j;, r1 |5 B6 W, L
 p=ar*br[j]; q=ai*bi[j];
) J7 K$ _6 v$ y1 } s=(ar+ai)*(br[j]+bi[j]);
/ p R( H2 Z! M/ | br=br-p+q;8 [3 A! V" g8 z
 bi=bi-s+p+q;
# x( o" p5 P. {' F# R6 E }
) U6 j3 @: U; K5 Z V1 I js[n-1]=n-1;
4 `. t, S9 e) L- t for (k=n-1;k>=0;k--)# Y1 s# L7 d. C2 o: g3 S5 n/ w
 if (js[k]!=k)" R) Q- S* B' J h
 { p=br[k]; br[k]=br[js[k]]; br[js[k]]=p;! [/ f7 Q+ D0 X3 k4 e
 p=bi[k]; bi[k]=bi[js[k]]; bi[js[k]]=p;
 N* F- |9 n. b5 P7 z3 }& p; z4 P }
, @# Z1 _; n( I! N, J, J free(js);
5 X0 |. U! @. y1 r' P9 } return(1);+ J" X# i% ?/ W- S& L& f
 }<>平方根法<>#include "math.h"
6 l: T% |8 P- ^/ T- t1 @/ U- G #include "stdio.h"
; G( h4 T3 I/ e int achol(a,n,m,d)
4 x) q4 o3 W0 V: z4 S int n,m;
6 S- t! G% f# j2 E9 B1 O# \ double a[],d[];
" h2 w6 q+ i* ~) l, |' A8 V { int i,j,k,u,v;) }! k! l% r; l. n. [) L+ i5 X
if ((a[0]+1.0==1.0)||(a[0]<0.0))
! ]5 {! C+ c" Q' j6 L6 W( O { printf("fail\n"); return(-2);}& @/ K4 U& U& D) b9 U
a[0]=sqrt(a[0]);8 X9 e" g, h" g) Q
for (j=1; j<=n-1; j++) a[j]=a[j]/a[0];+ h* D* p; ~) A/ w# H/ F+ k: T2 h
for (i=1; i<=n-1; i++)
+ |/ `+ O$ j# \, d9 N; o6 D; g { u=i*n+i;9 \' r8 J) R2 M0 i( ]
 for (j=1; j<=i; j++)
% Y7 n! r Q" W: G& ~( P3 m" L { v=(j-1)*n+i;7 r; R: E7 c4 _+ j
 a=a-a[v]*a[v];
$ X. p* }7 O0 b }
; M7 U8 m( |# ~' L if ((a+1.0==1.0)||(a<0.0))
7 j9 v' l: Q, T2 m* L { printf("fail\n"); return(-2);}
+ K1 }3 ?7 w! z2 A+ o4 I1 \4 r a=sqrt(a);
: K3 z5 e! E& ]( H1 J, c1 p if (i!=(n-1))
/ r+ M. p* O3 O { for (j=i+1; j<=n-1; j++)
7 q8 U0 T2 ~# M& x, Q { v=i*n+j;% z7 j: u- m0 y. j, t* B
 for (k=1; k<=i; k++)$ ]$ W% O- Q7 c7 M D
 a[v]=a[v]-a[(k-1)*n+i]*a[(k-1)*n+j];( ^$ o4 M. N1 S* Q8 G0 t
 a[v]=a[v]/a;
; r+ E/ e9 \# Q. j! K5 T: E8 j }' G0 o6 w$ }! K0 S! f
 }
" `+ @4 t, v1 n. A% w: O. g* | }8 ^7 ^& F7 L2 c! s4 V' k" d; A
for (j=0; j<=m-1; j++)
) L! ^) y! b5 l { d[j]=d[j]/a[0];4 l6 m- n/ Z; a% d. w: e
 for (i=1; i<=n-1; i++)
$ ? e, M# B0 D5 d { u=i*n+i; v=i*m+j;8 K. A! [8 q" b2 _& n/ h% _4 W
 for (k=1; k<=i; k++)& _6 ]/ z+ n! f9 V8 H
 d[v]=d[v]-a[(k-1)*n+i]*d[(k-1)*m+j];! j" F# k/ T+ i0 u- |+ {/ X. m* z
 d[v]=d[v]/a;1 X* ^4 |) j. n. E" F j2 X% b
 } A+ A0 _9 {2 b0 S- z) z- h' z
 }2 Z5 J" B, b# ~ [
for (j=0; j<=m-1; j++)" N( M0 O. H- K% Y: @; b
 { u=(n-1)*m+j;
, G8 p y4 e; f* @/ M# ?2 q, g d=d/a[n*n-1];
- }5 k6 a( x$ c* `) d: Q: Z9 Y for (k=n-1; k>=1; k--)
 S& q- G; O/ R* ?0 l0 U { u=(k-1)*m+j;' I7 l3 ?' i5 i% d
 for (i=k; i<=n-1; i++): y9 Z0 r2 ]3 E
 { v=(k-1)*n+i;
$ C" k: h/ P$ W5 ?+ [( y% V d=d-a[v]*d[i*m+j];
/ I, f3 M) l/ ? }
3 p6 ]% ?% _; V3 S" S& r- |) d v=(k-1)*n+k-1;0 D+ B6 x7 X, u% Z3 o! u8 D2 h
 d=d/a[v];. l# n/ j. I! A0 W$ d! B
 }5 b( J; t+ W4 k B8 d2 U7 r+ n: R
 }
3 W0 w' C5 X2 k, I: s* `, ? return(2);
( ]! m, g2 Z) E- y }<>牛顿向前，向后插值法<>double eelgr(x0,h,n,y,t)
4 s+ E6 P* i! } int n;
2 q8 {) m5 e! z% E3 v, F. E double x0,h,t,y[];* K9 M2 W# w* P1 c! M
 { int i,j,k,m;
8 W* ^! `, z: R: n& E& k' w" m double z,s,xi,xj;4 O: Y& F' W+ W4 h2 W" O5 W* W, q
float p,q;
/ _2 v1 V4 ]& X8 `) }) Q7 ` z=0.0;. g) {6 Y6 Q7 m6 s) Y3 E6 T) V; y
if (n<1) return(z);& E) j# O. _% A2 W! q I
if (n==1) { z=y[0]; return(z);}
1 R+ z3 E* G. J* c7 \ if (n==2)# m% X& e4 v/ @. m% G* B
 { z=(y[1]*(t-x0)-y[0]*(t-x0-h))/h;
 W: o% E- j) @, |- j- S8 a& L return(z);
, r/ r* m: Y% P" \: l }
+ G% }: q& I' J* f# s: R if (t>x0)
* x/ }4 H8 b& h: K: r' M4 y { p=(t-x0)/h; i=(int)p; q=(float)i;
* |1 I4 _, J3 k% i1 s if (p>q) i=i+1;# f1 h( s# C. w) [* ^9 v2 l5 @
 }4 h+ k8 T( \) i9 {% f
else i=0;
/ i& n; E- S/ t/ i& R; ? k=i-4;
1 G1 i. a" F6 D' `" ?) k9 w if (k<0) k=0;
$ M2 U& l7 B5 U6 g m=i+3;
. V8 d) h0 p a$ U2 A4 _4 ^ if (m>n-1) m=n-1;+ d$ {4 h* s( i v; b" J5 a( b+ \
for (i=k;i<=m;i++)
+ w- n, X0 ~, p [/ _ F { s=1.0; xi=x0+i*h;
( t7 Z# d* y% T4 z for (j=k; j<=m; j++); m8 ]3 h, y7 `" R
 if (j!=i)8 w" y/ u4 \& [
 { xj=x0+j*h;
* `: [# d: o2 W- o' P s=s*(t-xj)/(xi-xj);
2 w. M" P+ r, y! S$ G$ v }/ q- Z0 p- d+ }
 z=z+s*y;% X1 q3 D+ e( ?+ N( t
 }
( ^" H1 T$ i9 O# f4 i u return(z);0 Q' y+ z- Y* g- U3 F: a
 }
+ Z; K" X U( @4 z; S向前，向后是一样的思想！！<>加权最小二乘法<>#include "math.h"# E) v6 X' x+ _' M: h
 void hpir1(x,y,n,a,m,dt); G1 @' [# s0 ?1 f! p
 int n,m;
9 R( ^6 |* Y3 f/ u6 U double x[],y[],a[],dt[];0 \7 _& R# w; L4 x, x' R+ D/ ^
 { int i,j,k;
5 [( h4 ^. s, [+ p) Y double z,p,c,g,q,d1,d2,s[20],t[20],b[20];
) }: h, Q F) j( X for (i=0; i<=m-1; i++) a=0.0;* X# z7 B4 M( z
if (m>n) m=n;" {0 r7 {. o9 Z2 B# H8 G
if (m>20) m=20;
, h0 v% H! E0 J' w! P0 S. r: k z=0.0; l O' j% l6 n9 P; T1 T5 {
for (i=0; i<=n-1; i++) z=z+x/(1.0*n);
4 U4 @" K* p, {1 ] b[0]=1.0; d1=1.0*n; p=0.0; c=0.0;
* ~/ u7 x$ z7 F- M$ Z+ w0 a for (i=0; i<=n-1; i++)
# C# a3 \: Z5 M { p=p+(x-z); c=c+y;}9 J% x X7 L) F5 ?/ m
c=c/d1; p=p/d1;1 d2 o9 g: O4 }0 V! b/ f. `: y
a[0]=c*b[0];
4 |& R8 l K7 y/ H# ] if (m>1)" I* |) b3 E3 @' ?2 Y* G
 { t[1]=1.0; t[0]=-p;9 a/ r: F0 Q+ E' F. A
 d2=0.0; c=0.0; g=0.0;. t7 e6 Y6 u' R( f9 X$ z
 for (i=0; i<=n-1; i++)8 f4 h# \; X8 i" ~
 { q=x-z-p; d2=d2+q*q;
: Q. |, S0 a& |- _# {# \& v+ Z: ? c=c+y*q;
: L% E) ^& r/ W2 R, ] g=g+(x-z)*q*q;
3 \7 n0 L/ {1 x }& R i1 t! S) G3 v# ]& d
 c=c/d2; p=g/d2; q=d2/d1;
- E# _. E' f- P) M0 b2 _" R p. K d1=d2;
3 \5 D) P) s) d; z6 p a[1]=c*t[1]; a[0]=c*t[0]+a[0];
, w" B! |' h, }& v4 P3 Z9 ` }
& L% [- I6 N' k# g for (j=2; j<=m-1; j++)+ O! F$ y b: j# i7 m; |! P
 { s[j]=t[j-1];
7 y# G( W0 H$ m% @ s[j-1]=-p*t[j-1]+t[j-2];
9 ^4 L; L2 @! w p if (j>=3)
8 J3 R" N( k# e3 z! G$ r k4 J for (k=j-2; k>=1; k--)3 a7 r/ i! u: z' t+ [- e, g5 Z, y
 s[k]=-p*t[k]+t[k-1]-q*b[k];9 k$ \: l3 C! o* j& i
 s[0]=-p*t[0]-q*b[0];
; n" b1 _; S5 k- q8 x! u d2=0.0; c=0.0; g=0.0;
. V# s7 ^$ O4 `: v9 G for (i=0; i<=n-1; i++)
. _7 v* z5 |1 j0 C) v' ^ { q=s[j];! B9 Q& T2 W+ _2 t
 for (k=j-1; k>=0; k--)
6 m. o- |: t" ]/ y; S/ U$ w q=q*(x-z)+s[k];9 c% }9 f6 ]+ g; V& ~
 d2=d2+q*q; c=c+y*q;
2 g3 N. ?: e( @6 |1 K- ]- ` g=g+(x-z)*q*q;- p( |1 k2 _) v6 @- w* R H/ @
 }0 R, @/ x' k) k: [% p
 c=c/d2; p=g/d2; q=d2/d1;: f0 m) }$ {% S2 H( h) l1 e
 d1=d2;
0 m; f. y7 a. \. |# d2 ^ a[j]=c*s[j]; t[j]=s[j];
+ V% m' L* o4 ^- J for (k=j-1; k>=0; k--)
, {1 w- O7 _4 g9 P3 X$ b5 O { a[k]=c*s[k]+a[k];$ p9 X# w \6 [) q! t+ _ m
 b[k]=t[k]; t[k]=s[k];" k, }0 m/ E2 L2 Z. S1 o
 }! V8 K1 K& `3 P" }+ ?
 }0 e4 A ]6 s ~, N" ^
dt[0]=0.0; dt[1]=0.0; dt[2]=0.0;
8 g: Q. k6 |2 m' x7 t. \( k" a for (i=0; i<=n-1; i++), F9 l% z' e/ _2 k9 ^$ n
 { q=a[m-1];
% S. |4 k/ L* Z7 D: M. z for (k=m-2; k>=0; k--)
3 [! c) c' o' l9 ` q=a[k]+q*(x-z);1 U w$ Z& _: N1 s
 p=q-y;& u! c$ M. g! X) t- x7 M
 if (fabs(p)>dt[2]) dt[2]=fabs(p);' ?& I* S5 E2 q. E% v% Z3 j3 } Z
 dt[0]=dt[0]+p*p;
, L4 [0 N+ F# k) E7 a) Q dt[1]=dt[1]+fabs(p);
, u' X n0 J8 m' T8 `/ g. I }
: Y+ w3 f" R4 J/ ~7 s! U" M return;0 n+ Q3 d, ^( \8 ~# A7 L3 P- _
 }<>龙贝格积分法<>#include "math.h"
3 N7 K: f/ g% f3 h G double fromb(a,b,eps)
& u! V3 C. \, M double a,b,eps;
% C4 y- u# \6 ~ y! k* `: ]! ^ { extern double frombf();
* ^7 c5 a+ G$ T- {0 e- a4 j) @ int m,n,i,k;. Q' ` H* d( n
double y[10],h,ep,p,x,s,q;% u, B6 t4 |- m1 h; y% F' j+ \2 t
h=b-a;
. C" O% e5 E; S* y; U y[0]=h*(frombf(a)+frombf(b))/2.0;. {3 e- f) W' n4 d8 z
m=1; n=1; ep=eps+1.0;
/ v5 E6 G" V2 G while ((ep>=eps)&&(m<=9))
2 F" }! H0 ]7 J# O; k2 L3 J { p=0.0;
9 w2 c, n% K) i& T* K. G& c for (i=0;i<=n-1;i++)
; q% A) f/ a" i/ H { x=a+(i+0.5)*h;, i% k) V* |# p ]- n& i
 p=p+frombf(x);
" ]* ?% R# j; ~0 D& p2 l ` }' n5 O3 K* u1 B ~, e! {3 w! q8 b
 p=(y[0]+h*p)/2.0;0 l4 T$ A" v8 J: L# x N1 S
 s=1.0;* Q1 }" W4 b/ R
 for (k=1;k<=m;k++)8 o# d" [8 {- x
 { s=4.0*s;' h4 \ y$ H0 j
 q=(s*p-y[k-1])/(s-1.0);
$ L: Y$ m% z/ Q5 L* P- C y[k-1]=p; p=q;
' [4 T* o+ i% _. x) i }
/ U8 p9 Z8 b* m* M+ G; Y1 W( Z! A n ep=fabs(q-y[m-1]);
8 n# N" M @: c W6 ] Y( s: T m=m+1; y[m-1]=q; n=n+n; h=h/2.0;
7 K0 l" t# ?7 r9 u9 |7 K }
' j- `& p' O2 N, N% f$ z0 g- R9 {# s return(q);& l X; W; a" T4 h8 h# i
 }<>呵呵希望对你有用！！

我也说一句

发表

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

使用道具举报

返回列表

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

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

？

浏览过的版块

新人进步奖

元老勋章

QQ

电话咨询

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

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

？

浏览过的版块

新人进步奖

元老勋章

QQ

电话咨询

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

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