Re: [問題] 解三角函數數值問題
※ 引述《roy80110 (培特羅)》之銘言：
: 有一函數如下
: tanh(x)*tanh(x*y)+tan(x)*tan(x*y)=0
: 他解出來的圖如下
: http://images.plurk.com/ll7j-2U4e35BqWndVMtRQpWRub5.jpg
: (我將u設做x，x/L設做y)
: 但我用下面的指令去跑，跑出來的結果全為0
: syms y
: for x=[1:1:1000]
: Y(1,x) = solve(tanh(x/100)*tanh(x*y/100)+tan(x/100)*tan(x*y/100)==0)
: end
: 不知道為什會這樣QQ
: 還是說有其他方法可以解出y?
: 有請各位大大解惑，謝謝
首先,雙曲函數基本上是正數,也就是說
若x y 沒有負數的話是很難讓方程式可以成立的
我有先用圖形來看,非常醜
後面用數值解來解
-----------------------------------------------------
function pttex136
global x n
% 圖形分析(以下可以刪掉)
figure(1)
ezsurfc('tanh(x)*tanh(x*y)')
figure(2)
ezsurfc('tan(x)*tan(x*y)')
figure(3)
ezsurfc('tanh(x)*tanh(x*y)')
hold on
ezsurfc('tan(x)*tan(x*y)')
%
% 數值解
b = length(0.01:0.01:10);
F = zeros(b,2);
x = 0.01:0.01:10;
for n = 1:b
y = fsolve(@agly,1);
F(n,1) = x(n);
F(n,2) = y;
end
F
plot(F(:,1),F(:,2))
function f=agly(y)
global x n
a = x(n);
f = tanh(a).*tanh(a.*y)++tan(a)*tan(a*y);
------------------------------------------------
我覺得x與y條件可以以負數來討論比較好
--
1.MATLAB programming 2.ASPEN process simulation package
3.FORTRN programming 4.Advance Engineering Mathematics
5.Process Control Theory
6.Chemical Engineering Basic Theory(Kinetic.thermodynamics.transport)
7.Numerical Method and Analysis
8.MATLAB Toolbox.Simulink system basic design
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 113.163.21.45
→
11/12 23:39, , 1F
11/12 23:39, 1F
→
11/13 00:21, , 2F
11/13 00:21, 2F
討論串 (同標題文章)
本文引述了以下文章的的內容：
完整討論串 (本文為第 2 之 2 篇)：