作业帮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)
来源：学生作业帮助网 编辑：作业帮 时间：2024/12/01 20:07:39
![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)](/uploads/image/z/8954674-34-4.jpg?t=matlab%E8%BF%AD%E4%BB%A3%E9%97%AE%E9%A2%98%EF%BC%8D%EF%BC%8D50%E5%88%86%E7%94%A8%E5%BE%AA%E7%8E%AF%E8%BF%AD%E4%BB%A3%E6%B1%82%E5%8F%96c%E7%82%B9%E7%9A%84%E8%BD%A8%E8%BF%B9%2C%E5%92%8C%E7%A7%BB%E5%8A%A8%E7%9A%84g%E7%82%B9%E8%BD%A8%E8%BF%B9%E5%81%87%E8%AE%BExA%E4%B8%BA%E5%B7%B2%E7%9F%A5%2C%E5%85%88%E5%8F%96xC1%EF%BC%9DxA%EF%BC%9Bk%EF%BC%9D%5Bsqrt%EF%BC%88R0%5E2%2B%28e%5E2-1%29%2AxC1%5E2%EF%BC%89%5D%2FxC1%EF%BC%9B%E6%AD%A4%E4%B8%BAc%E7%82%B9%E7%9A%84%E6%B3%95%E7%BA%BF%E6%96%9C%E7%8E%87%2CxC2%EF%BC%9DxA-r%2Fsqrt%EF%BC%881%EF%BC%8Bk%5E2%EF%BC%89%EF%BC%9BzC2%EF%BC%9D%5BR0%EF%BC%8Dsqrt%EF%BC%88R0%5E2%2B%28e%5E2-1%29)
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)
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)
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)
看得不是很懂.一下程序只能作为例子参考!
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+(...
全部展开
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,'+')
收起