erfc with complex number

18 次查看（过去 30 天）

显示更早的评论

rahman 2011-5-10

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/7153-erfc-with-complex-number

评论： Dilip Jose 2020-2-13

采纳的回答： Matt Fig

Hi all

I want to compute erfc(.) of a complex number. How can I do that? Is there any way to use erfc(.) of matlab?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

采纳的回答

Matt Fig 2011-5-10

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/7153-erfc-with-complex-number#answer_9795

A quick search of the FEX turned up these two. There may be others...

File 1

File 2

14 个评论
显示 12更早的评论隐藏 12更早的评论

Dilip Jose 2020-1-30

在 MATLAB Online 中打开

close all
syms z;
gr=5;
gc=5;
M=1;
w=0.5;
pr=7;
sc=2.01;
t=0.4;
m=(M+(2*1i*w));
s1=z;
s2=(2*sqrt(t));
x=((s1)./(s2));
a1=(x*t);
a2=(2*m*(sqrt(pi)));
a=((m)./(pr-1));
b=((m)./(sc-1));
x=((s1)./(s2));
u1=((x^2+(m*t))*t);
u2=4*m;
d1=(x*(sqrt(t)))*(1-(4*m*t));
d2=(8*m^3/2);
v1=gr;
v2=(a*(1-pr));
k1=gc;
k2=(b*(1-sc));
h1=exp(a*t);
h2=2;
e1=(exp(2*x*(sqrt(m*t))));
e2=(erfc(x+(sqrt(m*t))));
e3=(exp(-2*x*(sqrt(m*t))));
e4=(erfc(x-(sqrt(m*t))));
b1=(exp(2*x*((sqrt((m+a)*t)))));
b2=(erfc(x+(sqrt((m+a)*t))));
b3=(exp(-2*x*(sqrt((m+a)*t))));
b4=(erfc(x-(sqrt((m+a)*t))));
b5=(exp(2*x*((sqrt((m+b)*t)))));
b6=(erfc(x+(sqrt((m+b)*t))));
b7=(exp(-2*x*(sqrt((m+b)*t))));
b8=(erfc(x-(sqrt((m+b)*t))));
e5=exp(-(x^2+(m*t)));
c1=(exp(2*x*(sqrt(pr*a*t))));
c2=(erfc(x*(sqrt(pr)+(sqrt(a*t)))));
c3=((exp(-2*x*(sqrt(pr*a*t)))));
c4=(erfc(x*(sqrt(pr)-(sqrt(a*t)))));
l1=(exp(2*x*(sqrt(sc*b*t))));
l2=(erfc(x*(sqrt(sc)+(sqrt(b*t)))));
l3=((exp(-2*x*(sqrt(sc*b*t)))));
l4=(erfc(x*(sqrt(sc)-(sqrt(b*t)))));
f1=(exp(b*t));
f2=2;
l=((a1)./(a2));
f=((f1)./(f2));
h=((h1)./(h2));
j1=((u1)./(u2));
j2=((d1)./(d2));
j3=((v1)./(v2));
j4=((k1)./(k2));
q1=(2*(j1*((e1*e2)+(e3*e4))+(j2*((e3*e4)-(e1*e2))-(l*e5))));
q2=((j3+j4)*(1/2)*((e1*e2)+(e3*e4)));
q3=((j3*(h*(b1*b2)+(b3*b4))));
q4=((j4*(f*(b5*b6)+(b7*b8))));
q5=(j3*(erfc(x*sqrt(pr))));
q6=(j4*erfc(x*sqrt(sc)));
q7=(j3*(h*((c1*c2)+(c3*c4))));
q8=(j4*(f*((l1*l2)+(l3*l4))));
q=(q1+q2-q3-q4-q5-q6+q7+q8);
fplot(z,q);
xlim([0 5])
ylim([0 1])
hold on
title('Axial velocity profiles for different values of gr & gc');
hold off

Walter Roberson 2020-2-5

在 MATLAB Online 中打开

To check, you have

d2=(8*m^3/2);

Is it possible that that should be

d2=(8*m^(3/2));

The current expression would be equivalent to the clearer

d2=4*m^3;

Or even 8/2 instead of 4 if that somehow made documentation sense. With the /2 where it is and with the leading multiplier already even, people are going to wonder.

Dilip Jose 2020-2-13