[問題] 解斜率已回收
題目如下
圖:http://ppt.cc/4i2i
數學式:http://ppt.cc/M;a~
L1=85.*(1-es1.^2);
c1=85.*es1;
a2=30./(1.0000000001-es2.^2);
c2=a2.*es2;
% 解聯立方程組
eq1='-L1/(2*c1)*x1+L1=y1';
eq2='((x1-(2*c1+c2))/(a2))^2+y1^2/(a2^2-c2^2)=0';
[x1 y1]=solve(eq1,eq2);
for m=1:1001
X1(1:2,m)=subs(x1,c1(:,m));
Y1(1:2,m)=subs(subs(y1,a2(:,m)),c2(:,m));
end
% [x1 y2] 其中x1 y2皆各有兩組解
% 給定L1
% 解Solv 表示為 (L1 es1 c d slope)
% C點座標(2*(c1+c2) , 0)
% slope為 第四象限解與C點之斜率
% slov = [L1 es1 x1 y1 slope]
Solv(:,1)=L1;
Solv(:,2)=es1;
Solv(:,3)=X1(2,:);
Solv(:,4)=Y1(2,:);
Solv(:,5)=(0-Y1(2,:))./(2*(c1+c2)-X1(2,:));
最後解出來只有L1和es1
數學上應該沒錯...只是不知道指令上是否有問題?
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.117.32.5
※ 編輯: sunproxima 來自: 140.117.32.5 (10/26 21:20)