用finverse反函数求不出来，想用插值算，请高手指点

银杏树之舞

求1.5*t-2*sin(0.05*t)*tan(0.05*t)=0的值

liuyalong008

首先用ezplot做出曲线，再选择合适的初值获取数值结果

andm

syms t  x

eq=1.5*t-2*sin(0.05*t)*tan(0.05*t);
x=solve(eq,t)
结果
x =

matrix([[0]])

银杏树之舞

andm 发表于 2012-4-9 13:14
syms t  x

eq=1.5*t-2*sin(0.05*t)*tan(0.05*t);

除了0以外的答案

银杏树之舞

求方程的值t=  (2* (1 –cos?( 0.1*t ) ))/y， 其中，y=[1.5  -0.4  -1.5  -1.3];求t？
算法一：syms t
>> y=2*(1-cos(0.1*t))/t;
finverse(y,t)

ans =

0

算法二：syms  p；
>> g=solve('1.5-2*p*(1-cos(0.1/p))=0')
Warning: Explicit solution could not be found.
> In solve at 140

g =

[ empty sym ]
没出错，但没结果。
算法三：syms  t ;
s=dsolve('1.5*t-2*(1-cos(0.1*t))=0','t');
??? Error using ==> maple
at offset 2, `}` unexpected

Error in ==> dsolve at 147
var_set = maple([var_set ' intersect ' var_set]);

银杏树之舞

银杏树之舞 发表于 2012-4-10 16:21
求方程的值t=  (2* (1 –cos?( 0.1*t ) ))/y， 其中，y=[1.5  -0.4  -1.5  -1.3];求t？
算法一：syms t
>>  ...

大虾帮忙看看

银杏树之舞

银杏树之舞 发表于 2012-4-10 16:21
求方程的值t=  (2* (1 –cos?( 0.1*t ) ))/y， 其中，y=[1.5  -0.4  -1.5  -1.3];求t？
算法一：syms t
>>  ...

三种算法都不行，请大虾指点

银杏树之舞

liuyalong008 发表于 2012-4-13 15:22
首先用ezplot做出曲线，再选择合适的初值获取数值结果

呵呵，迭代太麻烦了，我已经用taylor公式简化公式了

liuyalong008

银杏树之舞 发表于 2012-4-14 11:55
呵呵，迭代太麻烦了，我已经用taylor公式简化公式了

泰勒展开  你在哪个数值附近展开的？？   那样只会有一个值
我用EZPLOT测试过了，没有问题，基本都可以解出来

银杏树之舞

liuyalong008 发表于 2012-4-17 16:34
泰勒展开  你在哪个数值附近展开的？？   那样只会有一个值
我用EZPLOT测试过了，没有问题，基本都可以解 ...

我试出来了，呵呵，谢谢:)

银杏树之舞

liuyalong008 发表于 2012-4-17 16:34
泰勒展开  你在哪个数值附近展开的？？   那样只会有一个值
我用EZPLOT测试过了，没有问题，基本都可以解 ...

还有个问题请教一下：问题一：
z1=[0.00003 0.00033 0.00034];
for z=0:0.000001:0.00034
if z>=0&z<z1(1,1)
f(z)=s(1,1);
else if z>=z1(1,1)&z<z1(1,2)
f(z) = s(1,2);
else
f(z) =s(1,3);
end
end
plot(z,f(z));

这个程序没有运行结果，why？
问题二：

f=@(x)0*(x<=-1)+(2*x+1)*(x>-1&x<1)+1*(x>=1);
fplot(f,[-2 2])
为何要加@(x)？