Identify an unrecognized function or variable

Question

Mathew 2025-3-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2174926-identify-an-unrecognized-function-or-variable

评论： Mathew 2025-3-14

clear all; close all;
syms T U t w
a=0.3; M=1; h0=1; h1=1; k0=1; k1=1; p=0.5; d=0.5; Z=10; 
%t=0:1:10; %w=0.75;
H=@(s) (p*Z-d-exp(-s))*k1; 
K=@(r,g) p*Z*k1-(d+exp(-r))*k1-2*exp(a*r).*g;
for n=1:10
y=int(H(U),0,10);
yy= h0 + ((1-a)./M)*H(t).*h1 + (a./M)*y;
h1=yy;
end
for n=1:10
z=int(K(U,T),0,10);
zz= k0 + ((1-w)./M).*K(t,yy).*k1 + (w./M)*z;
k1=zz;
end
zzfcn = matlabFunction(zz)  
zzfcn = function_handle with value:
    @(T,t,w)w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0)+1.0)+1.0)+1.0)+1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0)+1.0)+1.0)+1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0)+1.0)+1.0)+1.0
t=0:0.1:10;                                                             
w = 0.1:0.1:0.75;                                                     
[T,W] = ndgrid(t,w);                                                   
figure
surf(T, W, zzfcn(T,W))                                                  
Not enough input arguments.

Error in sym/matlabFunction>@(T,t,w)w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0)+1.0)+1.0)+1.0)+1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0)+1.0)+1.0)+1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0)+1.0)+1.0)+1.0
xlabel('t')
ylabel('\omega')
zlabel('k(t)')
set(gca, 'ZScale','log')                                               
colormap(turbo)
colorbar

2 个评论
显示无隐藏无

Torsten 2025-3-8

If you mean your function "zzfcn" instead of "fcn": it has input arguments T,t and w. Thus you have to call it with three instead of two variables as you do in the surf-command.

Mathew 2025-3-8

@Torsten Thanks

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Sulaymon Eshkabilov 2025-3-8

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2174926-identify-an-unrecognized-function-or-variable#answer_1561439

在 MATLAB Online 中打开

Note that the function zzfcn has three input variables: @(T,t,w)

Here is the corrected code:

clear all; close all;

syms T U t w

a=0.3; M=1; h0=1; h1=1; k0=1; k1=1; p=0.5; d=0.5; Z=10;

%t=0:1:10; %w=0.75;

H=@(s) (p*Z-d-exp(-s))*k1;

K=@(r,g) p*Z*k1-(d+exp(-r))*k1-2*exp(a*r).*g;

for n=1:10

y=int(H(U),0,10);

yy= h0 + ((1-a)./M)*H(t).*h1 + (a./M)*y;

h1=yy;

end

for n=1:10

z=int(K(U,T),0,10);

zz= k0 + ((1-w)./M).*K(t,yy).*k1 + (w./M)*z;

k1=zz;

end

zzfcn = matlabFunction(zz)

zzfcn = function_handle with value:

@(T,t,w)w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0).*(w.*(exp(-1.0e+1)-T.*(exp(3.0).*(1.0e+1./3.0)-1.0e+1./3.0).*2.0+4.4e+1)+(w-1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0)+1.0)+1.0)+1.0)+1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0)+1.0)+1.0)+1.0).*(exp(-t)+exp(t.*(3.0./1.0e+1)).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-1.0e+1).*(3.0./1.0e+1)-(exp(-t).*(7.0./1.0e+1)-6.3e+1./2.0e+1).*(exp(-t).*(-7.0./1.0e+1)+exp(-1.0e+1).*(3.0./1.0e+1)+3.47e+2./2.0e+1)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0)+7.1e+1./5.0).*2.0-9.0./2.0)+1.0)+1.0)+1.0

t=0:0.1:10;

w = 0.1:0.1:0.75;

[T,W] = ndgrid(t,w);

figure

surf(T, W, zzfcn(T, T,W))

xlabel('t')

ylabel('\omega')

zlabel('k(t)')

set(gca, 'ZScale','log')

colormap(turbo)

colorbar

●

11 个评论
显示 9更早的评论隐藏 9更早的评论

Walter Roberson 2025-3-10

在 MATLAB Online 中打开

for i=1:length(t)
    y=int(H(U),0,10);
    yy= h0 + ((1-w)./M).*H(t).*h1 + (w./M).*y;
    h1=yy;
end

t is a vector.

yy involves H(t) where t is the entire t vector. So yy returns a vector the same size as the row vector t.

for i=1:length(t)
    z=int(K(yy,t),0,10);
    zz= k0 + ((1-w)./M).*K(yy,t).*k1 + (w./M).*z;
    k1=zz;
end

yy is a vector, and t is a vector, and these are the same size. K(yy,t) is working with vectors.

K=@(g,s) p*Z*k1-(d+exp(-s))*k1-2*exp(a.*s).*g;

Working with two vectors the same size is okay in K, with a vector the same size being returned.

So the int(K(yy,t),0,10) is a vector. The loop is not changing yy or t during the iterations, so z will be the same each iteration.

Then zz involves computing with the vector K(yy,t) together with the vector z, so zz is going to end up being a vector the same size as t.

Except...

K involves computation with the vector a and the vector s and those vectors are not the same size, so you would generate an error in the computation exp(a.*s). Unless, that is, you deliberately made a the same size as t by using something such as

a = linspace(0.1, 0.95, numel(t));

Still, the loops over length(t) that use the entire vector t are suspicious. It would make more sense if you were to be doing something such as

for i=1:length(t)
    y=int(H(U),0,10);
    yy= h0 + ((1-w)./M).*H(t(i)).*h1 + (w./M).*y;
    h1=yy;
end

and

for i=1:length(t)
    z=int(K(yy,t(i)),0,10);
    zz= k0 + ((1-w)./M).*K(yy,t(i)).*k1 + (w./M).*z;
    k1=zz;
end

Torsten 2025-3-10

编辑：Torsten 2025-3-10

在 MATLAB Online 中打开

a=linspace(0.1,0.95,10);

?

It's important that t and a are both row vectors of equal size (here:10). Otherwise your code will not work.

Mathew 2025-3-14

@Walter Robers and @Torsten Thanks

请先登录，再进行评论。

Identify an unrecognized function or variable

2 个评论
显示无隐藏无

回答（1 个）

11 个评论
显示 9更早的评论隐藏 9更早的评论

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

Identify an unrecognized function or variable

2 个评论 显示 无隐藏 无

回答（1 个）

11 个评论 显示 9更早的评论隐藏 9更早的评论

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

2 个评论
显示无隐藏无

11 个评论
显示 9更早的评论隐藏 9更早的评论