簡支梁及有限元分析過程-資料下載頁
2025-08-23 12:11本頁面
【正文】 includeincludedefine f(t) (1*pow(t,2)15*t+56)define eps pow(10,6)void sb(double *a,double *b){double t0,t1,t,h,alpha,f0,f1。int k=0。printf(請輸入初始點(diǎn)t0=)。scanf(%lf,amp。t0)。printf(\n請輸入初始步長h=)。scanf(%lf,amp。h)。printf(\n請輸入加步系數(shù)alpha(需大于1)=)。scanf(%lf,amp。alpha)。f0=f(t0)。t1=t0+h。f1=f(t1)。while(1){ printf(\nf1=%lf,f2=%lf,t0=%lf,t=%lf,h=%lf,k=%d,f0,f1,t0,t1,h,k)。 if(f1f0) { h=alpha*h。 t=t0。 t0=t1。 f0=f1。 k++。 } else { if(k==0) {h=h。t=t1。} else { *a=tt1?t:t1。 *b=tt1?t:t1。 break。 } } t1=t0+h。 f1=f(t1)。}}double hjfg(){double beta,t1,t2,t。double f1,f2。double a=0,b=0。double *c,*d。int k=0。c=amp。a,d=amp。b。sb(c,d)。printf(\n[a,b]=[%lf,%lf],a,b)。beta=(sqrt(5))/2。t2=a+beta*(ba)。f2=f(t2)。t1=a+bt2。f1=f(t1)。while(1){ printf(\n第%d次迭代的過程如下:,k+1)。 printf(\n[t1,t2]=[%lf,%lf],t1,t2)。 if(fabs(t1t2)eps) break。 else { if(f1f2) { t=(t1+t2)/2。 b=t2。 t2=t1。 f2=f1。 t1=a+bt2。 f1=f(t1)。 } else { a=t1。 t1=t2。 f1=f2。 t2=a+beta*(ba)。 f2=f(t2)。 } } k++。}t=(t1+t2)/2。printf(\nt=%lf,t)。return t。}main(){double t。 t=hjfg()。printf(\n函數(shù)的最優(yōu)值為f(%lf)=%lf,t,f(t))。} 50
電大資料相關(guān)推薦
鄂ICP備17016276號-1