青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

平凡的世界

神鷹忽展翅,頭頂青天飛
隨筆 - 10, 文章 - 0, 評(píng)論 - 34, 引用 - 0
數(shù)據(jù)加載中……

三維管線建模Matlab仿真程序(只做到插值)

進(jìn)程被槍斃了.心里不痛快.把做了一半的matlab仿真程序拿出來(lái)供大家拍拍磚吧.寫(xiě)得比較亂.有不明白的可以回復(fù)或是發(fā)郵件給我.我的郵件是
phoenix8848@gmail.com

 1clear;
 2A=[0 0 0;0.65 0.25 0.376;1 0 0];R=0.25*1.414;
 3%calculate line vectors
 4v=ones(1,3);%initial line vectors
 5for i=1:2
 6    v(i,:)=A(i+1,:)-A(i,:);
 7end
 8%calculate the interpolation angle
 9ang=acos(abs(v(1,1)*v(2,1)+v(1,2)*v(2,2)+v(1,3)*v(2,3))/(sqrt(v(1,1)^2+v(1,2)^2+v(1,3)^2)*sqrt(v(2,1)^2+v(2,2)^2+v(2,3)^2)));
10%the dis pB to pAi
11d=R*tan(ang/2);
12%direction cosine of V(Ai-1,Ai)
13cosa1=v(1,1)/sqrt(v(1,1)^2+v(1,2)^2+v(1,3)^2);
14cosb1=v(1,2)/sqrt(v(1,1)^2+v(1,2)^2+v(1,3)^2);
15cosc1=v(1,3)/sqrt(v(1,1)^2+v(1,2)^2+v(1,3)^2);
16%interpolation start point
17b1=[A(2,1)-d*cosa1,A(2,2)-d*cosb1,A(2,3)-d*cosc1];
18%direction cosine of V(Ai,Ai+1)
19cosa2=v(2,1)/sqrt(v(2,1)^2+v(2,2)^2+v(2,3)^2);
20cosb2=v(2,2)/sqrt(v(2,1)^2+v(2,2)^2+v(2,3)^2);
21cosc2=v(2,3)/sqrt(v(2,1)^2+v(2,2)^2+v(2,3)^2);
22%interpolation end point
23b2=[A(2,1)+d*cosa2,A(2,2)+d*cosb2,A(2,3)+d*cosc2];
24%vector of angular bisector
25v0=sqrt(v(1,1)^2+v(1,2)^2+v(1,3)^2);
26v1=[v(1,1) v(1,2) v(1,3)]/v0;
27v0=sqrt(v(2,1)^2+v(2,2)^2+v(2,3)^2);
28v2=[v(2,1) v(2,2) v(2,3)]/v0;
29L0=[(-v1(1)+v2(1))/2,(-v1(2)+v2(2))/2,(-v1(3)+v2(3))/2];
30%direction cosine of angular bisector
31cosa3=L0(1)/sqrt(L0(1)^2+L0(2)^2+L0(3)^2);
32cosb3=L0(2)/sqrt(L0(1)^2+L0(2)^2+L0(3)^2);
33cosc3=L0(3)/sqrt(L0(1)^2+L0(2)^2+L0(3)^2);
34%calculate the center point of interpolation circle
35t=R*sec(ang/2);
36p=[A(2,1)+cosa3*t,A(2,2)+cosb3*t,A(2,3)+cosc3*t];
37%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
38%calculate the normal vector of interpolation circle
39dx=[b1(1)-p(1) b1(2)-p(2) b1(3)-p(3)];
40dx=dx/sqrt(dx(1)^2+dx(2)^2+dx(3)^2);
41dy=[b2(1)-p(1) b2(2)-p(2) b2(3)-p(3)];
42dy=dy/sqrt(dy(1)^2+dy(2)^2+dy(3)^2);
43dz=[v(1,2)*v(2,3)-v(2,2)*v(1,3),v(1,1)*v(2,3)-v(2,1)*v(1,3),v(1,1)*v(2,2)-v(2,1)*v(1,2)];
44dz=dz/sqrt(dz(1)^2+dz(2)^2+dz(3)^2);
45
46T=[dx;dy;dz];
47%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
48P=ones(3,11);
49delta=pi*0.5/10;
50for i=0:10
51    m=[R*cos(i*delta) R*sin(i*delta) 0];
52    P(:,i+1)=T'*m'+p';
53end
54plot3(P(1,:),P(2,:),P(3,:),'r+-',A(:,1),A(:,2),A(:,3),'b+-');
55grid on;
56axis on;
57xlabel('X');ylabel('Y');zlabel('z');


