【在线求助】符号型变量转数值型变量问题(matlab 7.8.0(R2009a))！

nizp1982

『我的本意是求解和角度为alpha1相同的q15函数值的其它角度』
程序如下：
a=sym('a');
alpha1=20;
q15=1.0/(sin(alpha1*pi/180)^4*(33*cos(alpha1*pi/180)^4-30*cos(alpha1*pi/180)^2+5))
func15=1/(sin(a*pi/180.0)^4*(33*cos(a*pi/180.0)^4-30*cos(a*pi/180.0)^2+5))-q15
S=solve(func15,a)                    %到这里求解结果含有z和k，错误！
S=vpa(S,6)
S1=double(S)
[m,n]=size(S1)

运行提示：
alpha1 =
    20
q15 =
   17.2342
func15 =
1/(sin((pi*a)/180)^4*(33*cos((pi*a)/180)^4 - 30*cos((pi*a)/180)^2 + 5)) - 2425491604553925/140737488355328

S =
-(180*(pi*k + (i*log(z))/2))/pi
  
S =
- 180.0*k - 28.6479*i*log(z)

??? Error using ==> sym.double at 29
DOUBLE cannot convert the input expression into a double array.
If the input expression contains a symbolic variable, use the VPA function
instead.

Error in ==> main at 16
      S1=double(S)
调了很久，还是不行，希望大家给出出注意，谢谢！

swinder

不是很明白，要干什么，不过猜测尝试运行了哈，不知是否达到你本来的要求。
alpha1=20;
x=solve(1/(sin(a*pi/180.0)^4*(33*cos(a*pi/180.0)^4-30*cos(a*pi/180.0)^2+5))...
    -1.0/(sin(alpha1*pi/180)^4*(33*cos(alpha1*pi/180)^4-30*cos(alpha1*pi/180)^2+5)),a)

nizp1982

你的也是错的！
x =

-(180*(pi*k + (i*log(z))/2))/pi
此答案是错误的！！！！！！
我的本意是求解和角度为alpha1相同的q15函数值的其它角度

nizp1982

-(180*(pi*k + (i*log(z))/2))/pi
这个结果是错误的！结果至少包含20在内，还有其它几个角度的

messenger

可能是你的Matlab版本太低吧，高版本算出来的是正确的，-160 160 -20 20 ................

nizp1982

你的意思是我的程序没问题？   我的版本是matlab 7.8.0(R2009a)   
我也觉得没什么问题，但是就是提示出错！  
请问哪里可以下载到高精度或者更高的版本，谢谢版主！

messenger

R2009a符号计算用的是Mupad内核，其符号计算能力远远不行。你用R2008a及以下的版本算一下，应该行，不过版本也别太老了。

nizp1982

I see,3ku，谢谢版主

tjqhnu

messenger 发表于 2009-11-24 20:17
R2009a符号计算用的是Mupad内核，其符号计算能力远远不行。你用R2008a及以下的版本算一下，应该行，不过版 ...

R2009a中如何定义符号变量进行运算？用sym和syms进行定义不行。请版主进行解答，谢谢