An error occurred while propagating data type 'double' through...

127 次查看（过去 30 天）

Joe Jones 2022-7-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1761370-an-error-occurred-while-propagating-data-type-double-through

评论： Joe Jones 2022-7-26

Dear

I'm doing an iterative algorithm in MATLAB/Simulink/Matalb Function. However, when i run the simulation, the errors occurred.(please see the attached picture and code). If possible, please help me. Thank you very much for your help.

function gamma = cg(gamma0,voutx,tao,e1,e2,M)

beta1=gamma0(1);

beta2=gamma0(2);

beta3=gamma0(3);

beta4=gamma0(4);

beta5=gamma0(5);%gamma0 is a vector.

xd=20;%设定位置x

T=0.01;%积分步长

falf=fal(tao,0.5,0.1);

fal1=dfal(e1,0.7,0.1);

fal2=dfal(e2,0.95,0.1);

b0=1;

k=1;

N=100000;%times maximum

epsilon=1e-5;%parameter

gamma=gamma0;

while(k<N)

g=gradient(xd,voutx,fal1,fal2,falf,T,b0,M,gamma);%gradient

beta1=beta1-2*(xd-voutx)*M*(beta4*fal1*T*tao)/(1+(beta4*fal1*T+beta5*fal2)*T*b0);

beta2=beta2-2*(xd-voutx)*M*T*falf*(beta4*fal1*T+beta5*fal2)/(1+(beta4*fal1*T+beta5*fal2)*T*b0);

beta3=beta3-2*(xd-voutx)*M*((beta4*fal1*T+beta5*fal2)*T*T+T/b0)*falf/(1+(beta4*fal1*T+beta5*fal2)*T*b0);

beta4=beta4-2*(xd-voutx)*M*fal1/(1+(beta4*fal1*T+beta5*fal2)*T*b0);

beta5=beta5-2*(xd-voutx)*M*fal2/(1+(beta4*fal1*T+beta5*fal2)*T*b0);%the three parameters that need to be adjusted

gamma=[beta1;beta2;beta3;beta4;beta5];%vector

if(norm(g)<epsilon), break; end

end

19 个评论
显示 17更早的评论隐藏 17更早的评论

Joe Jones 2022-7-19

编辑：Joe Jones 2022-7-19

Dear

I modify the code according to your comments.I click "edit data" and change the 'DataType' to 'double'. Then the errors are solved. However, a new error appears:

"Output argument 'gamma' is not assigned on some execution paths."

My 'gamma' is the output of the matlab function, but my output is not in the if-else statement. If possible, please help me. Thank you. @Fangjun Jiang

function gamma=frcg(gamma0,voutx,tao,e1,e2,M)

k=0;

n=length(gamma0);

N = 10000;

epsilon = 1e-5;

xd=20;

T=0.01;

b0=1;

falf=fal(tao,0.5,0.1);

fal1=dfal(e1,0.7,0.1);

fal2=dfal(e2,0.95,0.1);

g0=zeros(5,1);

d0=zeros(5,1);

while(k<N)

g=tidu(xd,voutx,fal1,fal2,falf,T,b0,M,gamma0,tao);

itern=k-(n+1)*floor(k/(n+1));

itern=itern+1;

if(itern==1)

d=-g;

else

beta=(g'*g)/(g0'*g0);

d=-g+beta*d0; gd=g'*d;

if(gd>=0.0)

d=-g;

end

if(norm(g)<epsilon), break; end

alpha=-(g'*d)/(d'*d);

gamma0=gamma0+alpha*d;

g0=g; d0=d;

k=k+1;

gamma=gamma0;

end

Joe Jones 2022-7-20

Dear

I have added "gamma=gamma0" as the first line of the function, however, the function does't work. The error is the same as former. The change has been mentioned with '%' in the attached code. If possible, please help me. Thank you!@Fangjun Jiang

function gamma=frcg(gamma0,voutx,tao,e1,e2,M)

k=0;

n=length(gamma0);

N = 10000;

epsilon = 1e-5;

xd=20;

T=0.01;

b0=1;

falf=fal(tao,0.5,0.1);

fal1=dfal(e1,0.7,0.1);

fal2=dfal(e2,0.95,0.1);

g0=zeros(5,1);

