Question

我试图写一个计算每个步骤的固定组成部分价值的文件。职能F非常复杂。简言之,它包含了从tn-(n/2)到s=tn+(n/2)的一组术语。在融合之后,F仍然有变数。因此,你可以说F(t) = 整体(f(s,t)),在“F = @t......”栏中,f(s,t)是“int(int)”之后的主要术语。在最后一行,我想对F(t)在F计算 f(s,t)的整体性之后的tn进行评价。

然而,在进行这项工作后,我发现“不确定功能或变数”错误。

function [ bigTheta_n ] = Untitled( bigTheta_o, bigOmega_o )
nt=5001; %since (50-0)/.01 = 5000
dt = .01; % =H
H=.01;
l=.05;

bigTheta_n = ones(nt,1);
bigTheta_n(1)=bigTheta_o; %theta_o

bigOmega_n = ones(nt,1);
bigOmega_n(1)=bigOmega_o; %omega_o
littleOmega_n = ones(nt,1);

epsilon=10^(-6);
eta = epsilon*10;

t_o=0;

    function Keta = K(t)
        Keta = (422.11/eta)*exp((5*(4*((t-tn)^2)/eta^2)-1)^(-1))        
    end

    function F = F(t)
        F = int(Keta(eta*t-s)*littleOmega_n,s,tn-(n/2),tn+(n/2))
    end

for n=1:4999
    tn=t_o+n*dt;
    littleOmega_n(n) = bigOmega_n(n) - sin(bigTheta_n(n))*cos(2*pi*tn/epsilon)/(2*pi*l);
    F = @(t) int(Keta(eta*t-s)*littleOmega_n,s,tn-(n/2),tn+(n/2));
    bigTheta_n(n+1) = bigTheta_n(n) + H*F(tn)
end

end

即使我删除了“F = @(t.)”的字句,并以下文取代(并放在Keta = K(t)的定义下)

function F = F(t, ti, tf)
    F = int(Keta(eta*t-s)*littleOmega_n,s,ti,tf)
end

然后取代最后一行

bigTheta_n(n+1) = bigTheta_n(n) + H*F(tn)

iii

bigTheta_n(n+1) = bigTheta_n(n) + H*F(tn, tn-(n/2), tn+(n/2))

我仍然有同样的错误:“不确定的职能或变数”

Answer 1

引言

F = @(t) int(@(s) Keta(eta*t-s)*littleOmega_n,s,tn-(n/2),tn+(n/2));

这使综合成为匿名功能,输入s。

友情链接