二維插值算法與實現(xiàn)演示文稿

二維插值算法與實現(xiàn)演示文稿目前一頁\總數(shù)五十一頁\編于十六點(diǎn)第一個例題三角形的面積求解2目前二頁\總數(shù)五十一頁\編于十六點(diǎn)三角形面積求解可以采用一個界面的輸入方式，也可采用inputbox的輸入方式，但要注意，因為輸入的次數(shù)為三次，要用定數(shù)循環(huán)，不能沒有條件的循環(huán)。3目前三頁\總數(shù)五十一頁\編于十六點(diǎn)三角形面積求解PrivateSubForm_Load()Text2=""EndSubPrivateSubText1_Change(IndexAsInteger)IfNotIsNumeric(Text1(Index).Text)ThenText2.Text="輸入的不是數(shù)字，請重新輸入"Text1(Index).SetFocusText1(Index).Text=""ElseText2="輸入的數(shù)字是正確的"EndIfEndSub初始化判斷輸入的數(shù)據(jù)4目前四頁\總數(shù)五十一頁\編于十六點(diǎn)三角形面積求解PrivateSubCommand1_Click()‘計算DimA#,B#,C#,S#,Mianji#A=Val(Text1(0).Text)B=Val(Text1(1).Text)C=Val(Text1(2).Text)IfA+B>CAndB+C>AAndA+C>BThenS=(A+B+C)/2Mianji=Sqr(S*(S-A)*(S-B)*(S-C))Text2.Text=Format(Mianji,"三角形的面積為0.###")ElseText2.Text="輸入的三條邊不能構(gòu)成三角形"EndIfEndSub5目前五頁\總數(shù)五十一頁\編于十六點(diǎn)三角形面積求解方法26目前六頁\總數(shù)五十一頁\編于十六點(diǎn)三角形面積求解PrivateSubCommand4_Click()OnErrorResumeNextDimAA()AsStringIfText3(0).Text<>""OrText3(1).Text<>""ThenIfOption1(0).Value=TrueThenAA=Split(Text3(0).Text,",")IfOption1(1).Value=TrueThenAA=Split(Text3(1).Text,"")ElseText4="沒有數(shù)據(jù)輸入"EndIf下續(xù)7目前七頁\總數(shù)五十一頁\編于十六點(diǎn)三角形面積求解FunctionSplit(ExpressionAsString,[Delimiter],[LimitAsLong=-1],[CompareAsVbCompareMethod=vbBinaryCompare])VBA.Strings的成員Splitastringintoanarray解釋：split把一個字符串按分隔符分成一個數(shù)組，這個數(shù)組須是動態(tài)數(shù)組?？崭穹指羧菀壮霈F(xiàn)問題，但逗號分隔不會出現(xiàn)問題。例如：8目前八頁\總數(shù)五十一頁\編于十六點(diǎn)三角形面積求解多加幾個空格沒有問題，但不能有兩個逗號連著多加幾個空格出現(xiàn)問題，分解的值有問題9目前九頁\總數(shù)五十一頁\編于十六點(diǎn)代碼續(xù)DimA#,B#,C#,S#,Mianji#A=Val(AA(0))B=Val(AA(1))C=Val(AA(2))IfA+B>CAndB+C>AAndA+C>BThenS=(A+B+C)/2Mianji=Sqr(S*(S-A)*(S-B)*(S-C))Text4.Text="三角形的三邊分別為a="&A&",b="&B&",c="&C&vbCrLf&Format(Mianji,"三角形的面積為0.###")ElseText4.Text="三角形的三邊分別為a="&A&",b="&B&",c="&C&vbCrLf&"輸入的三條邊不能構(gòu)成三角形"EndIfEndSub10目前十頁\總數(shù)五十一頁\編于十六點(diǎn)103.75,-1,38,7,15,1;96.75,1,86,52,30,1;420.77,-1,39,18,30,1;208.64,1,62,31,30,1;436.79,1,25,56,30,-1;56.48,-1,87,1,15,-1;"98.8,-1,53,22,0,-1;"68.62,1,36,38,0,-1;"95.1,1,59,29,0,-1;""207.41,-1,51,40,45,-1"11目前十一頁\總數(shù)五十一頁\編于十六點(diǎn)PublicSubStrtoShuzhu(ByValStr1AsString,A()AsDouble)DimRowDataAsVariantDimColDataAsVariantDimRowtmpDatAsVariantDimColTmpDatAsVariantDimrAsInteger,CAsIntegerRowData=Split(Str1,";")ReDimA(UBound(RowData),1)r=-112目前十二頁\總數(shù)五十一頁\編于十六點(diǎn)ForEachRowtmpDatInRowDataColData=Split(RowtmpDat,",")ReDimPreserveA(UBound(RowData),UBound(ColData))r=r+1C=-1ForEachColTmpDatInColDataC=C+1A(r,C)=ColTmpDatNextNext13目前十三頁\總數(shù)五十一頁\編于十六點(diǎn)第二個例題含水率求解14目前十四頁\總數(shù)五十一頁\編于十六點(diǎn)含水率試驗要求兩組計算一個平均值，當(dāng)兩個值求平均時，差值大于規(guī)范要求時，用紅色背景顯示該數(shù)據(jù)15目前十五頁\總數(shù)五十一頁\編于十六點(diǎn)含水率試驗PrivateSubKJSZ(ByValShuAsInteger)OnErrorResumeNextDimiAsIntegerIfShu>=10ThenText55.BackColor=vbRedText55="加載的試驗組數(shù)過多，本程序暫時不能處理"ExitSubEndIfIfShu>1AndShu<10ThenFori=2ToShuLoadHSBH(i-1)HSBH(i-1).Visible=TrueHSBH(i-1).Top=HSBH(i-2).Top+350

