function zhuhanshu
clear all;clc
format short
global e
L0=[0.5 0.5];
e1=linspace(0.07,0.1,10);
for j=1:length(e1)
    e=e1(j);
    L=fsolve(@zihanshu,L0);
    L1(j)=L(1);
    L2(j)=L(2);
end
[e1' L1' L2']


function f=zihanshu(x)
global e
A=-0.30714;
B=0.12966;
alpha=0.0206;
L1=x(1);
L2=x(2);
f(1)=A*alpha*sin(alpha*(L2-L1))-B*alpha*cos(alpha*(L2-L1))+e;
f(2)=A*cos(alpha*(L2-L1))-B*sin(alpha*(L2-L1));

12楼: Originally posted by waterque at 2014-06-13 18:23:44
设x = L2 - L1,
f(x) = Acos(ax) + Bsin(ax) % a ---> alpha
则f'(x) = -Aasin(ax) + Bacos(ax)
那么你的问题就变成：
f(x) = 0;
f'(x) = epsilon_m;

并不能得到具体的L1, L2.

13楼: Originally posted by waterque at 2014-06-13 18:54:15
matlab代码：
clear;
A = -0.30714;
B = 0.12966;
a = 0.0206;
epm = linspace(0.07,0.1,10);
x = -1500:0.01:1500;
f = A*cos(a*x)+B*sin(a*x);
h = -B*a*cos(a*x)+A*a*sin(a*x);
plot(x,f)
hold on;
...

11楼: Originally posted by 月只蓝 at 2014-06-13 14:42:24
function zhuhanshu
clear all;clc
format short
global e
L0=;
e1=linspace(0.07,0.1,10);
for j=1:length(e1)
    e=e1(j);
    L=fsolve(@zihanshu,L0);
    L1(j)=L(1);
    L2(j)=L(2);
end



...

10楼: Originally posted by dingd at 2014-06-13 09:34:14
方程、系数等都对吗？似乎没有实数解。

14楼: Originally posted by Littlecaesar at 2014-06-14 00:08:57
非常感谢大神的回复，一个方程正是另外一个求导所得。刚开始我以为是解方程有问题所以给出了其中一组特殊值。但是我的方程是对应的一个具体问题。我也想画一个你画的那个图，但是里面有个参数B是跟待求的一个未知数 ...

17楼: Originally posted by waterque at 2014-06-14 13:39:10
这其实就是解一个二元的方程的组，有第二个式子解出B的表达式，然后代入第一个式子中。然后得到 x = L2 - L1, 在利用B的定义解L1.

matlab codes:

function [] = functionsolver( G )


t_max = 0.52E-6;
...

17楼: Originally posted by waterque at 2014-06-14 13:39:10
这其实就是解一个二元的方程的组，有第二个式子解出B的表达式，然后代入第一个式子中。然后得到 x = L2 - L1, 在利用B的定义解L1.

matlab codes:

function [] = functionsolver( G )


t_max = 0.52E-6;
...

19楼: Originally posted by Littlecaesar at 2014-06-15 10:06:42
大神我还想请教下，这种函数不可以用fsolve直接解吗？...