如何求解一个方程的实根

tianssss

已知t1=pi/2;t2=pi/2;
      a= 0.4787;
      b= 0.1005;


要求解的方程为
Q=4*L^2*((L^2-1)*sin(t1)^2*sin(t2)^2+sin(L*(t1-t2))^2-L^2*sin(t1-t2)^2)*a^2+4*L^2*(sin(t1)^2*(sin(L*t2)^2-L^2*sin(t2)^2)+sin(t2)^2*( sin(L*t1)^2-L^2*sin(t1)^2))*a*b+4*(sin(L*t1)^2-L^2*sin(t1)^2)*(sin(L*t2)^2-L^2*sin(t2)^2)*b^2+2*(sin(L*t1)^2-sin(L*t2)^2-L^2*sin(t1)^2+L^2*sin(t2)^2-2*L^2*(sin(t2)^2*sin(L*t1)^2-sin(t1)^2*sin(L*t2)^2))*a+4*L^2*(sin(t2)^2*sin(L*t1)^2-sin(t1)^2*sin(L*t2)^2)*b+ sin(L*(t1+t2))^2-L^2*sin(t1+t2);

这个方程看似复杂，其实a，b，t1,t2已经知道,只是L的函数。大侠帮忙看看，能不能画出Q关于L在 -1到0  区间的曲线。或者，解出L。

mhkmars

1# tianssss

1. 
   
2. t1=pi/2;t2=pi/2;a=0.4787;b=0.1005;
   
3. syms L
   
4. Q=4*L^2*((L^2-1)*sin(t1)^2*sin(t2)^2+sin(L*(t1-t2))^2-L^2*sin(t1-t2)^2)*a^2+4*L^2*(sin(t1)^2*(sin(L*t2)^2-L^2*sin(t2)^2)+sin(t2)^2*( sin(L*t1)^2-L^2*sin(t1)^2))*a*b+4*(sin(L*t1)^2-L^2*sin(t1)^2)*(sin(L*t2)^2-L^2*sin(t2)^2)*b^2+2*(sin(L*t1)^2-sin(L*t2)^2-L^2*sin(t1)^2+L^2*sin(t2)^2-2*L^2*(sin(t2)^2*sin(L*t1)^2-sin(t1)^2*sin(L*t2)^2))*a+4*L^2*(sin(t2)^2*sin(L*t1)^2-sin(t1)^2*sin(L*t2)^2)*b+ sin(L*(t1+t2))^2-L^2*sin(t1+t2);
   
5. Q=char(Q);
   
6. Q=inline(Q);
   
7. fplot(Q,[-1,0])
   

复制代码

图如下；

qibbxxt

2# mhkmars
不错，不过还是建议你用匿名函数比较好，内联函数是7以前的用法，为了顾及以前的用户，Matlab仍然保留这个用法，但是从效率上参数的共享角度来说，还是匿名函数更具优势些

1. 
   
2. t1=pi/2;t2=pi/2;a=0.4787;b=0.1005;
   
3. Q=@(L)4*L^2*((L^2-1)*sin(t1)^2*sin(t2)^2+sin(L*(t1-t2))^2-L^2*sin(t1-t2)^2)*a^2+4*L^2*(sin(t1)^2*(sin(L*t2)^2-L^2*sin(t2)^2)+sin(t2)^2*( sin(L*t1)^2-L^2*sin(t1)^2))*a*b+4*(sin(L*t1)^2-L^2*sin(t1)^2)*(sin(L*t2)^2-L^2*sin(t2)^2)*b^2+2*(sin(L*t1)^2-sin(L*t2)^2-L^2*sin(t1)^2+L^2*sin(t2)^2-2*L^2*(sin(t2)^2*sin(L*t1)^2-sin(t1)^2*sin(L*t2)^2))*a+4*L^2*(sin(t2)^2*sin(L*t1)^2-sin(t1)^2*sin(L*t2)^2)*b+ sin(L*(t1+t2))^2-L^2*sin(t1+t2);
   
4. fplot(Q,[-1,0])
   

复制代码

mhkmars

3# qibbxxt
恩，谢谢斑竹的提醒。最近也发现用匿名函数可以使代码紧凑很多，参数传递也很方便，准备系统的学习一下匿名函数的用法，\(^o^)/~

tianssss

多谢mhkmars。谢谢您的帮助。我是使用点乘的方法解决了这个问题。你们提供的方法更好。多谢多谢。也算结贴了。

allenfieldin

:hug: 5# tianssss
跟版主学匿名函数，哈哈，以前觉得这东西很神奇:hug:

sxgywjj2009

:victory:学习一下，看来高手很多的。