加載控件數(shù)組16目前十六頁\總數(shù)五十一頁\編于十六點(diǎn)含水率試驗LoadHShh(i-1)HShh(i-1).Visible=TrueHShh(i-1).Top=HShh(i-2).Top+350LoadHSmh(i-1)HSmh(i-1).Visible=TrueHSmh(i-1).Top=HSmh(i-2).Top+350LoadHSmhst(i-1)HSmhst(i-1).Visible=TrueHSmhst(i-1).Top=HSmhst(i-2).Top+350LoadHSmhgt(i-1)HSmhgt(i-1).Visible=TrueHSmhgt(i-1).Top=HSmhgt(i-2).Top+350

LoadHSmst(i-1)HSmst(i-1).Visible=True

'HSmst(i-2).Enabled=FalseHSmst(i-1).Top=HSmst(i-2).Top+350

LoadHSmgt(i-1)HSmgt(i-1).Visible=True

'HSmgt(i-2).Enabled=FalseHSmgt(i-1).Top=HSmgt(i-2).Top+350

LoadHSw(i-1)HSw(i-1).Visible=True

'HSw(i-2).Enabled=FalseHSw(i-1).Top=HSw(i-2).Top+350

LoadHSpjw(i-1)HSpjw(i-1).Visible=True

'HSpjw(i-2).Enabled=FalseHSpjw(i-1).Top=HSpjw(i-2).Top+350NextiEndIfHSQK_ClickEndSub17目前十七頁\總數(shù)五十一頁\編于十六點(diǎn)含水率試驗PrivateSubCommand16_Click()OnErrorResumeNextDimiAsIntegerFori=HSBH.UBoundTo1Step-1UnloadHSBH(i)UnloadHShh(i)UnloadHSmh(i)UnloadHSmhst(i)UnloadHSmhgt(i)UnloadHSmst(i)UnloadHSmgt(i)UnloadHSw(i)UnloadHSpjw(i)Nexti控件是不能重復(fù)加載的，如需重復(fù)加載，可以用錯誤捕捉也可以先缷載要加載的控件18目前十八頁\總數(shù)五十一頁\編于十六點(diǎn)含水率試驗KJSZ(Val(Text54))'UShu=Val(Text54)Label66.ForeColor=vbGreenLabel66.FontBold=TrueLabel66.Caption="計算成功，請保存數(shù)據(jù)"ExitSuberr1:Label66.ForeColor=vbRedLabel66.FontBold=TrueLabel66.Caption="計算出錯了,請確認(rèn)是否輸入數(shù)據(jù)有誤，注意試樣體積和含水率的輸入"EndSub函數(shù)調(diào)用，有兩種方法19目前十九頁\總數(shù)五十一頁\編于十六點(diǎn)含水率試驗PrivateSubHSJS_Click()OnErrorGoToerr1DimiAsIntegerIfVal(Text54)>1ThenFori=0ToVal(Text54)-1HSmst(i)=Format(Val(HSmhst(i))-Val(HSmh(i)),"#####0.000")HSmgt(i)=Format(HSmhgt(i)-HSmh(i),"#####0.000")HSw(i)=Format((HSmhst(i)-HSmhgt(i))*100/HSmgt(i),"#####0.000")NextFori=1ToVal(Text54)IfiMod2=0ThenIfAbs(Val(HSw(i-1))-Val(HSw(i-2)))<2ThenHSpjw(i-2).ForeColor=vbBlackElseHSpjw(i-2).BackColor=vbRedHSpjw(i-2).FontBold=TrueEndIf

