苏州工业园区发布,上海seo优化,线上商城app,如何推广新产品的方法使用了Matlab代码#xff0c;化简平时遇到的计算问题#xff0c;也可以用于验算结果来自211工科专业2学分复变函数与积分变换课程求复角主值sym(angle(待求复数))%公式
sym(angle(1sqrt(3)*i))%举例代入化简将 代入关于z的函数f(z)中并化解#xff0c;用于公式法计算无穷远点…使用了Matlab代码化简平时遇到的计算问题也可以用于验算结果来自211工科专业2学分复变函数与积分变换课程求复角主值sym(angle(待求复数))%公式
sym(angle(1sqrt(3)*i))%举例代入化简将 代入关于z的函数f(z)中并化解用于公式法计算无穷远点的留数syms f z
f1/(z*(z1));%某一函数f(z)
fsubs(f,z,1/z);%可以把1/z换成任意关系式
fsimplify(f/(z^2));一维积分计算简单的积分计算matlab中不加即可显示符号型变量使用double()函数可以转化为数值型显示。syms x z(x) %定义关于x的函数z
u(x)3*(i*x^2-x^2);%某函数
jf1int(u,x,[0,i]);%这里可以是定积分、不定积分、变上限积分只需要在[]内更改
jf1%符号型显示
double(jf1)%数值型显示级数求和用于计算级数是否收敛、条件收敛、绝对收敛。对目标级数s(n)判断是否收敛直接使用symsum()级数求和函数如果需要判断是否绝对收敛时使用abs()函数取模即可。syms s n
s(n)(8^n)/factorial(n);%factorial()阶乘函数
s(n)abs(((-1)^n)/n);
qhsymsum(s,n,1,inf)
double(qh)二次积分连续计算两次积分的基础计算详细见「概率论」部分低薪小明「概率论」基础计算代码zhuanlan.zhihu.comsyms x y z(x,y)
u(x,y)6-x-y;
jf1int(u,y,[2,3-x]);
jf2int(jf1,x,[0,1]);
%jfint(z,x,[0,inf]);
jf1
jf2
%double(jf)一阶、二阶微分计算经典题型由u(x,y)构造解析函数f(z)首先对u相应的求导然后由于解析函数的性质对v积分最终化简得到f(z).clc
syms x y u
u(x,y)x^2-y^2x*y;
uxdiff(u,x)%u对x求导
uydiff(u,y)%u对y求导
vxdiff(u,x,2)
vydiff(u,y,2)
vjfyint(ux,y)
vxdiff(vjfy,x)普通求导并代入计算结合高阶求导公式计算结果先求导再将相应的z0代入求值非常简便。syms f z
f1/sin(z);
dfdiff(f,z)
subs(df,z,0)普通极限显而易见这用来求一个函数的极限用于判断极限是否存在判断是否收敛判断奇点类型求留数等等syms f n
f(n)(11/n)^n;
limit(f,n,inf)%这里也可以计算在某一点的极限不一定是无穷∞三角极限用于含三角函数的留数计算不常用syms f z
for k0:1:10f(z)(z-k*pi)/sin(z);limit(f,z,k*pi)
end对圆积分在圆上的积分用于第三章、第五章的积分计算不过使用柯西积分公式、高阶求导公式或留数公式等等也能算出。syms seita
u(seita)exp(seita.*j);%圆的指数是方程
f(seita)exp(-1*j*seita);%被积函数
Iint(f*diff(u,seita),seita,0,pi);保形映射公式计算针对第六章知道几对映射点求保形映射后的点只是简单的化简计算在代码中输入初始点就行。clc
syms z w(z);
w(z2)*i/(z-2);
w1subs(w,z,0)
scatter(real(w1),imag(w1));
hold on
w2subs(w,z,1)
scatter(real(w2),imag(w2));
hold on
w3subs(w,z,-1)
scatter(real(w3),imag(w3));
hold on
w4subs(w,z,-i)
scatter(real(w4),imag(w4));
hold on
%zuoexpand((wi)/((w3*i))*4)
%youexpand(z/(2*(z-1)))保形映射画图针对第六章对复平面某一区域求保形映射后的区域在第一行的x和y分别输入21个表示该区域边界的离散点并在循环中输入映射函数即可得到映射后图形。不过这里有一个限制条件初始区域只能是方形区域对圆形区域或者三角区域无法处理欢迎大家提出宝贵的改进意见。x-10:1:10;ylinspace(-pi,0,21);wzeros(1,441);
for j1:20for k1:21w((j-1)*21k)exp(x(j)y(k)*i);end
end
realwreal(w);
imagwimag(w);
scatter(realw,imagw);
