切換到寬版
  • 廣告投放
  • 稿件投遞
  • 繁體中文
    • 3256閱讀
    • 0回復

    [求助]請教:求一根光線與平面交點問題 [復制鏈接]

    上一主題 下一主題
    離線bss6
     
    發(fā)帖
    5
    光幣
    194
    光券
    0
    只看樓主 倒序閱讀 樓主  發(fā)表于: 2010-10-29
    我想請教各位大俠,我在matlab 中使用maple('solve({ y=-2e-6,x+5e-6>=0,1e-6-x>=0,z-4.4e-5>=0,5.2e-5-z>=0,e1,e2,e3},{t,x,y,z})')語句時所遇到的難題。 IW~q,X+`V  
    我原本的目的是想解六個方程(其中有好幾個不等式方程),所以我搜尋了好久找到了maple語句,但是其中的三個方程 p10i_<J]=  
    x=xt+reflect(1)*t ]K^#'[  
    y=yt+reflect(2)*t f:!b0j  
    z=zt+reflect(3)*t B=,j$uH  
    其中xt,yt,zt,reflect(1),reflect(2),reflect(3)是嵌套在循環(huán)中,每次循環(huán)得到的數(shù)值不一樣,所以這三個方程就變化了,我設想是否可以將e1,e2,e3代表這三個方程,問題出現(xiàn)了,e1,e2,e3怎么也賦不進去。希望大俠們幫幫我,我將不勝感激!我已經(jīng)困擾了好長時間! {*lRI  
    {~=Z%Cj2Q  
    注:我做這個目的是:想求一條直線與一個有邊界條件的三維曲面的交點,所以涉及到了六個不等式,也可以不用maple,由于我是個初學者,懂得太少,所以我不知道是否還有其它辦法,也請大俠們給個建議! =0" Zse,  
    下面是我寫的代碼(是錯的) Y`tv"v2  
    xt=0 t:N3k ;k  
    yt=0 e5HHsR6  
    zt=yt+4.6e-5 XW2{I.:in>  
    reflect= [  0   -0.3965    0.9180] 9NPOdt:@  
    syms t x y z j_#oP  
    x=xt+reflect(1)*t '(S@9%,aK1  
    y=yt+reflect(2)*t 0;6 ^fiSY;  
    z=zt+reflect(3)*t ,|Xibfw  
    maple('solve({ y=-2e-6,x+5e-6>=0,1e-6-x>=0,z-4.4e-5>=0,5.2e-5-z>=0,x,y,z},{t,x,y,z})')。
     
    分享到