我试图创建一个等于卷积积分的函数 F(t),然后通过取值 F(tn) 来计算 theta_n(n+1)。但是,我收到错误“‘char’类型的输入参数的未定义函数‘F’”。有什么问题吗?
function [ theta ] = Untitled( theta_o,omega_o )
nt=5001; %since (50-0)/.01 = 5000
dt = .01; % =H
H=.01;
theta_n = ones(nt,1);
theta_n(1)=0; %theta_o
omega_n = ones(nt,1);
omega_n(1)=-0.4; %omega_o
epsilon=10^(-6);
eta = epsilon*10;
t_o=0;
for n=1:4999
tn=t_o+n*dt;
F := int((422.11/eta)*exp((5*(4*((eta*t-s-tn)^2)/eta^2)-1)^(-1))*omega, s,tn- (n/2),tn+(n/2))
theta_n(n+1) = theta_n(n) + h*F(tn);
end
end
那
F = @(tn) ...;
呢?
您可能还混淆了这段代码中的其他一些内容,例如
H
与 h
等等...