请高手帮解个方程式，急，谢谢！

beetle521

arccos[(l*sinb-1)/(l+1)]-arccos[(l*sinb-1)/(l-1)]-c=0
l为所求
b、c输入数据

[ 本帖最后由 beetle521 于 2009-1-5 11:22 编辑 ]

ch_j1985

看看此帖
http://forum.vibunion.com/forum/ ... mp;page=1#pid394719

不爱陌生人

b=input('b=');
c=input('c=');
l=solve('acos[(l*sinb-1)/(l+1)]-acos[(l*sinb-1)/(l-1)]-c=0')

beetle521

我算出的结果是2.168459……
而人家给的结果是4.53635

请问这是怎么回事？
谢谢

beetle521

我算出的结果是2.168459……
而人家给的结果是4.53635

请问这是怎么回事？
谢谢

ch_j1985

b、c贴出来大家试试

beetle521

b=40°
c=12.5581°

ChaChing

我试过LZ的, 怎解不出!? 版本问题吗? 我的 v6.5

ch_j1985

将度数转化为弧度

beetle521

十分感谢！

beetle521

一个新问题，我想做一个界面：

K，α已知
θ=180*(K-1)/(K+1)
arccos[(t*sinα-1)/(t+1)]-arccos[(t*sinα-1)/(t-1)]-θ=0求出t.


一个push button为计算，我编写程序为


global K theta
theta=180*(K-1)/(K+1);
theta=num2str(theta);
set(handles.theta,'string',theta);


global phi t theta
theta=get(handles.theta,'string');
theta=str2num(theta);
t=solve('acos((t*sin(phi*pi/180)-1)/(t+1))-acos((t*sin(phi*pi/180)-1)/(t-1))-theta*pi/180');
t=num2str(t);
set(handles.t,'string',t);

运行后，提示内存不足，请问应该如何解决，谢谢！

ch_j1985

对界面不熟悉
上面的两段代码放在同一个push button的CallBack里面吗？

beetle521

是的，就是计算这个push button下

beetle521

想问一下关于解方程函数的问题

例如：
temp2=40*pi/180;
temp3=12.5581*pi/180;
solve('acos((t*sin(temp2)-1)/(t+1))-acos((t*sin(temp2)-1)/(t-1))-temp3',0)

感觉我的问题就是这个solve函数，有没有什么函数允许式子中有参数呢？

beetle521

程序修改后为：
temp2=phi/180*pi;
temp3=theta/180*pi;
[email=f=@(t)acos((t*sin(temp2)-1)/(t+1))-acos((t*sin(temp2)-1)/(t-1))-temp3]f=@(t)acos((t*sin(temp2)-1)/(t+1))-acos((t*sin(temp2)-1)/(t-1))-temp3[/email];
t=fzero(f,[4,5])

因为用solve函数时容易丢根，所以求解前可以先作图，然后指定解出某区间的的一个解

找了一个作图求解的例子，可以参考：
下面的程序作出了[0，10]上的图像，通过观察，指定解出了[2，3]内的一个解，注意fzero函数每次只能找到一个解

X=[0:.1:10];Y=zeros(1,length(X));
plot(X,1.5*sin(X)-sin(pi/3+X/3),X,Y,'r');
grid on
axis equal

f=@(X)1.5*sin(X)-sin(pi/3+X/3);
z=fzero(f,[2,3])