过圆外一点与圆相切的直线

文章目录


1.问题描述

过定点与圆相切的直线,一般分为三种情况,a 定点在圆上(一个解)。b 定点在圆内(无解)。c 定点在圆外(两个解)。前面一节讨论了过圆上一点与圆相切的直线问题,该文接着讨论过圆外的一点,与圆相切的情况。

2.已知条件

圆心C,半径r,圆外点P。求过点P与圆相切的直线。Q为切点。

3.求解

4.code

matlab 复制代码
clc;
close all;
clear;

%% Point P on the line is outside circle
rng default;
C = rand(1,2);
r = rand()*3;
theta = rand()*2*pi;
r_out = r + rand();
P = C+[r_out*cos(theta), r_out*sin(theta)];
u = P - C;
r_2 = r*r;
u_2 = u.^2;
delta = sqrt(-r_2^2*u_2(1) + r_2*u_2(1)^2 + r_2*u_2(1)*u_2(2));
x1 = (r_2 - r_2*u_2(2)/sum(u_2) - u(2)*delta/sum(u_2))/u(1);
y1 = (r_2*u(2) + delta)/sum(u_2);

x2 = (r_2 - r_2*u_2(2)/sum(u_2) + u(2)*delta/sum(u_2))/u(1);
y2 = (r_2*u(2) - delta)/sum(u_2);



u = linspace(0,2*pi,50);
circle = C + [r*cos(u'), r*sin(u')];
t = linspace(-5,5,50);
n = [-y1, x1];
L = P + t'*n;
n2 = [-y2, x2];
L2 = P + t'*n2;


figure;
plot(circle(:,1),circle(:,2));
hold on
plot(L(:,1), L(:,2));
hold on
plot(L2(:,1), L2(:,2));
hold on
scatter(P(1), P(2),'filled','o')
axis equal
grid on
相关推荐
Eric.Lee202119 分钟前
数据集-目标检测系列- F35 战斗机 检测数据集 F35 plane >> DataBall
人工智能·算法·yolo·目标检测·计算机视觉
Gsen281930 分钟前
AI大模型从0到1记录学习 数据结构和算法 day20
数据结构·学习·算法·生成对抗网络·目标跟踪·语言模型·知识图谱
一定要AK30 分钟前
天梯——L1-110 这不是字符串题
数据结构·c++·算法
2401_858286112 小时前
E47.【C语言】零散的练习题(1)
c语言·数据结构·算法·指针
凯子坚持 c2 小时前
深度解析之算法之分治(快排)
算法·leetcode·职场和发展
✿ ༺ ོIT技术༻4 小时前
笔试强训:Day2
开发语言·c++·笔记·算法
飞桨PaddlePaddle6 小时前
Wan2.1和HunyuanVideo文生视频模型算法解析与功能体验丨前沿多模态模型开发与应用实战第六期
人工智能·算法·百度·音视频·paddlepaddle·飞桨·deepseek
Starry_hello world6 小时前
C++ 快速幂算法
c++·算法·有问必答
石去皿7 小时前
力扣hot100 91-100记录
算法·leetcode·职场和发展
SsummerC9 小时前
【leetcode100】组合总和Ⅳ
数据结构·python·算法·leetcode·动态规划