matlab解如下方程并画图x^2+y^2=625;y0^2+z0^2=900;a=arcsin(50/56.43);tan(a)=z0/y0;(y-y0)*(z-z1)=(y-y1)*(z-z0);(z-z0)*(60-y)=(z-z1)*sqrt((y-y0)^2+z0^2);plot(x,y1,z1,'r*')
来源:学生作业帮助网 编辑:作业帮 时间:2024/05/05 20:01:50
matlab解如下方程并画图x^2+y^2=625;y0^2+z0^2=900;a=arcsin(50/56.43);tan(a)=z0/y0;(y-y0)*(z-z1)=(y-y1)*(z-z0);(z-z0)*(60-y)=(z-z1)*sqrt((y-y0)^2+z0^2);plot(x,y1,z1,'r*')
matlab解如下方程并画图
x^2+y^2=625;
y0^2+z0^2=900;
a=arcsin(50/56.43);
tan(a)=z0/y0;
(y-y0)*(z-z1)=(y-y1)*(z-z0);
(z-z0)*(60-y)=(z-z1)*sqrt((y-y0)^2+z0^2);
plot(x,y1,z1,'r*')
matlab解如下方程并画图x^2+y^2=625;y0^2+z0^2=900;a=arcsin(50/56.43);tan(a)=z0/y0;(y-y0)*(z-z1)=(y-y1)*(z-z0);(z-z0)*(60-y)=(z-z1)*sqrt((y-y0)^2+z0^2);plot(x,y1,z1,'r*')
代码供参考:
n = 1.5;% x1 = φ,x2 = Rdx = @(z,x)[1/(x(2)*(n+(n+1)*(2*cos(x(1))+cos(x(1))*cos(x(1))))); -1/tan(x(1))];[t,x]=ode45(dx,[0 1],[0.5 0.5]);plot(t,x(:,2))
你现在的问题有很多不清楚的地方:
1、第一个表达式很乱,我不确定是不是正确理解了你的意图;
2、没给初始条件;
3、常数n是多大?这个微分方程似乎属于刚性的,不同的n可能会有影响,而使用不同的ode求解器结果可能也不一样.