作业帮matlab迭代问题－－50分
来源：学生作业帮 编辑：作业帮 分类：数学作业 时间：2024/05/18 04:08:50
matlab迭代问题－－50分
用循环迭代求取c点的轨迹,和移动的g点轨迹
假设xA为已知,先取xC1＝xA；
k＝[sqrt（R0^2+(e^2-1)*xC1^2）]/xC1；此为c点的法线斜率,
xC2＝xA-r/sqrt（1＋k^2）；
zC2＝[R0－sqrt（R0^2+(e^2-1)xC2^2）] /(e^2-1)；
zA=zC2+ k(Xa- xC2)；
xG=xA＋q；
zG= zA+L；
已知
R0=8000;e^2=1.1;r=2.5;q=125;l=300;i=1;
xA为什么为1啊？补充条件：abs(b)>=10^(-8)
用循环迭代求取c点的轨迹,和移动的g点轨迹
假设xA为已知,先取xC1＝xA；
k＝[sqrt（R0^2+(e^2-1)*xC1^2）]/xC1；此为c点的法线斜率,
xC2＝xA-r/sqrt（1＋k^2）；
zC2＝[R0－sqrt（R0^2+(e^2-1)xC2^2）] /(e^2-1)；
zA=zC2+ k(Xa- xC2)；
xG=xA＋q；
zG= zA+L；
已知
R0=8000;e^2=1.1;r=2.5;q=125;l=300;i=1;
xA为什么为1啊？补充条件：abs(b)>=10^(-8)
看得不是很懂.一下程序只能作为例子参考!
clc;clear;
R0=8000;e2=1.1;r=2.5;q=125;l=300;
xA=1;
xC(1)=xA;
for i=1:100;
k(i)=sqrt(R0^2+(e2-1)*xC(i)^2)./xC(i);
xC(i+1)=xA-r/sqrt(1+k(i)^2);
zC(i+1)=(R0-sqrt(R0^2+(e2-1)*xC(i+1)^2))/(e2-1);
zA(i)=zC(i+1)+k(i)*(xA-xC(i+1));
xG(i)=xA+q;
zG(i)=zA(i)+l;
end
figure(1)
plot(xC,zC,'-')
figure(2)
plot(xG,zG,'+')
clc;clear;
R0=8000;e2=1.1;r=2.5;q=125;l=300;
xA=1;
xC(1)=xA;
for i=1:100;
k(i)=sqrt(R0^2+(e2-1)*xC(i)^2)./xC(i);
xC(i+1)=xA-r/sqrt(1+k(i)^2);
zC(i+1)=(R0-sqrt(R0^2+(e2-1)*xC(i+1)^2))/(e2-1);
zA(i)=zC(i+1)+k(i)*(xA-xC(i+1));
xG(i)=xA+q;
zG(i)=zA(i)+l;
end
figure(1)
plot(xC,zC,'-')
figure(2)
plot(xG,zG,'+')