作业帮 > 综合 > 作业

如何在matlab下求解三角函数方程

来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/05/12 04:25:48
如何在matlab下求解三角函数方程
matlab 三角函数方程求解
10 - 离问题结束还有 13 天 13 小时
碰到一个问题,想在matlab下求解这样一个方程:
tan(k*(1.57E-6))/(k*(1.57E-6))=1/0.2159,
我用了solve命令,结果出现错误,
我的用法是:
syms k
k=solve('tan(k*(1.57E-6))/(k*(1.57E-6))=1/0.2159')
不知道这样的方程应该用什么命令?
%%下面方法比较笨,但是可行.
%%先看一下为什么你解不出来
%%把你的k*1.57e-6看成一个整体x
%%变成求解tan(x)./x-1/0.2159==0
x=-2:0.0001:2;
y=tan(x)./x-1/0.2159;
plot(x,y,x,zeros(size(x)))
title('[-2,2]')
figure
x=-1.5:0.0001:1.5;
y=tan(x)./x-1/0.2159;
plot(x,y,x,zeros(size(x)))
title('[-1.5,1.5]')
%%从上面两张图可以看出
%%第一幅图中由于存在-1.57(-pi/2)时,tan变成无穷大
%%这里精确解我也不知道怎么取到就是了.
%%不过放大第二幅图(放大那个交点),解应当在-1.420到-1.419之间
figure
x=-1.420:0.0000001:-1.419;
y=tan(x)./x-1/0.2159;
plot(x,y,x,zeros(size(x)))
for i=1:length(x)-1
if (y(i)*y(i+1)