matlab学习笔记009之solve函数

Matlab中solve函数主要是用来求解线性方程组的解析解或者精确解。对于得出的结果是符号变量，可以通过vpa()得出任意位数的数值解！

solve函数的语法定义主要有以下四种：

solve(eq)
solve(eq, var)
solve(eq1, eq2, …, eqn)
g = solve(eq1, eq2, …, eqn, var1, var2, …, varn)

eq代表方程，var代表的是变量。

但在我的使用中，solve却不能帮我解决问题：

我的目的：已知控制对象传函G0和期望截止频率wc，求在任意幅值处的频率

一开始是这样用的：

    wc=solve('20*log10(abs(G0（wc）)=10*log10(3)','wc'); %G0（wc）是关于wc的具体传函，10*log10(3)就是具体幅值

结果调试时显示wc是空矩阵，提示：方程（20*log10(abs(G0（wc）)=10*log10(3)）无效。后来在前面加这个：sym wc；

还是不行。

然后又使用这样的方法：

      for w=0.1:0.1:20
           L=20*log10(abs(wn^2/((w*i)*((i*w)^2+2*zeta*w*wn*i+6^2))));
             if L==(10*log10(3)-0.001)
               wc=w;
            end
      end

L就是具体的传函.调试时还是不行，wc依旧是空。发现if选择语句根本没运行。

是不是w的精度小了，无法使if的条件（L==10*log10(3)）满足？

然后我把w精度改为0.01,0.001.都不行。

接着我想是不是这样的精度都满足不了条件，于是乎我把条件改了：L>=(10*log10(3)-0.001)&&L<=(10*log10(3)+0.001)

最终还是不行！！！

还有一个问题：我测试了abs（）这个取模值函数和solve（）函数，发现如果是解实数方程，OK，没有问题；但若是复数方程就不行了：

      

>> clear
>> syms x y
>> [x,y]=solve('x^2+y^2=1','x','y')
 
x =
 
 -(1 - z^2)^(1/2)
  (1 - z^2)^(1/2)
 
 
y =
 
 z
 z

>> clear
>> syms x y
>> [x,y]=solve('abs(x*i+y^2)=1','x','y')
 
x =
 
- exp(2*u*i)*i + z^2*i
 
 
y =
 
z

这是为什么呢？