附上以上程序在matlab6.5環(huán)境下的運(yùn)行結(jié)果:


posted on 2007-10-31 10:03 西門有悔 閱讀(1823) 評(píng)論(2)  編輯 收藏 引用

評(píng)論

# re: 三維管線建模Matlab仿真程序(只做到插值)  回復(fù)  更多評(píng)論   

02 年的數(shù)學(xué)建模~ 比較sure~
2007-11-02 12:11 | jerry

# re: 三維管線建模Matlab仿真程序(只做到插值)  回復(fù)  更多評(píng)論   

@jerry

我也不清楚是不是02年的數(shù)學(xué)建模題.

我不是純學(xué)數(shù)學(xué)的,感覺(jué)這個(gè)問(wèn)題從思路上不是特別復(fù)雜.

關(guān)鍵是一些數(shù)學(xué)知識(shí)的運(yùn)用,比如四元數(shù)\空間坐標(biāo)框架的仿射變換等等.

然后就是一些線性代數(shù)方程組的求解與矩陣的運(yùn)算
2007-11-02 20:57 | 西門有悔

只有注冊(cè)用戶登錄后才能發(fā)表評(píng)論。
網(wǎng)站導(dǎo)航: 博客園   IT新聞   BlogJava   博問(wèn)   Chat2DB   管理


