用Matlab實(shí)現(xiàn)直線插補(bǔ)計(jì)算程序_第1頁
用Matlab實(shí)現(xiàn)直線插補(bǔ)計(jì)算程序_第2頁
用Matlab實(shí)現(xiàn)直線插補(bǔ)計(jì)算程序_第3頁
用Matlab實(shí)現(xiàn)直線插補(bǔ)計(jì)算程序_第4頁
用Matlab實(shí)現(xiàn)直線插補(bǔ)計(jì)算程序_第5頁
全文預(yù)覽已結(jié)束

下載本文檔

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡介

1、姓名：學(xué)號： 班級：用 Matlab 實(shí)現(xiàn)直線插補(bǔ)計(jì)算程序clear; F=0;x=0;y=0; dx=0.5;dy=0.6;xe=-5;ye=-6; K=ye/xe;xx(1)=x;yy(1)=y; index=2;while (1)if (F>=0) if (xe>0) x=x+dx; else x=x-dx;end y=y;F=abs(xe*y)-abs(ye*x);else if (ye>0) y=y+dy;else y=y-dy;end x=x;F=abs(xe*y)-abs(ye*x); endif (abs(F)<=0.01)F=0;endxx(index

2、)=x;yy(index)=y; index=index+1;if (abs(x)+abs(y)>=(abs(xe)+abs(ye) break ;endend if (xe>0)xxx=0:0.01:xe; elsexxx=0:-0.001:xe; end yyy=K*xxx;if (xe>0)plot(xx,yy,'g>' ,xxx,yyy,'b' );elseplot(xx,yy,'g-<' ,xxx,yyy,'b' );end grid on;圖表1當(dāng) dx=0.5;dy=0.6;xe=-5;

3、ye=-6;時，結(jié)果如圖表 1 所示。當(dāng) dx=0.6;dy=0.4; xe=6;ye=4;原因是matlab在倒數(shù)第二步計(jì)算果不為0，仍需繼續(xù)走一步?？梢詫Τ绦蜃鲞M(jìn)一步修改：加入條件if (abs(x)+abs(y)>=(abs(xe)+abs(ye)if (abs(abs(x)+abs(y)-abs(xe)-abs(ye)<dmin)走rh.步。if圖表2時結(jié)果如圖標(biāo)2所示，結(jié)果顯示多走了一步，(abs(x)+abs(y)>=(abs(xe)+abs(ye)時，結(jié)并把dmi n=mi n(dx,dy)/2;修改成結(jié)果如圖3所示，這樣就不會多圖表3采用不等步長和簡化的偏差計(jì)算公式，結(jié)果如圖表 4所示；我們發(fā)現(xiàn)對于，簡化的 偏差計(jì)算不適用于不等步長的情況。圖表4在土表所采用的程序基礎(chǔ)上對偏差計(jì)算公式做一個修改，F(xiàn)=F-

溫馨提示

  • 1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間，僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論