double(limit) trouble

조회 수: 4 (최근 30일)
Maxim Prokhorov
Maxim Prokhorov 2016년 6월 14일
답변: Debarati Banerjee 2016년 6월 23일
Good day. I have tried execute this code:
clear; clc;
ln = @(x) log(x);
scale = 1e-10; hbar = 1.055*1e-34; me = 9.11*1e-31; eV2J = 1.6*1e-19; J2eV = eV2J^(-1);
ma = 0.067; mb = 0.15;
b = 17; a = 45; c = 100;
dz = 1;
eV = 0;
deV = eV/(b+a+b);
rangeU = [
0
1 - deV*b/2
0 - deV*(b + a/2)
1 - deV*(b + a + b/2)
0 - eV
];
rangeM = [ ma, mb, ma, mb, ma ];
U( 1 : c ) = rangeU(1); m( 1 : c ) = rangeM(1)*me;
U( c+1 : c+b ) = rangeU(2); m( c+1 : c+b ) = rangeM(2)*me;
U( c+b+1 : c+b+a ) = rangeU(3); m( c+b+1 : c+b+a ) = rangeM(3)*me;
U( c+b+a+1 : c+b+a+b ) = rangeU(4); m( c+b+a+1 : c+b+a+b ) = rangeM(4)*me;
U( c+b+a+b+1 : c+b+a+b+c ) = rangeU(5); m( c+b+a+b+1 : c+b+a+b+c ) = rangeM(5)*me;
L = c + b + a + b + c;
z = 0:dz:L;
syms Zi Ui mi Ez t r;
ki = (2*mi*(Ez-Ui)*eV2J).^(0.5)./hbar;
Mi = [
exp( 1i*ki*Zi ), exp(-1i*ki*Zi )
ki./mi.*exp( 1i*ki*Zi ), -ki./mi.*exp(-1i*ki*Zi )
];
invMi = inv(Mi);
productM = [1 0; 0 1];
for index = length(z)-1:-1:2
localProduct = subs(invMi, {mi, Ui, Zi}, {m(index), U(index), z(index)*scale})*subs(Mi, {mi, Ui, Zi}, {m(index-1), U(index-1), z(index)*scale});
productM = productM*localProduct;
end
outWave = [t, 0]'; inWave = [1, r]';
[t, r] = solve( outWave == productM*inWave, t, r);
% outAmp = double(subs(t, Ez, 0)); Division by zero.
outAmp = limit(t, Ez, 0)
outAmp = double(outAmp);
And I have got:
outAmp =
limit((((exp(-(conj(((2315911741694958125423557058355*Ez)/118571099379011784113736688648896417641748464297615937576404566024103044751294464)^(1/2))*683982766270...... Output truncated. Text exceeds maximum line length of 25,000 characters for Command Window display.
Error using symengine
DOUBLE cannot convert the input expression into a double array.
Error in sym/double (line 616)
Xstr = mupadmex('symobj::double', S.s, 0);
Error in test (line 54)
outAmp = double(outAmp);
What can I do to get double value of outAmp?

답변 (1개)

Debarati Banerjee
Debarati Banerjee 2016년 6월 23일
When the input argument to DOUBLE function is a symbolic expression, MATLAB will output the error message you received.
You can workaround the issue by substituting suitable values to the variable Ez and then applying 'double' to the result.
>> t1 = subs(t,Ez,.001)
>> double(t1)
Cheers,
Debarati

카테고리

Help CenterFile Exchange에서 Conversion Between Symbolic and Numeric에 대해 자세히 알아보기

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by