银杏树之舞 发表于 2012-4-9 09:52:50

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

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

liuyalong008 发表于 2012-4-9 09:52:51

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

andm 发表于 2012-4-9 13:14:47

syms tx

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

matrix([])

银杏树之舞 发表于 2012-4-9 14:44:03

andm 发表于 2012-4-9 13:14 static/image/common/back.gif
syms tx

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


除了0以外的答案

银杏树之舞 发表于 2012-4-10 16:21:36

求方程的值t=(2* (1 –cos?( 0.1*t ) ))/y, 其中,y=;求t?
算法一:syms t
>> y=2*(1-cos(0.1*t))/t;
finverse(y,t)

ans =

0

算法二:symsp;
>> 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 ]
没出错,但没结果。
算法三:symst ;
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();

银杏树之舞 发表于 2012-4-10 16:23:14

银杏树之舞 发表于 2012-4-10 16:21 static/image/common/back.gif
求方程的值t=(2* (1 –cos?( 0.1*t ) ))/y, 其中,y=;求t?
算法一:syms t
>>...

大虾帮忙看看

银杏树之舞 发表于 2012-4-10 16:36:03

银杏树之舞 发表于 2012-4-10 16:21 static/image/common/back.gif
求方程的值t=(2* (1 –cos?( 0.1*t ) ))/y, 其中,y=;求t?
算法一:syms t
>>...

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

银杏树之舞 发表于 2012-4-14 11:55:07

liuyalong008 发表于 2012-4-13 15:22 static/image/common/back.gif
首先用ezplot做出曲线,再选择合适的初值获取数值结果

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

liuyalong008 发表于 2012-4-17 16:34:34

银杏树之舞 发表于 2012-4-14 11:55 static/image/common/back.gif
呵呵,迭代太麻烦了,我已经用taylor公式简化公式了

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

银杏树之舞 发表于 2012-4-18 09:10:58

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

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

银杏树之舞 发表于 2012-4-18 09:32:05

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

还有个问题请教一下:问题一:
z1=;
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)?

页: [1]
查看完整版本: 用finverse反函数求不出来,想用插值算,请高手指点