matlab 急初始条件为θ=0,x=-d(θ范围在0—90°之间吧)求大哥大姐们帮忙啊,本人急~用matlab求解出x和θ的关系,也可以是一系列数值解,需要程序!没有贴图吗?我已经上传啦
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/29 01:18:43
![matlab 急初始条件为θ=0,x=-d(θ范围在0—90°之间吧)求大哥大姐们帮忙啊,本人急~用matlab求解出x和θ的关系,也可以是一系列数值解,需要程序!没有贴图吗?我已经上传啦](/uploads/image/z/8558280-0-0.jpg?t=matlab+%E6%80%A5%E5%88%9D%E5%A7%8B%E6%9D%A1%E4%BB%B6%E4%B8%BA%CE%B8%3D0%2Cx%3D-d%EF%BC%88%CE%B8%E8%8C%83%E5%9B%B4%E5%9C%A80%E2%80%9490%C2%B0%E4%B9%8B%E9%97%B4%E5%90%A7%EF%BC%89%E6%B1%82%E5%A4%A7%E5%93%A5%E5%A4%A7%E5%A7%90%E4%BB%AC%E5%B8%AE%E5%BF%99%E5%95%8A%2C%E6%9C%AC%E4%BA%BA%E6%80%A5%7E%E7%94%A8matlab%E6%B1%82%E8%A7%A3%E5%87%BAx%E5%92%8C%CE%B8%E7%9A%84%E5%85%B3%E7%B3%BB%2C%E4%B9%9F%E5%8F%AF%E4%BB%A5%E6%98%AF%E4%B8%80%E7%B3%BB%E5%88%97%E6%95%B0%E5%80%BC%E8%A7%A3%2C%E9%9C%80%E8%A6%81%E7%A8%8B%E5%BA%8F%21%E6%B2%A1%E6%9C%89%E8%B4%B4%E5%9B%BE%E5%90%97%3F%E6%88%91%E5%B7%B2%E7%BB%8F%E4%B8%8A%E4%BC%A0%E5%95%A6)
matlab 急初始条件为θ=0,x=-d(θ范围在0—90°之间吧)求大哥大姐们帮忙啊,本人急~用matlab求解出x和θ的关系,也可以是一系列数值解,需要程序!没有贴图吗?我已经上传啦
matlab 急
初始条件为θ=0,x=-d(θ范围在0—90°之间吧)
求大哥大姐们帮忙啊,本人急~
用matlab求解出x和θ的关系,也可以是一系列数值解,需要程序!
没有贴图吗?我已经上传啦
matlab 急初始条件为θ=0,x=-d(θ范围在0—90°之间吧)求大哥大姐们帮忙啊,本人急~用matlab求解出x和θ的关系,也可以是一系列数值解,需要程序!没有贴图吗?我已经上传啦
没有方程怎么解
快给出方程
你先等等,这个东西需要迭代求解,速度比较慢.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function dx=nonlinear(t,x)
persistent count
if isempty(count)
count=0;
end
count=count+1;
fprintf('正在运行第%d次\n',count);
fprintf('运行到角度:%d\n',t);
theta1=85*pi/180;
r0=2000*tan(pi/6);
l=2000;
h=8;
d=15;
n=1.586;
%% 取使约束方程r0*sin(t)/sin(theta1)=n*l*sin(w)/sqrt(1-n^2*sin(w)^2)+(h+d+x)*tan(t)-x*tan(w)
%% 等式两边之差小于shredhold的最优w,如果没有满足精度要求的,则进一步精确划分W可能值,直到得到满足要求
%% 的W为止.
w_about=0;
calError=1;
shredhold=0.0005;
w=-pi;
percesion=0.001;
midVar=0.01;
flag=0; %精化等级
while w_about==0
while calError>shredhold
if w>=pi
break;
end
w=w+pi/180*percesion;
calError=abs(3172*sin(w)/sqrt(1-2.5154*sin(w)^2)+(23+x)*tan(t)-x*tan(w)-1159.1*sin(t));
if calError<midVar
midVar=[midVar calError];
w_about=[w_about w];
end
end
w=w_about((midVar==min(midVar)));
fprintf('截断误差:%f',min(midVar));
percesion=percesion/10;
flag=flag+1;
fprintf('flag:%d\n',flag);
end
dx=(h+d+x)/( -(n*cos(w)-cos(t))/(n*sin(w)-sin(t))-tan(t) )/cos(t)^2;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
求theta在[0 PI/4]之间的关系
tic,[T3,Y3] = ode45(@nonlinear,[0 pi/4],-15);toc
有些误差: