【正文】 c=c(y[i1]y[i2])/(x[i1]x[i2])。 則： ZA=(1p)(1q)Zi,j+(1p)qZi,j+1+p(1q)Zi+1,j+pqZi+1,j+1 —————————— 二元線性插值綜合式 ? （ 2） 非 正交網(wǎng)格的二元線性插值 Y X Yj+1 Yj Xi,j Xi+1,j J+1 J 1 2 YA XA A X1 X2 Zi,j Zi+1,j Zi,j+1 Zi+1,j+1 Z1 Z2 ? 插值計(jì)算方法 （ 1）沿斜線單方向插值 Z1=(Zi,j+1Zi,j)/(Yj+1Yj)*(Xi,j+1Xi,j)+Zi,j Z2=(Zi+1,j+1Zi,j+1)/((Yj+1Yj)*(Xi,j+1Xi,j)+Zi,j X1=(Xi,j+1Xi,j)/(Yj+1Yj)*(Xi,j+1Xi,j)+Xi,j X2=(Xi+1,j+1Xi,j+1)/((Yj+1Yj)*(Xi,j+1Xi,j)+Xi,j (2) 沿水平單方向插值 ZA=(Z2Z1)/(X2X1)*(XAX1)+X1 ? 二元線性插值程序的設(shè)計(jì) 用二維數(shù)組如數(shù)法存儲(chǔ)二維數(shù)表或二元曲線：X(I),y(j),z(I,j) 給定插點(diǎn)坐標(biāo) :Xa,ya 用 X(I)XaX(I+1)判斷插點(diǎn)相鄰序號(hào) i0,i0+1。 擬合式： ? ?? ??mjnijiij yxaz0 1? 1）直接擬合：在曲線上直接取點(diǎn)建立擬合式 。PESET(X,Y),COLOR在 X,Y位置畫一點(diǎn) 39。 : 39。 ? A$=D20R30U30L20:DRAW A$ ? C語言基本繪圖指令 1 : 39。RECTANGLE(INT X1,INT Y1,INT X2,INT Y2)。CHAR S[10] 定義一字符串 ? 39。i— 各分割曲線序號(hào)； J— 曲線點(diǎn)序號(hào)； 3) 把數(shù)組 x(I,j),y(I,j)輸入計(jì)算機(jī)。 /*輸入圓參數(shù) :園心與半徑 */ int xc[4]={100,100,200,200}。 /*輸入 直線參數(shù) :起點(diǎn)與終點(diǎn)坐標(biāo) */ int x1[9]={80,100,140,140,280,280,170,70,40}。 cleardevice()。i=8。 20 FOR I = 1 TO 4: READ XC(I), YC(I), RC(I): NEXT I 30 DATA 100,200,10,100,200,20,200,200,10,200,200,20 40 39。 220 39。 300 FOR I = 1 TO 9: LINE (X1(I), Y1(I))(X2(I), Y2(I)): NEXT I 310 39。輸入水平閥門及其坐標(biāo) 39。屏幕顯示圖形 39。畫豎向閥門 39。 } void famen2(x1,y1) /*水平閥門定義 */ { moveto(x1,y1)。 linerel(15,15)。 moveto(x4,y4)。linerel(60,0)。linerel(90,0)。linerel(3,8)。 /*畫水平閥門 2*/ famen2(350,250)。 ? 通用繪圖程序的設(shè)計(jì) 通用數(shù)據(jù)接口 ： 曲線條數(shù) m。 /*圖形初始化 */ initgraph(graphdriver,graphmode,)。ymax=y[0]。i=m。 } for(i=0。 } for(j=0。 outtextxy(dl80,hy5,s)。 if(yy0) yy=(int)yy/。 lx=lo/5*i+dl。 /繪制坐標(biāo)軸線 / rectangle(dl,dh,lo+dl,ho+dh)。 if(xmaxx[i]) xmax=x[i]。 y[i]=sin(x[i])+。 float x[31],y[31],xmin,xmax,ymin,ymax,l[31],h[31]。 /*畫管路 1*/ pipe2(180,254)。 cleardevice()。moverel(0,10)。 } pipe2(x5,y5) /*管路 2定義 */ { moveto(180,254)。linerel(60,0)。 } void qushuikou(x3,y3) /*取水口定義 */ { moveto(x3,y3)。 linerel(10,8)。 linerel(8,0)。畫濾水器 39。輸入管道 39。輸入園的參數(shù) 39。 250 FOR I = 1 TO 4: CIRCLE (XC(I), YC(I)), RC(I), , , , : NEXT 260 39。 80 FOR I = 1 TO 9: READ X1(I), Y1(I), X2(I), Y2(I): NEXT 90 DATA 80,150,100,120,100,120,140,120,140,120,140,150 100 DATA 140,150,280,150,280,150,280,200,280,200,230,200 110 DATA 170,200,130,200, 70,200, 40,200, 40,200, 40,170 180 39。i++) line(x[i],y[i],x[i+1],y[i+1])。 /*畫園弧 */ for(i=0。 /*輸入自由曲線 :x(i),y(i)*/ int x[5]={40,50,60,70,80}。 int ya[2]={200,200}。i— 各分割曲線序號(hào)； J— 曲線點(diǎn)序號(hào)； 3) 把數(shù)組 x(I,j),y(I,j)輸入計(jì)算機(jī)。在 X,Y位置顯示數(shù)字S(即 X) ? 39。以 X,Y為園心 ,R為半徑 ,XT1為起始角 ,39。LINE(INT X1,INT Y1,INT X2,INT Y2)。 B:走筆不畫線 :BL20向左空走 20單位 。CIRCLE(X,Y),R以 X,Y為園心 ,以 R為半徑畫園 39。 ? 程序設(shè)計(jì) 離散曲線： Xk,yk,zk。 ? 曲線失真 2 1 n 完全多項(xiàng)式：Y=a0+a1X+a2X2+…+anXn 3 ???nkkk xay0? （ 2）最小而乘法擬合 ? 已知點(diǎn)： xi,yi。 m=a0+a1*(zx[i2])+a2*(zx[i2])*(zx[i1])。i6。 ? 實(shí)用式： y=a0+a1(xx1)+a2(xx1)(xx2) ? 1點(diǎn) ； y1=a0+a1(x1x1)+a2(x1x1)(x1x2) ? 得： y1=a0 ? 2點(diǎn)： y2=a0+a1(x2x1)+a2(x2x1)(x2x2) ? ? 得： a1= ? 3點(diǎn)： y3=a0+a1(x3x1)+a2(x3x1)(x3x2) ? 得： a3= 112xyy??1312121313xxxxyyxxyy??????? 歸納 131212131331212110xxxxyyxxyyaxxyyaya???????????? 程序設(shè)計(jì) 1）曲線離散： x(i),y(i)。 ? scanf(%f,z)。 printf(%,e1[i][j])。j5。j++) ? fprintf(fp,%,e1[i][j])。i5。 fprintf(fp,\n)。 } fclose(fp)。 int i,j, float e1[5][5]={{.1,.2,.3,.4,.5},{.2,.3,.4,.5,.6},{.3,.4,.5,.6,.7}, {.4,.5,.6,.7,.8},{.3,.4,.5,.6,.7}}。 2) 打開文件的方式 ： fp=open(“ abc” ,使用方式）； 使用方式 :r 讀 W寫 r+ 讀，寫等 3) 寫數(shù)據(jù)： int a=5。 用數(shù)據(jù)文件方式存放數(shù)據(jù) 存（寫文件）：數(shù)據(jù) — 計(jì)算機(jī)內(nèi)存（變量） — 磁盤文件 ?。ㄗx文件）：磁盤文件 — 計(jì)算機(jī)內(nèi)存 — （變量） —— 適用大量需要成塊變更的數(shù)據(jù)，共享數(shù)據(jù)。 int c=2023, int b[4]={1,2,3,4}。 n11— n11(j)。 ? 一、一元曲線的處理 曲線的離散與存儲(chǔ) i i+1 1 n yi xi 1) 把曲線離散位 n隔離散點(diǎn)： xi,yi 2) 用兩個(gè)一維數(shù)組分別表示自變量余函數(shù)： 3) X(I)— Xi； Y（ I） — Yi 4) 查處個(gè)點(diǎn) X(I),Y(I)值 ,把數(shù)組 X(I),Y(I)輸入計(jì)算機(jī)。 2）按照對(duì)應(yīng)關(guān)系輸入兩個(gè)數(shù)組； 3）引用：給定 z1(i),檢索 K(I). ? 不規(guī)則一維數(shù)表的處理 （ 1）不規(guī)則數(shù)表：自變量或函數(shù)是一個(gè)區(qū)間的數(shù)表 轉(zhuǎn)輪直徑D1(m) ~ ~ 導(dǎo)葉數(shù)目 z0 12 16 24 32 ? （ 2） 處理法 ? 數(shù)表改造： 轉(zhuǎn)輪直徑下限 D11（ I) 0 轉(zhuǎn)輪直徑上限 D12(I) 導(dǎo)葉數(shù)目 z0（ I） 12 16 24 32 ?用兩個(gè)一維數(shù)組分別表示自變量的山下限，用兩個(gè)一維數(shù)組 表示函數(shù) D11[i]={0,}。 ２、數(shù)據(jù)引用： 按數(shù)據(jù)的函數(shù)關(guān)系輸入，用檢索或插值引用。 ? 一、技術(shù)資料的種類 及在工程中的應(yīng)用 ? 文字型 ：用文字描述的技術(shù)資料：設(shè)計(jì)規(guī)程，設(shè)計(jì)方法說明、指令性文件。 ? 三、軟件系統(tǒng) ? 操作平臺(tái)： dos,Windows\Unix 等 ? 語言環(huán)境： QB、 C、 VB、 Fortran等 ? 輔助軟件：支持繪圖儀、掃描儀、打印機(jī)的軟件 ? 四、水電站動(dòng)力設(shè)備 CAD的特點(diǎn) ? 復(fù)雜系統(tǒng)的設(shè)計(jì)：設(shè)備選擇、計(jì)算、系統(tǒng)設(shè)計(jì)、設(shè)備設(shè)置等，過程復(fù)雜。 輸入設(shè)備： 1）文字輸入：鍵盤、手寫、語音。 ? 通用、專用應(yīng)用程序的設(shè)計(jì)。設(shè)計(jì)依據(jù)（參數(shù)選擇）。 處理方法： ? 數(shù)表化 水輪機(jī)水頭 H(m) 蝸殼壓力上升允許值△H/H(%) 40 40~100 100 =70 =50 =30 水輪機(jī)水頭下限H1(I)(m) 水輪機(jī)水頭上限H2(I)(m) 蝸殼壓力上升允許值 DH(I) 0 41 101 40 100 700 數(shù)組標(biāo)準(zhǔn)化 數(shù)表輸入計(jì)算機(jī) ? 第二講 技術(shù)資料的計(jì)算機(jī)處理 — 數(shù)表、曲線的處理 一、數(shù)表的處理 數(shù)表的種類與特征 種類：函數(shù)數(shù)表、經(jīng)驗(yàn)數(shù)表、試驗(yàn)數(shù)表； 特征 1：規(guī)則數(shù)表、不規(guī)則數(shù)表。 ? 不規(guī)則二維數(shù)表的處理 ? 不規(guī)則二維數(shù)表的形式 水輪機(jī)轉(zhuǎn)子估重公式： G=krDi2Lt G轉(zhuǎn)子帶軸重量（ t)； Di轉(zhuǎn)子鐵心外徑； Lt鐵心長(zhǎng)度； Kr=f(Di,n) 發(fā)電機(jī)轉(zhuǎn)速 n(r/min) Di(m) Kr 備注 214~750 150~240 —— =5 5~9 —— ~ ~ —— 高速機(jī)組 Di小 kr取上限 ? 處理原則與方法 原則：模擬人工方法，便于計(jì)算機(jī)輸入輸出。 n11— n11(J)。 3）共享數(shù)據(jù)：不同程序共同使用的數(shù)據(jù)。 C: scanf(“ %d%f” 。 ? 標(biāo)準(zhǔn)的接口： 如上： m,n,Q11[i]，n11[j],E[i][j] 可實(shí)現(xiàn)共享的數(shù)據(jù)表達(dá)方式 ：二進(jìn)、八進(jìn)、十六進(jìn)等、 ASCII碼， ASCII碼可供不同語言的程序調(diào)用。fprintf(fp,” %c” ,c)。j5