- 积分
- 0
- 注册时间
- 2006-11-13
- 仿真币
-
- 最后登录
- 1970-1-1
|
本帖最后由 messenger 于 2009-11-24 20:36 编辑
『我的本意是求解和角度为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)
调了很久,还是不行,希望大家给出出注意,谢谢! |
|