这段代码有很多问题,我认为它会更好,如果你现在在数学符号函数D。
多元函数的泰勒展开
42视图(30天)
显示旧的评论
大家好,
你能帮我找到这个非线性多元函数的一阶泰勒展开式的3 * M变量。我尝试matlab函数泰勒关于多变量函数,但我得到了错误,“没有足够的输入参数”。
这个函数是matlab代码:
w = [1 0 0 0];
M =长度(w);
x =符号(“x”[1]);
y =符号(“y”[1]);
z =符号(“z”[1]);
p0 = 0(3米);
a_f = @(θ,φ)0;
为m = 1: m
f = @(θ,φ,x, y, z) (exp (1 j *ρ* (x (m)。* sin(θ)。* cos(φ)+…
y (m)。* sin(θ)。*罪(φ)+ z (m)。* cos(θ))))。* w (m);
a_f = @(θ,φ,x, y, z) a_f(θ,φ)+ f(θ,φ);
结束
%定义被积函数
被积函数= @(θ,φ,x, y, z) abs (a_f(θ,φ))。^ 2。* sin(θ);
%定义集成的限制
theta_limits =(0,π);
phi_limits =(0, 2 *π);
%计算积分使用Matlab的integral2函数
Denominator_phi = @(φ,x, y, z) int(@(θ)被积函数(θ,φ,x, y, z), theta_limits (1) theta_limits (2));
分母= @ (x, y, z) int(@(φ)Denominator_phi(φ,x, y, z), phi_limits (1) phi_limits (2));
%计算最大的数组的因素
分子= @ (x, y, z) abs (a_f (theta_max phi_max)) ^ 2;
%计算方向性
D = @ (x, y, z) 4 *π*分子/分母(x, y, z)。(x, y, z);
L_D =泰勒(D,…
[x (1) (2) x (3) (4) (1) y y y (2) (3) (4) (1) z z z (2) (3) (4)],…
重塑(p0 ', 1, []),“秩序”,1)
5个评论
Torsten
约3小时前
编辑:Torsten
约2小时前
信谊θφρtheta_max phi_max真实
w = [1 0 0 0];
M =长度(w);
x =符号(“x”[1]);
y =符号(“y”[1]);
z =符号(“z”[1]);
假设(x,“真实”的)
假设(y,“真实”的)
假设(z,“真实”的)
p0 = 0(3米);
房颤(θ,φ)= ((exp (1 j *ρ*和(x (1: M)。* sin(θ)。* cos(φ)+…
y (1: M)。* sin(θ)。*罪(φ)+ z (1: M)。* cos(θ))))。* w (1: M))
被积函数= AF(θ,φ)* AF(θ,φ)。*罪(θ)
%定义集成的限制
theta_limits =(0,π);
phi_limits =(0, 2 *π);
分母= int (int(被积函数,θ,theta_limits (1) theta_limits(2)),φ,phi_limits (1) phi_limits (2))
分子= 4 *π* AF (theta_max phi_max) * AF (theta_max phi_max) '
D =分子/分母