HSpjw(i-2).Text=(Val(HSw(i-1))+Val(HSw(i-2)))/2UnloadHSpjw(i-1)

EndIfNext計算含水率的程序濕土計算干土計算含水率計算含水率差值計算20目前二十頁\總數(shù)五十一頁\編于十六點(diǎn)含水率試驗Elsei=0HSmst(i)=Format(Val(HSmhst(i))-Val(HSmh(i)),"#####0.000")HSmgt(i)=Format(HSmhgt(i)-HSmh(i),"#####0.000")HSw(i)=Format((HSmhst(i)-HSmhgt(i))*100/HSmgt(i),"#####0.000")HSpjw(i)=HSw(i)EndIfLabel66.ForeColor=vbGreenLabel66.FontBold=TrueLabel66.Caption="計算成功，請保存數(shù)據(jù)"ExitSuberr1:Label66.ForeColor=vbRedLabel66.FontBold=TrueLabel66.Caption="計算出錯了,請確認(rèn)是否輸入數(shù)據(jù)有誤"EndSub提示計算成功還是失敗的信息21目前二十一頁\總數(shù)五十一頁\編于十六點(diǎn)表格法實現(xiàn)方法DimStr1AsString,AA()AsDoubleSpreadsheet1.Cells(1,1)="盒號"Spreadsheet1.Cells(1,2)="m盒(g)"Spreadsheet1.Cells(1,3)="m盒＋濕土(g)"Spreadsheet1.Cells(1,4)="m盒＋干土(g)"Spreadsheet1.Cells(1,5)="含水率(%)"Spreadsheet1.Cells(1,6)="平均含水率(%)"Str1="314,14.414,84.355,75.614;345,14.397,82.576,73.323"StrtoShuzhuStr1,AA()DimIAsInteger,JAsIntegerForI=0ToUBound(AA,1)ForJ=0ToUBound(AA,2)Spreadsheet1.Cells(I+2,J+1)=AA(I,J)'把字符串分開填入表格中，以備以后測試用

NextJNext22目前二十二頁\總數(shù)五十一頁\編于十六點(diǎn)含水率計算DimIAsInteger,JAsInteger,AA()ReDimAA(I-2,2)ForI=0ToSpreadsheet1.ActiveSheet.UsedRange.Rows.Count-2ForJ=0To2AA(I,J)=Spreadsheet1.Cells(I+2,J+2)‘讀入表格數(shù)據(jù)

NextJSpreadsheet1.Cells(I+2,5)=Format((AA(I,1)-AA(I,2))/(AA(I,2)-AA(I,0))*100,"#.000")'含水率

IfI>=1AndIMod2=1ThenSpreadsheet1.Cells(I+2,6)=(Spreadsheet1.Cells(I+2,5)+Spreadsheet1.Cells(I+1,5))/2

EndIf

Next23目前二十三頁\總數(shù)五十一頁\編于十六點(diǎn)作業(yè)1添加Microsoftofficexp控件用表格法編寫含水率或密度計算程序?qū)W習(xí)用split()和strtoshuzhu()函數(shù)把一個字符串轉(zhuǎn)變成數(shù)組。有圖的請抓圖，作業(yè)可以是手寫版也可以是打印版。24目前二十四頁\總數(shù)五十一頁\編于十六點(diǎn)密度試驗第三個例題密度求解25目前二十五頁\總數(shù)五十一頁\編于十六點(diǎn)mdmd1md2md3md4md5md6md7md826目前二十六頁\總數(shù)五十一頁\編于十六點(diǎn)加載試驗組數(shù)PrivateSubText5_Change()DimiAsIntegerIfVal(Text5)<=0orVal(Text5)>10ThenMsgBox"輸入數(shù)字出錯",vbInformation,"請重新輸入"ExitSubEndIfIIfVal(Text5)>0Then'加載試驗組數(shù)

IfVal(Text5)>MD1.CountThenFori=MD1.ubound+1ToVal(Text5)-1LoadMD1(i)：LoadMD2(i)：LoadMD3(i)：LoadMD4(i)LoadMD5(i)：LoadMD6(i)：LoadMD7(i)LoadMD8(i)：MD1(i)=i+1：MD2(i)=423+i動態(tài)加載控件時，每個控件的index屬性設(shè)置為027目前二十七頁\總數(shù)五十一頁\編于十六點(diǎn)加載控件，控件顯示位置MD1(i).Top=MD1(i-1).Top+350'加載后的控件位置

MD2(i).Top=MD2(i-1).Top+350

MD3(i).Top=MD1(i-1).Top+350MD4(i).Top=MD2(i-1).Top+350MD5(i).Top=MD1(i-1).Top+350MD6(i).Top=MD2(i-1).Top+350MD7(i).Top=MD1(i-1).Top+350MD8(i).Top=MD2(i-1).Top+35028目前二十八頁\總數(shù)五十一頁\編于十六點(diǎn)讓控件顯示在窗體上MD1(i).Visible=True'顯示加載后的控件

MD2(i).Visible=True

MD3(i).Visible=TrueMD4(i).Visible=TrueMD5(i).Visible=TrueMD6(i).Visible=TrueMD7(i).Visible=TrueMD8(i).Visible=TrueNext29目前二十九頁\總數(shù)五十一頁\編于十六點(diǎn)當(dāng)界面上的控件多時，控件的卸載ElseIfVal(Text5)<MD1.CountThen'如果加載的控件比目前顯示的少，卸載控件

Fori=MD1.uboundToVal(Text5)Step-1Ifi>0ThenUnloadMD1(i)UnloadMD2(i)：UnloadMD3(i)：UnloadMD4(i)UnloadMD5(i)：UnloadMD6(i)：UnloadMD7(i)UnloadMD8(i)EndIfNextEndIfEndIfMD1(0)=1'這些數(shù)據(jù)不重要，隨便編個號MD2(0)=423EndSub30目前三十頁\總數(shù)五十一頁\編于十六點(diǎn)2演示數(shù)據(jù)PrivateSubMDYS_Click()DimiAsIntegerText5=3’設(shè)只有三組試驗Text5_Change‘事件過程也可直接調(diào)用，加載試驗組數(shù)Fori=0To2’試驗自動編號與環(huán)刀自動編號

MD1(i)=i+1MD2(i)=423+iNextMD3(0)=44.089：MD4(0)=161.605：MD5(0)=60MD6(0)=23.01‘第一行數(shù)據(jù)MD3(1)=41.089：MD4(1)=169.605：MD5(1)=60MD6(1)=22.6MD3(2)=43.089：MD4(2)=169.607：MD5(2)=60MD6(2)=20.6EndSub31目前三十一頁\總數(shù)五十一頁\編于十六點(diǎn)計算過程1PrivateSubMDJS_Click()DimiAsInteger,KAsIntegerDimtmpAsDouble,StrtmpAsStringFori=0ToVal(Text5)-1tmp=(Val(MD4(i))-Val(MD3(i)))/Val(MD5(i))MD7(i)=Format(tmp,"0.####")‘濕密度

IfVal(MD6(i))>0Thentmp=Val(MD7(i))/(1+Val(MD6(i))/100)MD8(i)=Format(tmp,"0.####")‘干密度

EndIfNext32目前三十二頁\總數(shù)五十一頁\編于十六點(diǎn)計算過程2Strtmp="試樣編號,環(huán)刀號,刀質(zhì)量(g),環(huán)刀+濕土質(zhì)量(g),試樣體積(cm3),含水率(%),濕密度(g/cm3),干密度(g/cm3)"&vbCrLfFori=0ToVal(Text5)-1Strtmp=Strtmp&MD1(i)&","&MD2(i)&","&MD3(i)&","&MD4(i)&","&MD5(i)&","&MD6(i)&","&MD7(i)&","&MD8(i)&vbCrLfNexttmp=0Fori=0ToMD1.uboundtmp=tmp+MD7(i)Next把輸入數(shù)據(jù)和計算結(jié)果以逗號隔開，以便導(dǎo)入Excel中33目前三十三頁\總數(shù)五十一頁\編于十六點(diǎn)計算過程3tmp=Format(tmp/MD1.Count,"0.####")Strtmp=Strtmp&"-----------------"&vbCrLfStrtmp=Strtmp&"濕密度平均值為,"&tmp&"g/cm3,"tmp=0Fori=0ToMD1.ubound‘獲得控件數(shù)組的下標(biāo)最大值

tmp=tmp+MD8(i)Nexttmp=Format(tmp/MD1.Count,"0.####")Strtmp=Strtmp&"干密度平均值為,"&tmp&"g/cm3"Text70=Text68&StrtmpEndSub顯示到text70中，text70的mutiline屬性應(yīng)設(shè)置成true34目前三十四頁\總數(shù)五十一頁\編于十六點(diǎn)數(shù)據(jù)格式轉(zhuǎn)換35目前三十五頁\總數(shù)五十一頁\編于十六點(diǎn)以逗號隔開的數(shù)據(jù)分割PrivateSubCommand31_Click()OnErrorResumeNextDimAA()AsString,bb()AsString,iAsIntegerIfText6<>""Then'分割數(shù)據(jù)

