发布网友
共1个回答
热心网友
简单写如下:
t0的初值理论上可随便给,意味着求t0附近的零点,
step是为了求到正的t,不用也行,见后面关于t0和step的说明
V0=10:5:100;
Alpha=0:5:90;
h0=1.8;
g=10;
function t=f(t0,V0,Alpha,h0,g,step)
N1=length(V0);
N2=length(Alpha);
N=N1*N2; %没对参数对应一个方程,一共有N种组合
t=zeros(N,1); %初始化
k=1; %数组t的序号
for i=1:N1
for j=1:N2
v0=V0(i);
alpha=Alpha(j); %每对参数组合 (V0(i),Alpha(j))
t_k=fzero(@nested_f,t0); %试解nested_f方程在t0附近的零点
while t_k