Matlab求解三角函数的周期
来源:学生作业帮 编辑:搜搜考试网作业帮 分类:综合作业 时间:2024/07/28 10:41:23
Matlab求解三角函数的周期
如题,f(x)=a1*sin(b1*x+c1) + a2*sin(b2*x+c2) + a3*sin(b3*x+c3) + a4*sin(b4*x+c4) + a5*sin(b5*x+c5) + a6*sin(b6*x+c6) + a7*sin(b7*x+c7) + a8*sin(b8*x+c8).令f(x)=f(x+T).求T,即f(x)的周期.全部参数都有具体的值,恕不列出.请高手给出Matlab程序,附上T的表达式.
请看仔细:f(x)是由8个sin函数组成的~b1----b8分别是4172,19690,5403,21250,14900,39460,27360,44120.
如题,f(x)=a1*sin(b1*x+c1) + a2*sin(b2*x+c2) + a3*sin(b3*x+c3) + a4*sin(b4*x+c4) + a5*sin(b5*x+c5) + a6*sin(b6*x+c6) + a7*sin(b7*x+c7) + a8*sin(b8*x+c8).令f(x)=f(x+T).求T,即f(x)的周期.全部参数都有具体的值,恕不列出.请高手给出Matlab程序,附上T的表达式.
请看仔细:f(x)是由8个sin函数组成的~b1----b8分别是4172,19690,5403,21250,14900,39460,27360,44120.
![Matlab求解三角函数的周期](/uploads/image/z/16993548-36-8.jpg?t=Matlab%E6%B1%82%E8%A7%A3%E4%B8%89%E8%A7%92%E5%87%BD%E6%95%B0%E7%9A%84%E5%91%A8%E6%9C%9F)
a1*sin(b1*x+c1) 周期为 T1=2*pi/b1;
...
a8*sin(b8*x+c8) 周期为 T8=2*pi/b8;
f的周期为T1,...,T8的最小公倍数
即为2*pi/gcd(b1,...,b8)
gcd(b1,...,b8)=1
故T=2*pi;
a=[4172,19690,5403,21250,14900,39460,27360,44120];
g=a(1);
for i=1:8
g=gcd(g,a(i))
end
disp(2*pi/g)
...
a8*sin(b8*x+c8) 周期为 T8=2*pi/b8;
f的周期为T1,...,T8的最小公倍数
即为2*pi/gcd(b1,...,b8)
gcd(b1,...,b8)=1
故T=2*pi;
a=[4172,19690,5403,21250,14900,39460,27360,44120];
g=a(1);
for i=1:8
g=gcd(g,a(i))
end
disp(2*pi/g)