如何使用Matlab程序求解极值点（函数表达式有两个符号变...
显示 更早的评论
clear all
close all
clc
syms sita1 sita2;
Z=[1 0 0
0 exp(i*2*pi/3) 0
0 0 exp(i*4*pi/3)];
X=[0 0 1
1 0 0
0 1 0];
XZ=X*Z;
[A,V1]=eig(Z);
[B,V2]=eig(X);
[C,V3]=eig(XZ); %求解特征向量；
c=[1/sqrt(3);1/sqrt(3)*exp(i*sita1);1/sqrt(3)*exp(i*sita2)];
CC=[c,c,c];
z=sum(abs(dot(B,CC)))+sum(abs(dot(C,CC)));
%做内积，此后，要对z进行求解极值点，z为sita1,sita2的函数。
%----------------------------------------------------------------------
dzsita1=diff(z,sita1)
dzsita2=diff(z,sita2)
S=solve(dzsita1==1,dzsita2==1,sita1,sita2); %S求解不出来，程序在此中遇到的难题了
采纳的回答
sihexa
2022-11-25
0 个投票
你的表达式有点复杂，用 solve 无法求出符号解。
可以考虑用 fsolve 求解数值解，不过，需要将方程组定义为函数句柄。你可以将最后一句改为：
f = matlabFunction([dzsita1;dzsita1]);
f = @(x) f(x(1),x(2));
fsolve(f,[1 1])
最后一句的 [1 1] 是方程组求解的初值，你要是能知道解的大致范围，就换做相应的范围。数值求解的一个缺点就是对初值敏感，必须指定一个合理的初值
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 编程 的更多信息
