Matlab 算数学题 -- 空间曲线求切线和法向量
clc; clear;
syms x y z t
%t = linspace(0,2*pi,100);
x = t-sin(t);
y = 1-cos(t);
z = 4.*sin(t./2);
S = jacobian([x,y,z],t) % 曲线方程的切向量(得到x,y,z关于t的导数方程)
t = pi/2;
S0 = subs(S) % 切点处的x,y,z方向分别的切向量
x0 = t-sin(t) % 切点(x0,y0,z0)
y0 = 1-cos(t)
z0 = 4.*sin(t./2)
F = -[x;y;z]+[x0;y0;z0]+S0*t %切线方程
G = [x-x0,y-y0,z-z0]*S0 %法平面方程
figure
%axis square %坐标轴尺寸正方形
%set(gca,'xtick',[],'ytick',[],'xcolor','w','ycolor','w') %隐藏坐标轴
t = linspace(0,2*pi,100);
x = eval(x);
y = eval(y);
z = eval(z);
plot3(x,y,z)
grid on
xlabel('x')
ylabel('y')
zlabel('z')
hold on
u = -1:0.01:1;
x_Tangent = u.*S0(1)+x0;
y_Tangent = u.*S0(2)+y0;
z_Tangent = u.*S0(3)+z0;
plot3(x_Tangent,y_Tangent,z_Tangent) % 切线
xlab = -1:0.1:1;
[x_normalPlane,y_normalPlane] = meshgrid(xlab);
z_normalPlane = double((-S0(1).*(x_normalPlane-x0)-S0(2).*(y_normalPlane-y0))/S0(3))+z0;
mesh(x_normalPlane,y_normalPlane,z_normalPlane)
Matlab 算数学题 -- 空间曲线求切线和法向量相关推荐
- 大一MATLAB课程设计——空间曲线的切线与法平面
大一MATLAB课程设计 空间曲线的切线与法平面 课时娱也,课毕遂卒.念及淘宝:忽明.奈何贵也.呜呼,吾魂兮毋求乎永生,当竭尽人事之所能. 先生要求的最简单的gui页面 one,two-seven: ...
- 【原创】求空间曲线的切线,法平面
求空间曲线的切线,法平面:归结为求空间曲线的切向量 进而用点向式直线方程表示出切线,点法式方程表示出法平面 情况一:空间曲线以参数式给出,求切向量时直接求导即可,如下题 情况二:空间曲线以一般式形式给 ...
- matlab 曲线法线,平面曲线的切线及法线.ppt
平面曲线的切线及法线 返回 后页 前页 在本节中所讨论的曲线和曲面, 由于它们 的方程是以隐函数(组)的形式出现的, 因此 在求它们的切线或切平面时, 都要用到隐函 数(组)的微分法. §3 几 何 ...
- 导数和梯度,切线和法向量
转自:http://www.cnblogs.com/jerrylead/archive/2011/03/09/1978280.html 作者讲得很清楚 记得在高中做数学题时,经常要求曲线的切线.见到形 ...
- 空间曲线的切线和法平面与曲面的切平面和法线
(一)空间曲线的切线和法平面 1. 参数方程的形式 理解和记忆如下公式: 参数方程在知道偏导数的情况下,得到该点的切线以及法平面的公式,笔者可以理解但是无法证明. 2. 可以转换为参数方程的第二种形式 ...
- matlab中利用xy求取多项式z,matlab基础练习题
3. 求有理分式()()()()3323230.522521x x x R x x x x ++=+-++的商多项式和余多项式 4. 一元多项式42234p x x x =-+,写出表示p 的MATL ...
- Matlab计算微分方程曲线求导及过曲线上点的切线方程
Matlab计算微分方程曲线求导及过曲线上点的切线方程 求解f(x)=x^2一元二次方程上某点的切线方程并绘制出方程的切线图.点(4,f(4))是曲线方程f(x)上的一个点,求出该点的切线并绘制出来. ...
- matlab怎么对sinx求导,用matlab程序求y=ln(sinx 1)的导数
哪位高手会画函数f(x,y)=(1+e^y)*cosx-y*e^y的图像啊?求真相,最好是用Matlab并附上程序! clearall;clc;%画函数f(x,y)=(1+e^y)*cosx-y*e^ ...
- 递归 算例一(求一个简单嵌套字典的深度)
递归 算例一(求一个简单嵌套字典的深度) c=[]def dcc(dic):for key in dic.keys():print (key) c.append(key) vv=dic[key]#判断 ...
最新文章
- Gold Code,Gold Sequence
- logsoftmax(dim=1)是对批量样本中的每个样本取概率,而不是针对每个特征
- WSGI,uWSGI,uwsgi,Nginx
- Linux-grep命令
- 蓝桥杯java第六届决赛第一题--分机号
- LINUX中常用操作命令
- virtualbox中Centos6.6和windows共享文件夹设置
- Gemstones(牛客第八场多校)
- (软件工程复习核心重点)第三章需求分析习题
- Windows Server 2008 Server Core模式
- UI-148xLE-M相机参数
- python学到什么程度可以找到工作-月薪2万+的Python Web岗,学到什么程度能找到工作?...
- 视频码率[百科词条]
- java图表分析_怎样用图表分析数据(一)
- 利用 HttpModule,基于输出,统一控制、干预、处理(例如: 过滤关键字、AntiXSS) ASP.Net WebForm Control 展现属性的方案原型...
- 软件项目管理课程设计-数字化校园学工信息系统
- 计算机网络中NTFS概念及功能,什么是NTFS
- 创建邮件服务器(メールサーバー)Radish
- 安防大数据在智慧城市建设中的地位与深度应用
- 如何使用计算机打印机共享的打印机驱动,怎样设置打印机共享_不同系统电脑设置打印机共享的方法 - 驱动管家...