d0=zeros(5,1);

g=g0;

while (k<N)

gamma=gamma0;%I have added this statement to the first line of 'while'.

g=tidu(xd,voutx,fal1,fal2,falf,T,b0,M,gamma0,tao);

itern=k-(n+1)*floor(k/(n+1));

itern=itern+1;

if(itern==1)

d=-g;

else

beta=(g'*g)/(g0'*g0);

d=-g+beta*d0; gd=g'*d;

if(gd>=0.0)

d=-g;

end

if norm(g)>epsilon, break; end

alpha=-(g'*d)/(d'*d);

gamma0=gamma0+alpha*d;

g0=g; d0=d;

k=k+1;

%gamma=gamma0; This statement has been moved to the firat line.

end

Joe Jones 2022-7-21

Dear

I am sorry for misunderstanding your comments.

I modified the code according to your comments. The previous error was indeed resolved, but there is a new one.(please see the attached figure)

I have an ideal: Is it because it may exponentiate negative numbers that it produces complex value?(e.g. a^b: a<0, b is fraction)

I have put the 'tidu' functionn, 'fal' function and 'dfal' function behind main function. If possible, please help me. Thank you ! @Fangjun Jiang @Walter Roberson

function gamma=frcg(gamma0,voutx,tao,e1,e2,M)

gamma=gamma0;%I have added this statement to the first line.

k=0;

n=length(gamma0);

N = 10000;

epsilon = 1e-5;

xd=20;

T=0.01;

b0=1;

falf=fal(tao,0.5,0.1);

fal1=dfal(e1,0.7,0.1);

fal2=dfal(e2,0.95,0.1);

g0=zeros(5,1);

d0=zeros(5,1);

while (k<N)

g=tidu(xd,voutx,fal1,fal2,falf,T,b0,M,gamma0,tao);

itern=k-(n+1)*floor(k/(n+1));

itern=itern+1;

if(itern==1)

d=-g;

else

beta=(g'*g)/(g0'*g0);

d=-g+beta*d0; gd=g'*d;

if(gd>=0.0)

d=-g;

end

if(norm(g)<epsilon), break; end

alpha=-(g'*d)/(d'*d);

gamma0=gamma0+alpha*d;

g0=g; d0=d;

k=k+1;

gamma=gamma0;

end

%tidu function

function g=tidu(xd,voutx,fal1,fal2,falf,T,b0,M,gammak,tao)

g1k=(gammak(4)*fal1*T*tao)/(1+(gammak(4)*fal1*T+gammak(5)*fal2)*T*b0);

g2k=T*falf*(gammak(4)*fal1*T+gammak(5)*fal2)/(1+(gammak(4)*fal1*T+gammak(5)*fal2)*T*b0);

g3k=((gammak(4)*fal1*T+gammak(5)*fal2)*T*T+T/b0)*falf/(1+(gammak(4)*fal1*T+gammak(5)*fal2)*T*b0);

g4k=fal1/(1+(gammak(4)*fal1*T+gammak(5)*fal2)*T*b0);

g5k=fal2/(1+(gammak(4)*fal1*T+gammak(5)*fal2)*T*b0);

g = -2*(xd-voutx)*M*[g1k;g2k;g3k;g4k;g5k];

%dfal functionn

function y=dfal(x,alpha,delta)

if x>delta

y=alpha*x^(alpha-1);

else if x<-delta

y=(-1)^(alpha+1)*alpha*x^(alpha+1);

else

y=delta^(alpha-1);

end

%fal function

function y=fal(x,alpha,delta)

if abs(x)>delta

y=abs(x)^alpha*sign(x);

else

y=x/(delta^(1-alpha));

end

Fangjun Jiang 2022-7-25

This is probably where the help from this community ends. The simulation of the Simulink model runs but it is extremely slow and it stopped due to instability of the simulation. You need to consult the expert in your field to resolve the problem. Most likely the parameters of the ADRC is too aggressive. Avoid directly derivative.

Joe Jones 2022-7-26

Dear

Thank you for your time and effort you have put into my problems. Your help is of great significance for me to solve the problem. I will continue to do my best to resolve these errors. Thank you again!@Fangjun Jiang

Wish you every success!

请先登录，再进行评论。

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