青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            红桃视频国产一区| 免费亚洲电影| 国产亚洲美州欧州综合国| 欧美在线亚洲| 久久青青草综合| 99国产精品久久久久久久| 亚洲日韩欧美一区二区在线| 欧美精品久久久久久久免费观看 | 性欧美精品高清| 先锋影音国产精品| 亚洲国产99精品国自产| 亚洲人体影院| 国产欧美一区二区三区视频| 欧美成人按摩| 国产精品vip| 久久综合色8888| 欧美日韩精品综合| 久久精品99国产精品日本| 麻豆久久精品| 香蕉乱码成人久久天堂爱免费 | 欧美视频中文字幕| 久久高清免费观看| 欧美成人黑人xx视频免费观看| 亚洲天堂网在线观看| 欧美在线观看网站| 亚洲午夜精品| 久久资源在线| 欧美一级在线播放| 欧美高清在线观看| 久久免费国产精品1| 欧美日韩一区二区三区在线看 | 亚洲少妇中出一区| 亚洲国产导航| 亚洲欧美日韩国产一区| 亚洲精品中文字幕在线| 欧美在线二区| 亚洲欧美在线一区二区| 欧美高清一区| 欧美mv日韩mv国产网站| 国产欧美日韩一区二区三区在线| 亚洲大胆人体视频| 禁久久精品乱码| 亚洲欧美成人网| 亚洲一区二区欧美日韩| 欧美激情综合色| 在线观看亚洲专区| 亚洲欧美日韩在线综合| 亚洲伊人久久综合| 欧美日韩国产成人在线| 亚洲国产精品t66y| 亚洲第一网站免费视频| 欧美专区一区二区三区| 久久国产视频网站| 国产精品永久免费观看| 中文av一区特黄| 亚洲图片欧洲图片日韩av| 久久精品国产免费观看| 久久国产免费看| 欧美午夜电影完整版| 亚洲美女黄网| 亚洲视频在线免费观看| 欧美日韩一区三区| 99re6热只有精品免费观看| 亚洲国产精品视频| 裸体歌舞表演一区二区| 亚洲高清毛片| 一本久道综合久久精品| 欧美日韩三级| 亚洲一区二区三区视频播放| 校园春色国产精品| 国产日韩欧美在线播放| 久久久久久9| 欧美激情片在线观看| 99爱精品视频| 国产精品久久久久久久7电影| 在线亚洲伦理| 久久久久免费视频| 亚洲高清久久网| 欧美了一区在线观看| 亚洲精品久久嫩草网站秘色| 亚洲中午字幕| 国内一区二区在线视频观看| 美女黄网久久| 亚洲另类自拍| 欧美一区二区三区视频| 伊人久久大香线蕉综合热线| 欧美成人小视频| 亚洲午夜av| 免费观看成人鲁鲁鲁鲁鲁视频| 亚洲精品一区二区在线| 欧美日韩高清在线播放| 欧美一区二区三区啪啪| 亚洲高清网站| 欧美专区在线播放| 亚洲精品欧美日韩专区| 国产精品日韩在线观看| 老鸭窝91久久精品色噜噜导演| 日韩一区二区精品| 美女诱惑黄网站一区| 国产精品99久久不卡二区| 国产日韩视频| 欧美日本韩国一区| 久久久999精品免费| 亚洲精品影视| 女人香蕉久久**毛片精品| 亚洲综合国产精品| 91久久国产综合久久| 国产欧美在线看| 欧美日韩国产不卡在线看| 久久精品一区二区国产| 亚洲午夜精品久久久久久浪潮| 欧美成人四级电影| 久久国产婷婷国产香蕉| 一本到12不卡视频在线dvd| 激情亚洲网站| 国产区精品在线观看| 欧美视频久久| 欧美成人国产| 欧美大尺度在线观看| 欧美专区日韩视频| 亚洲自拍电影| 99这里只有久久精品视频| 狠狠综合久久| 国产亚洲欧美一区二区| 欧美视频第二页| 欧美精品久久99| 欧美91大片| 免费成人在线观看视频| 久久人人爽国产| 久久都是精品| 欧美在线关看| 久久av一区二区| 亚洲欧美综合网| 亚洲欧美另类国产| 亚洲一区二区三区午夜| 亚洲一二区在线| 宅男噜噜噜66一区二区| 亚洲最黄网站| 一区二区三区高清在线| 999在线观看精品免费不卡网站| 亚洲人成高清| 亚洲日本免费电影| 亚洲毛片一区二区| 夜夜嗨av一区二区三区中文字幕| 最近中文字幕mv在线一区二区三区四区| 欧美激情精品久久久久久免费印度 | 久久最新视频| 欧美成人蜜桃| 欧美剧在线观看| 欧美性大战久久久久久久| 国产精品欧美日韩一区| 国产视频精品网| 一区二区亚洲| 亚洲精选久久| 午夜久久福利| 久久亚洲一区二区| 欧美电影打屁股sp| 亚洲久久视频| 午夜精品久久久久久久99水蜜桃| 久久国产精品99久久久久久老狼| 久久女同互慰一区二区三区| 欧美高清视频在线| 国产精品男女猛烈高潮激情| 国产视频在线观看一区二区三区| 国外视频精品毛片| 一区电影在线观看| 久久成人资源| 亚洲电影激情视频网站| 亚洲五月六月| 免费欧美在线| 国产欧美日本一区二区三区| 伊人天天综合| 亚洲尤物在线视频观看| 你懂的国产精品| 一本大道久久a久久精品综合| 久久精品2019中文字幕| 欧美.www| 一本一道久久综合狠狠老精东影业| 欧美主播一区二区三区| 欧美精选一区| 激情久久综合| 亚洲在线视频网站| 亚洲电影观看| 性欧美大战久久久久久久久| 欧美另类一区| 怡红院精品视频在线观看极品| 一区二区三区蜜桃网| 蜜桃av综合| 午夜精品久久久久久久白皮肤| 欧美~级网站不卡| 国内精品伊人久久久久av影院| 国产精品99久久久久久久久| 免费欧美电影| 欧美亚洲日本网站| 欧美三级电影一区| 亚洲精品视频一区| 国产精品网站在线播放| 国产精品麻豆成人av电影艾秋| 国产精品v亚洲精品v日韩精品 | 欧美日韩精品在线视频|