【正文】 以 X1,Y1X2,Y2為對(duì)角線畫(huà)矩形 畫(huà)園 : CIRCLE(INT X,INT Y,INT R)。從現(xiàn)行位置到增量 DX,DY畫(huà)一線段 : 39。在 X1,Y1X2,Y2間畫(huà)一直線 39。 : 39。PUTPIXEL(INT X,INT Y,COLOR)。使用 :DRAW D20R30U20L30畫(huà)一 20*30的矩形 。BM+10,10:向 +10,10方向空走 。 R,L向右 ,左畫(huà) Mx,y:向任意方向（ x,y)畫(huà)； M+20,10:從當(dāng)前位置畫(huà)到 +20,10處 。LOCATE 行 ,列 :PRINT X:在某行 ,列處顯示 X的值 ? Basic基本繪圖指令 2 ? . DRAW指令 : U:向上畫(huà) :U20向上畫(huà) 20單位 。XT2為終止角 (弧度 )畫(huà)園弧 . 39。 : 39。 : 39。 : 39。 : 39。 選擇擬合式及擬合式次數(shù) m,n： 計(jì)算矩陣 M,T 解聯(lián)立方程組求 a00~amn 代入 x值求 y值。,...2,1,0,0 ???? ?? 以 m=n=2展開(kāi)： P=0,q=0 a00∑xk0 yk0+a10∑xk1yk0+a20∑xk2yk0+ a01∑xk0 yk1+a11∑xk1yk1+a21∑xk2yk1+ a02∑xk0 yk2+a12∑xk1yk2+a22∑xk2yk2 P=0,q=0 時(shí)000???as寫(xiě)成矩陣形式： ∑1 ∑x ∑x2 ∑y ∑xy ∑x2y ∑y2 ∑xy2 ∑x2y2 ∑x………………………………………………………… ∑x2……………………………………………………………………………………. ∑y……………………∑xy 3…………………………… …………………………………………………………… ∑x2y2 ………………………………………………………………….. ∑x2y2 M= ? T=[∑Z ∑xZ ∑X2Z ∑yZ ∑XyZ ∑X2yZ …. ∑X 2y2Z ]T a= [a00 a10 a20 a01 a11 a21 a02 a12 a22]T 解聯(lián)立方程組，求出 a00~a22,得擬合適。 以 m=n=2展開(kāi)擬合式 222221220222111012202300yxaxyayayxaxyayaxaxaaz?????????在曲線上取 9點(diǎn)建立擬合式，組成聯(lián)立方程組，解之求出待定系數(shù) a00~a22即獲得擬合式。 ? k=1 k=2 k=3 Z k Yk Xk 二元曲線（曲面）的擬合 元曲線：曲面的平面表示 —以等高線形式表達(dá)的面。 ? 選擇擬合式（次數(shù) n) ? 各型值點(diǎn)的擬合誤差： ? 誤差平方和： inkkiki yxa ?? ?? 0? ???nkkk xay0??mis12?? 最佳擬合條件 ? S— min: 0).().(.2))(()(1 01 01 021 1 022?????????????? ?? ?? ?? ? ?? ??? ?? ?? ? ?nimkjiijkiknimkjiikiknimkikikjjninimkikikixyxaxyxayxaaasyxaS即：?mjasj,...2,1,0,0 ????? ?????????????????????????????????????????????????????nijiimjimnijinijinijijniiimimniiniiniiniiimimniiniiniijxyxaxaxaxaasxyxaxaxaxaasxyxaxaxaxaasmjas1122111101011321211101101221111000....:0....:0....0,...2,1,0,0展開(kāi)展開(kāi)展開(kāi)：? 寫(xiě)成矩陣形式 :=T a=[a0 a1 a2 …a m] ???????????????????????? ? ? ?? ? ? ?? ? ? ?? ? ? ???? ? ? ??? ? ? ?ninininimimimimininininimiiiininininimiiiixxxxxxxxxxxxM1 1 1 12211 1 1 113211 1 1 1210........ ........ .......... .......... ...........???????????????????????????????????nimiiniiiniiixyxyxyT11010....... .....? 求待定系數(shù) ? 解聯(lián)立方程組，求 a0 a1 a2 …am。 2）建立聯(lián)立方程組： a0+a1X1+a2X12+…+anX1n=y1 a0+a1X2+a2X22+…+anX2n=y2 ……………………………….. a0+a1Xm+a2Xm2+…+anXmn=ym 3) 解聯(lián)立方程組求待定系數(shù) a0~an 缺點(diǎn)：點(diǎn)數(shù)多時(shí)，失真。 用 y(j)yay(j+1)判斷插點(diǎn)相鄰序號(hào) j0,j0+1。 q= (Yj+1Yj)*(YaYj)。 printf(開(kāi)度 a0=%f(mm)\n飛逸轉(zhuǎn)速 =%f(rpm),z,m)。 a0=y[i2]。 a2=c/(x[i]x[i1])。 c=(y[i]y[i2])/(x[i]x[i2])。i++) if(z=x[i]) break。 for(i=1。 printf(請(qǐng)輸入開(kāi)度 a0的值 \n)。 float z,a1,a2,a3,m,x[6]={0,10,20,30,40,50}。 數(shù)據(jù)輸入計(jì)算機(jī)； 2）給出插值變量 x。 ? 缺點(diǎn)：需要解聯(lián)立方程組。 ? printf(\n進(jìn)口斷面流速為 \n%,a)。i++) ? if(z=x[i]) a=y[i1]+(y[i]y[i1])/(x[i]x[i1])。 ? for(i=1。 ? printf(請(qǐng)輸入設(shè)計(jì)水頭 hr的值 \n)。 ? float z,a。 數(shù)據(jù)輸入計(jì)算機(jī)； 2）給出插值變量 x。}} fclose(fp)。j++) {fscanf(fp,%,e1[i][j])。i++) for(j=0。} for(i=0。j++) {fscanf(fp,“ %5d” ,q1[j])。} for(j=0。i++) { fscanf(fp,“ %5d” ,n1[i])。 for (i=0。 ? fclose(fp)。j5。i5。j++) ? fprintf(fp,%5d,q1[j])。i++) ? for (j=0。 ? for (i=0。 ? 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}}。 ? int i,j,n1[5]={50,60,70,80,90}。 } fclose(fp)。j++) fscanff(fp,%,e1[i][j])。i++) {for(j=0。 for(i=0。 ? 讀數(shù)據(jù)程序 include main() {FILE *fp。 fprintf(fp,\n)。j5。i5。 fp=fopen(,w)。 ? 寫(xiě)數(shù)據(jù)程序 include main() {FILE *fp。fprintf(fp,” %c” ,c)。 fprintf(fp,” %f” ,b)。 fprintf(fp,” %d” ,a)。 輸出文件緩沖區(qū) 輸入文件緩沖區(qū) 磁盤(pán)文件 ? 基本指令 1)文件指針變量 ： FILE *fp。 ? 標(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碼可供不同語(yǔ)言的程序調(diào)用。 N11[]={60， 70， 80， ……..120}。 ? 三、數(shù)據(jù)存取的基本要求 統(tǒng)一的數(shù)據(jù)結(jié)構(gòu)形式： i=1 2 3 ……………..m J=n J=2 J=1 m=15?！?少量會(huì)話(huà)數(shù)據(jù)。 C: scanf(“ %d%f” 。 鍵盤(pán)輸入： basic:input “ x” 。—— 適用于少量數(shù)據(jù)。 ? 二、數(shù)據(jù)的存放方式 直接為變量賦值 ： a=。 3）共享數(shù)據(jù)：不同程序共同使用的數(shù)據(jù)。 數(shù)學(xué)特征： 單數(shù)據(jù)、一維數(shù)組、二維數(shù)組、多維數(shù)組 應(yīng)用特征： 1）公用數(shù)據(jù)：程序中不變更的數(shù)據(jù)，如、水輪機(jī)選型程序中的型譜參數(shù)。 EM(I,j)— η i,j 3）把三個(gè)數(shù)組輸入計(jì)算機(jī)； 4）用二元插值引用曲線 ? 第四講 數(shù)據(jù)的計(jì)算機(jī)存取 講授內(nèi)容：計(jì)算機(jī)數(shù)據(jù)存儲(chǔ)的基本方法及程序設(shè)計(jì) 講授時(shí)間： 2學(xué)時(shí) 講授方式：課堂教學(xué) 重點(diǎn)：數(shù)據(jù)文件的操作 學(xué)習(xí)要求：掌握數(shù)據(jù)存取程序的設(shè)計(jì) ? 一、數(shù)據(jù)的作用及類(lèi)型 作用： 為程序中的變量賦值。 ? 非正交網(wǎng)格離散法 Φ=50 Φ=00 Φ=+50 Φ=+100 Q11(I,j) N11(j) η i,j ? 非正交網(wǎng)格離散法 1）用一非正交網(wǎng)格覆蓋曲線，用網(wǎng)格的節(jié)點(diǎn)離散曲線； 2）用兩個(gè)一維數(shù)組何一個(gè)二維數(shù)組表示曲線坐標(biāo)參數(shù)，用一個(gè)二維數(shù)組表示函數(shù)： Q11— Q11(I， j)。 n11— n11(J)。 直接離散法的優(yōu)缺點(diǎn) 優(yōu)點(diǎn)：保持圖形幾何特征，便于圖形復(fù)原； 缺點(diǎn)：插值計(jì)算復(fù)雜。 5) 曲線的引用：用一元插值計(jì)算引用曲線，給定 x,求 y. ? 二、二元曲線的離散處理 ? 二元曲線的特點(diǎn)：曲面的平面表示 n11 Q11 η=f(Q11,n11) ? 直接離散 :在曲線上直接取點(diǎn)離散 i=1 i=2 i=3 J=1 J=2 J=3 η i,j ni,j Qi,j ? 離散步驟： 1）在曲線上直接取點(diǎn)離散 2）等值線序號(hào) i=1,2,3,…m。 處理原則： 離散化（數(shù)字化）、便于存儲(chǔ)，便于引用。 ? 不規(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ī)輸入輸出。{,},…}。 D