matlab求解倒数第三个解,matlab求解数值解的问题

求大神帮忙分析一下，求方程组的数值解的问题，程序如下，目的就是让L(w)=|TAO|^2-1;G(w)=αL(w)/w

求解方程组L(w)=0,G(w)=0时的变量q1,q2的值。希望有大神帮忙，感激万分。

clc;clear all;close all;

syms q1 q2 w real;

% q1=0.8;

% q2=0.8;

%--------------------------%参数初始化-------------------------------%

a=0.8;

b=0.7;

r=1;

I=eye(2);

f=1.57;

A1=[0 -1; 0 0];

A2=[0 1; 0 0];

A3=[0 -1;a*f -a-b];

A4=[0 1;0 b];

B1=[0 1];

B11=(1/r)*B1*(B1');

PI11=[sqrt(q1*q2+2*sqrt(q1^3*r)) -sqrt(q1*r);-sqrt(q1*r) sqrt(q2*r+2*sqrt(q1*r^3))];

a1=sqrt(q1*r);

b1=-sqrt(q2/r+2*sqrt(q1/r));

vecPI11=PI11(:);

m=-inv((kron(I,(A1'-PI11*B11))+kron(A3',I)))*(kron(A4',I));

m0=-inv((kron(I,(A1'-PI11*B11))+kron(A3',I)))*(kron(A2',I));

vecPI12=m0*vecPI11; %PI12

PI12=reshape(vecPI12,2,2);

a2=-1/r*PI12(2,1);

b2=-1/r*PI12(2,2);

vecPI13=m*m0*vecPI11; %PI13

PI13=reshape(vecPI13,2,2);

a3=-1/r*PI13(2,1);

b3=-1/r*PI13(2,2);

vecPI14=m^2*m0*vecPI11; %PI14

PI14=reshape(vecPI14,2,2);

a4=-1/r*PI14(2,1);

b4=-1/r*PI14(2,2);

vecPI15=m^3*m0*vecPI11; %PI15

PI15=reshape(vecPI15,2,2);

a5=-1/r*PI15(2,1);

b5=-1/r*PI15(2,2);

s=1i*w;

G1=s^2-b1*s+a1;

T0=(b*s+a*f)/(s^2+(a+b)*s+a*f);

F6=(a5+b*b5)*s+PI15(1,1)+b*PI15(1,2);

%-----------------------%transfer function------------------------------------%

TAO=(a1*T0^4+(s*F6)/G1+(a2+(b2*s-a2)*T0)*T0^3+(a3+(b3*s-a3)*T0)*T0^2+(a4+(b4*s-a4)*T0)*T0+(a5+(b5*s-a5)*T0))/G1;  %transfer function

Q1=abs(TAO)-1;%TAO的模

Q2=diff(Q1,w);%TAO的模一阶倒数

w=0.5

fsolve(@(q) [Q1 Q2],[1 1])%求解参数 q1,q2