[問題] 哪裡錯了?已回收
題目如下
http://ppt.cc/M;a~
clear all
clc
%範圍給定
for m=1:1001
es1=(m-0.999999999999)/1000; % 表示es1=0:0.001:1
es2=(m-0.999999999999)/1000; % 表示es2=0:0.001:1
%條件與參數
L1=85.*(1-es1.^2);
c1=85.*es1;
a2=30./(1.0000000001-es2.^2);
c2=a2.*es2;
% 解聯立方程組
% L1必須大於27.5
% 當L1>27.5 另其解(x1,y1)=(0,0) 方便看
eq1='-L1/(2*c1)*x1+L1';
eq2='((x1-(2*c1+c2))/(a2))^2+y1^2/(a2^2-c2^2)=0';
[x1 y1]=solve(eq1,eq2);
X1=subs(x1,c1);
Y1=subs( subs(y1,a2),c2);
if L1>27.5
xx1=[0 0];
yy1=[0 0];
else
xx1=X1;
yy1=Y1;
end
% [x1 y2] 其中x1 y2皆各有兩組解
% 第一組解(x1,y1)=(a,b)=( xx1(1),yy1(1) ) 第二象限的解
% 第二組解(x1,y1)=(c,d)=( xx1(2),yy1(2) ) 第四象限的解
% 給定L1
% 解Solv 表示為 (L1 es1 c d slope)
% C點座標(c1+c2 , 0)
% slope為 第四象限解與C點之斜率
% slov = [L1 es1 x1 y1 slope ]
Solv(m,1)=L1;
Solv(m,2)=es1;
Solv(m,3)=xx1(2);
Solv(m,4)=yy1(2);
slope=-(0-yy1(2))/((c1+c2)-xx1(2));
Solv(m,5)=slope;
end
數據跑出來都是L1<27.5才有值,L1>27.5的斜率全部都是零
和題目設定的限制條件好像矛盾了
請問我哪裡錯了 還是題目有問題?
懇請各位幫忙
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.117.32.5
※ 編輯: sunproxima 來自: 140.117.32.5 (10/18 20:01)
→
10/18 21:29, , 1F
10/18 21:29, 1F