图形几何——三角形：三角函数小角度近似的成立条件

悠以悠

已于 2024-04-16 13:02:49 修改

阅读量1.7k

点赞数 3

分类专栏：图形几何文章标签：线性代数

于 2023-08-13 11:29:09 首次发布

本文链接：https://blog.csdn.net/chaser30/article/details/132257869

版权

图形几何专栏收录该内容

7 篇文章 0 订阅

订阅专栏

文章目录

小角度近似成立边界(误差小于1%)

小角度近似成立边界(误差小于1%)

$cos\theta \approx 1$ ， $\theta <= 8.07 °$ 时成立

对比结果图

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-ocIrArLE-1691896890813)(https://note.youdao.com/yws/res/62410/WEBRESOURCE52ce5dcf92c7e0e2799926284d0ebe4d)]

matlab代码

clc;clear;close all
x = 0:0.01:15;
y_real = cosd(x);
y_hat = ones(1,size(x,2));
err = (y_hat - y_real)./y_real.*100;
index = find(err > 0.997 & err< 1.002);
figure;plot(x,err);hold on;plot(x(index),err(index),'ro');
hold on; plot(x,y_hat);title('1 - cos x');
ylabel('Relat error, %');xlabel('Angle，°');
figure;plot(x,y_real,'b-');hold on;
plot(x,y_hat,'r-');legend('y = cos x','y = 1','Location','best')
xlabel('Angle，°');ylabel('y');title('y = cosx')

$\theta \approx \theta$ ， $\theta <= 9.91 °$ 时成立

对比结果图

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-yEgUuKDc-1691896890814)(https://note.youdao.com/yws/res/62420/WEBRESOURCE4a96f1016f2d116a870ead8399cabf48)]

matlab代码

clc;clear;close all
x = 0:0.01:45;
y_real = tand(x);
y_hat = deg2rad(x);
figure;plot(x,y_real,'b-');hold on;
plot(x,y_hat,'r-');legend('y = tan x','y = x(rad)','Location','best')
xlabel('Angle，°');ylabel('y');title('y = tanx')
err = (y_real - y_hat)./y_real.*100;
index = find(err > 0.998 & err< 1.002);
figure;plot(x,err);hold on;plot(x(index),err(index),'ro');
hold on; plot(x,ones(1,size(x,2)));title('tan x - x');
ylabel('Relat error, %');xlabel('Angle，°');

$\theta \approx \theta$ ， $\theta <= 13.99 °$ 时成立

对比结果图

在这里插入图片描述

matlab代码

clc;clear;close all
x = 0:0.01:90;
y_real = sind(x);
y_hat = deg2rad(x);
figure;plot(x,y_real,'b-');hold on;
plot(x,y_hat,'r-');legend('y = tan x','y = x(rad)','Location','best')
xlabel('Angle，°');ylabel('y');title('y = sinx')
err = (y_hat - y_real)./y_real.*100;
index = find(err > 0.998 & err< 1.002);
figure;plot(x,err);hold on;plot(x(index),err(index),'ro');
hold on; plot(x,ones(1,size(x,2)));title('x - sinx');
ylabel('Relat error, %');xlabel('Angle，°');

$\theta \approx 1 - \theta^2/2$ ，` $\theta <= 37.93 °$ 时成立

对比结果图

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-LFIwoQaJ-1691896890815)(https://note.youdao.com/yws/res/62433/WEBRESOURCE485bdefcaa25b97524f273620789e96d)]

matlab代码

clc;clear;close all
x = 0:0.01:90;
y_real = cosd(x);
y_hat = 1 - deg2rad(x).^2./2;
figure;plot(x,y_real,'b-');hold on;
plot(x,y_hat,'r-');
legend('y = cos x','y = 1-x^2/2(rad)','Location','best')
xlabel('Angle，°');ylabel('y');title('y = cosx')
err = abs((y_hat - y_real)./y_real.*100);
index = find(err > 0.999 & err< 1.002);
figure;plot(x,err);hold on;plot(x(index),err(index),'ro');
hold on; plot(x,ones(1,size(x,2)));title('cosx(°) - (1-x^2/2(rad))');
ylabel('Relat error, %');xlabel('Angle，°');
ylim([0,3])