【正文】 outtextxy(dl80,hy5,s)。 if(yy0) yy=(int)yy/。 } for(j=0。 lx=lo/5*i+dl。 } for(i=0。 /繪制坐標(biāo)軸線 / rectangle(dl,dh,lo+dl,ho+dh)。i=m。 if(xmaxx[i]) xmax=x[i]。ymax=y[0]。 y[i]=sin(x[i])+。 /*圖形初始化 */ initgraph(graphdriver,graphmode,)。 float x[31],y[31],xmin,xmax,ymin,ymax,l[31],h[31]。 ? 通用繪圖程序的設(shè)計(jì) 通用數(shù)據(jù)接口 ： 曲線條數(shù) m。 /*畫管路 1*/ pipe2(180,254)。 /*畫水平閥門 2*/ famen2(350,250)。 cleardevice()。linerel(3,8)。moverel(0,10)。linerel(90,0)。 } pipe2(x5,y5) /*管路 2定義 */ { moveto(180,254)。linerel(60,0)。linerel(60,0)。 moveto(x4,y4)。 } void qushuikou(x3,y3) /*取水口定義 */ { moveto(x3,y3)。 linerel(15,15)。 linerel(10,8)。 } void famen2(x1,y1) /*水平閥門定義 */ { moveto(x1,y1)。 linerel(8,0)。畫豎向閥門 39。畫濾水器 39。屏幕顯示圖形 39。輸入管道 39。輸入水平閥門及其坐標(biāo) 39。輸入園的參數(shù) 39。 300 FOR I = 1 TO 9: LINE (X1(I), Y1(I))(X2(I), Y2(I)): NEXT I 310 39。 250 FOR I = 1 TO 4: CIRCLE (XC(I), YC(I)), RC(I), , , , : NEXT 260 39。 220 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。 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。i++) line(x[i],y[i],x[i+1],y[i+1])。i=8。 /*畫園弧 */ for(i=0。 cleardevice()。 /*輸入自由曲線 :x(i),y(i)*/ int x[5]={40,50,60,70,80}。 /*輸入 直線參數(shù) :起點(diǎn)與終點(diǎn)坐標(biāo) */ int x1[9]={80,100,140,140,280,280,170,70,40}。 int ya[2]={200,200}。 /*輸入圓參數(shù) :園心與半徑 */ int xc[4]={100,100,200,200}。i— 各分割曲線序號(hào)； J— 曲線點(diǎn)序號(hào)； 3) 把數(shù)組 x(I,j),y(I,j)輸入計(jì)算機(jī)。i— 各分割曲線序號(hào)； J— 曲線點(diǎn)序號(hào)； 3) 把數(shù)組 x(I,j),y(I,j)輸入計(jì)算機(jī)。在 X,Y位置顯示數(shù)字S(即 X) ? 39。CHAR S[10] 定義一字符串 ? 39。以 X,Y為園心 ,R為半徑 ,XT1為起始角 ,39。RECTANGLE(INT X1,INT Y1,INT X2,INT Y2)。LINE(INT X1,INT Y1,INT X2,INT Y2)。 ? A$=D20R30U30L20:DRAW A$ ? C語言基本繪圖指令 1 : 39。 B:走筆不畫線 :BL20向左空走 20單位 。 : 39。CIRCLE(X,Y),R以 X,Y為園心 ,以 R為半徑畫園 39。PESET(X,Y),COLOR在 X,Y位置畫一點(diǎn) 39。 ? 程序設(shè)計(jì) 離散曲線： Xk,yk,zk。 擬合式： ? ?? ??mjnijiij yxaz0 1? 1）直接擬合：在曲線上直接取點(diǎn)建立擬合式 。 ? 曲線失真 2 1 n 完全多項(xiàng)式：Y=a0+a1X+a2X2+…+anXn 3 ???nkkk xay0? （ 2）最小而乘法擬合 ? 已知點(diǎn)： xi,yi。 則： 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。 m=a0+a1*(zx[i2])+a2*(zx[i2])*(zx[i1])。 c=c(y[i1]y[i2])/(x[i1]x[i2])。i6。 float y[6]={0,100,116,123,128}。 ? 實(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)。 ? printf(設(shè)計(jì)水頭 hr為 \n%,z)。 ? scanf(%f,z)。 3）判斷插點(diǎn)位置： k,k+1 4）應(yīng)用插值公式： y=(yk+1yk)/(Xk+1Xk)*(XXk)+YK ? 程序例： ? /*一元插值以計(jì)算混凝土蝸殼進(jìn)口斷面流速為例 */ ? include ? main() ? {int i,x[6]={10,20,30,40,50,60}。 printf(%,e1[i][j])。i5。j5。i5。j++) ? fprintf(fp,%,e1[i][j])。 ? for(i=0。i5。 ? int q1[5]={500,600,700,800,900}。 fprintf(fp,\n)。i5。 } fclose(fp)。i++) {for(j=0。 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}}。 char c=“ word” 。 2) 打開文件的方式 ： fp=open(“ abc” ,使用方式）； 使用方式 :r 讀 W寫 r+ 讀，寫等 3) 寫數(shù)據(jù)： int a=5。 E[][]={,….},{,},…}。 用數(shù)據(jù)文件方式存放數(shù)據(jù) 存（寫文件）：數(shù)據(jù) — 計(jì)算機(jī)內(nèi)存（變量） — 磁盤文件 ?。ㄗx文件）：磁盤文件 — 計(jì)算機(jī)內(nèi)存 — （變量） —— 適用大量需要成塊變更的數(shù)據(jù)，共享數(shù)據(jù)。x。 int c=2023, int b[4]={1,2,3,4}。 2）可替換數(shù)據(jù)：通用程序中針對(duì)不同情況要更換的數(shù)據(jù)，如水輪機(jī)特性曲線計(jì)算程序中的模型曲線參數(shù)。 n11— n11(j)。 ? 間接里散法 n11(j) Q11(I) η i,j ? 正交網(wǎng)格離散法 1）用一正交網(wǎng)格覆蓋曲線，用網(wǎng)格的 節(jié)點(diǎn)離散曲線； 2）用兩個(gè)一維數(shù)組表示曲線坐標(biāo)參數(shù)，用一個(gè)二維數(shù)組表示函數(shù)： Q11— Q11(I)。 ? 一、一元曲線的處理 曲線的離散與存儲(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) 把三個(gè)數(shù)組 Q11(i)、 n11(J)、 Em(I,j)按照對(duì)應(yīng)關(guān)系輸入計(jì)算機(jī)； 3) 引用：用二元插值計(jì)算法應(yīng)用。 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,}。 4) 數(shù)組輸入計(jì)算機(jī) 5）引用：給出 T(I), 檢索 Uf1(I),Uf2(I) ? 其他例子：水輪機(jī)調(diào)節(jié)保證計(jì)算 水輪機(jī)蝸殼壓力上升允許值：水輪機(jī)蝸殼壓力上升最大允許值：額定水頭在40米以下時(shí)不得大于 70%，額定水頭在 40~100米時(shí)不得大于 50%，額定水頭在 100米以上時(shí)不得大于 30%。 ２、數(shù)據(jù)引用： 按數(shù)據(jù)的函數(shù)關(guān)系輸入，用檢索或插值引用。來源于函數(shù)、統(tǒng)計(jì)（經(jīng)驗(yàn)）、試驗(yàn)。 ? 一、技術(shù)資料的種類 及在工程中的應(yīng)用 ? 文字型 ：用文字描述的技術(shù)資料：設(shè)計(jì)規(guī)程，設(shè)計(jì)方法說明、指令性文件。 ? 人的設(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ì)