Question

求教：用matlab解下面的非线性参联立方程

方程1：(log((alpha - 1)*(c + omega + k*r - a*f(k)))*((alpha - 1)*(c + omega + k*r - a*f(k)))^((Y*alpha)/(Y*alpha - X*(alpha - 1)))*(Y/(Y*alpha - X*(alpha - 1)) + (Y*alpha*(X - Y))/(Y*alpha - X*(alpha - 1))^2) + (Y*alpha*((alpha - 1)*(c + omega + k*r - a*f(k)))^((Y*alpha)/(Y*alpha - X*(alpha - 1)) - 1)*(c + omega + k*r - a*f(k)))/(Y*alpha - X*(alpha - 1)))/(c + omega - alpha*(c + omega + k*r - a*f(k)))^((X*(alpha - 1))/(Y*alpha - X*(alpha - 1))) - ((log(c + omega - alpha*(c + omega + k*r - a*f(k)))*(X/(Y*alpha - X*(alpha - 1)) + (X*(X - Y)*(alpha - 1))/(Y*alpha - X*(alpha - 1))^2))/(c + omega - alpha*(c + omega + k*r - a*f(k)))^((X*(alpha - 1))/(Y*alpha - X*(alpha - 1))) - (X*(alpha - 1)*(c + omega + k*r - a*f(k)))/((Y*alpha - X*(alpha - 1))*(c + omega - alpha*(c + omega + k*r - a*f(k)))^((X*(alpha - 1))/(Y*alpha - X*(alpha - 1)) + 1)))*((alpha - 1)*(c + omega + k*r - a*f(k)))^((Y*alpha)/(Y*alpha - X*(alpha - 1)))=0

方程2：(X*((alpha - 1)*(c + omega + k*r - a*f(k)))^((Y*alpha)/(Y*alpha - X*(alpha - 1)))*(alpha - 1)^2)/((Y*alpha - X*(alpha - 1))*(c + omega - alpha*(c + omega + k*r - a*f(k)))^((X*(alpha - 1))/(Y*alpha - X*(alpha - 1)) + 1)) + (Y*alpha*((alpha - 1)*(c + omega + k*r - a*f(k)))^((Y*alpha)/(Y*alpha - X*(alpha - 1)) - 1)*(alpha - 1))/((Y*alpha - X*(alpha - 1))*(c + omega - alpha*(c + omega + k*r - a*f(k)))^((X*(alpha - 1))/(Y*alpha - X*(alpha - 1))))=0

要联立上述两个方程解出alpha和c用其他参数表示的值，应该怎么做？？求高手指点！！万分感谢！！！！！

Answer 1

matlab的solve函数具有解联立参数方程的功能，调用语法参见下例：

[x,y]=solve('x^2+y^2=a','x+c*y=b','x,y')

不过看你这两个方程的复杂程度，想求出解析解估计很难，建议先尽可能化简或者求数值解