AA=Split(Text6.Text,vbCrLf)‘先分行，以回車符分隔EndIfText5.Text=UBound(AA)+1：Text5_Change‘分配控件數(shù)組IfOption2(0).Value=TrueThen’逗號隔開的情況

Fori=0ToUBound(AA)IfAA(i)<>""Then‘不是空行的情況

bb=Split(AA(i),",")MD1(i)=bb(0)：MD2(i)=bb(1)：MD3(i)=bb(2)MD4(i)=bb(3)：MD5(i)=bb(4)：MD6(i)=bb(5)Else’是空行的情況

Text5=Text5.Text-1Text5_ChangeEndIfNext36目前三十六頁\總數(shù)五十一頁\編于十六點(diǎn)其它分隔符情況可以是單空格分隔，也可以是Excel數(shù)據(jù)復(fù)制過來的，其中Excel復(fù)制的分隔符是chr(9),是個特殊字符。其余的請仿此程序進(jìn)行設(shè)計。37目前三十七頁\總數(shù)五十一頁\編于十六點(diǎn)第五個問題線性插值計算引言一維線性插值原理及程序設(shè)計二維線性插值原理及程序設(shè)計結(jié)論38目前三十八頁\總數(shù)五十一頁\編于十六點(diǎn)一、引言如表1所示，當(dāng)勘察出來的e指標(biāo)為0.75，IL指標(biāo)為0.625時，且有很多組數(shù)據(jù)需要查表插值時，則計算量較大。在列表中，e最小值為0.5,如果是0.45，則計算者更難著手計算。表1粘性土承載力基本值f0

ILe0.50.750.72652400.822020039目前三十九頁\總數(shù)五十一頁\編于十六點(diǎn)二、一維線性插值原理及程序設(shè)計40目前四十頁\總數(shù)五十一頁\編于十六點(diǎn)2.1一維線性插值原理(1)情況１（e1≤e3≤e2內(nèi)部插值情況）從圖1(a)中可得： △ABC≌△ADE

當(dāng)f2≥f1時：

當(dāng)f2<f1時，同理，作出相應(yīng)圖可得：

(1)(2)41目前四十一頁\總數(shù)五十一頁\編于十六點(diǎn)2.1一維線性插值原理情況2（e3<e1外部插值情況）從圖1(b)中可得：△ABC≌△ADE當(dāng)f2≥f1時：

當(dāng)f2<f1時，同理，作出相應(yīng)圖可得：

(3)(4)42目前四十二頁\總數(shù)五十一頁\編于十六點(diǎn)2.1一維線性插值原理第三種情況的外插情況，作出相應(yīng)圖形，同理可得(5)(6)43目前四十三頁\總數(shù)五十一頁\編于十六點(diǎn)2.2線性插值計算的流程圖44目前四十四頁\總數(shù)五十一頁\編于十六點(diǎn)2.3一維線性插值程序PrivateFunctionChazhi(ByVale1AsDouble,ByValf1AsDouble,ByVale2AsDouble,ByValf2AsDouble,ByVale3AsDouble)’注：(e1,f1)與（e2,f2）已知，求e3對應(yīng)值f3DimTemp1AsDoubleIfe1>=e2Then

Temp1=e1e1=e2e2=Temp1Temp1=f1f1=f2f2=Temp1EndIf‘注：設(shè)e2>e1，否則e1與e2交換，f1與f2交換45目前四十五頁\總數(shù)五十一頁\編于十六點(diǎn)2.3一維線性插值程序Ife3>=e1Ande3<=e2Then

Iff2>=f1ThenChazhi=(e3-e1)*(f2-f1)/(e2-e1)+f1ElseChazhi=(e2-e3)*(f1-f2)/(e2-e1)+f2EndIfElseIfe3<e1Then

Iff2>=f1ThenChazhi=((e1-e3)*f2-(e2-e3)*f1)/(e1-e2)

ElseChazhi=(f1-f2)*(e2-e3)/(e2-e1)+f2

EndIf‘注：內(nèi)部插值計算,情況1‘注：f2>=f1的計算方法(1)’

注：f1>＝f2的計算方法(2)’‘注：外部插值計算，情況2’注：計算方法(3)’注：計算方法(4)46目前