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 算数学题 -- 空间曲线求切线和法向量相关推荐

  1. 大一MATLAB课程设计——空间曲线的切线与法平面

    大一MATLAB课程设计 空间曲线的切线与法平面 课时娱也,课毕遂卒.念及淘宝:忽明.奈何贵也.呜呼,吾魂兮毋求乎永生,当竭尽人事之所能. 先生要求的最简单的gui页面 one,two-seven: ...

  2. 【原创】求空间曲线的切线,法平面

    求空间曲线的切线,法平面:归结为求空间曲线的切向量 进而用点向式直线方程表示出切线,点法式方程表示出法平面 情况一:空间曲线以参数式给出,求切向量时直接求导即可,如下题 情况二:空间曲线以一般式形式给 ...

  3. matlab 曲线法线,平面曲线的切线及法线.ppt

    平面曲线的切线及法线 返回 后页 前页 在本节中所讨论的曲线和曲面, 由于它们 的方程是以隐函数(组)的形式出现的, 因此 在求它们的切线或切平面时, 都要用到隐函 数(组)的微分法. §3 几 何 ...

  4. 导数和梯度,切线和法向量

    转自:http://www.cnblogs.com/jerrylead/archive/2011/03/09/1978280.html 作者讲得很清楚 记得在高中做数学题时,经常要求曲线的切线.见到形 ...

  5. 空间曲线的切线和法平面与曲面的切平面和法线

    (一)空间曲线的切线和法平面 1. 参数方程的形式 理解和记忆如下公式: 参数方程在知道偏导数的情况下,得到该点的切线以及法平面的公式,笔者可以理解但是无法证明. 2. 可以转换为参数方程的第二种形式 ...

  6. matlab中利用xy求取多项式z,matlab基础练习题

    3. 求有理分式()()()()3323230.522521x x x R x x x x ++=+-++的商多项式和余多项式 4. 一元多项式42234p x x x =-+,写出表示p 的MATL ...

  7. Matlab计算微分方程曲线求导及过曲线上点的切线方程

    Matlab计算微分方程曲线求导及过曲线上点的切线方程 求解f(x)=x^2一元二次方程上某点的切线方程并绘制出方程的切线图.点(4,f(4))是曲线方程f(x)上的一个点,求出该点的切线并绘制出来. ...

  8. 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^ ...

  9. 递归 算例一(求一个简单嵌套字典的深度)

    递归 算例一(求一个简单嵌套字典的深度) c=[]def dcc(dic):for key in dic.keys():print (key) c.append(key) vv=dic[key]#判断 ...

最新文章

  1. Gold Code,Gold Sequence
  2. logsoftmax(dim=1)是对批量样本中的每个样本取概率,而不是针对每个特征
  3. WSGI,uWSGI,uwsgi,Nginx
  4. Linux-grep命令
  5. 蓝桥杯java第六届决赛第一题--分机号
  6. LINUX中常用操作命令
  7. virtualbox中Centos6.6和windows共享文件夹设置
  8. Gemstones(牛客第八场多校)
  9. (软件工程复习核心重点)第三章需求分析习题
  10. Windows Server 2008 Server Core模式
  11. UI-148xLE-M相机参数
  12. python学到什么程度可以找到工作-月薪2万+的Python Web岗,学到什么程度能找到工作?...
  13. 视频码率[百科词条]
  14. java图表分析_怎样用图表分析数据(一)
  15. 利用 HttpModule,基于输出,统一控制、干预、处理(例如: 过滤关键字、AntiXSS) ASP.Net WebForm Control 展现属性的方案原型...
  16. 软件项目管理课程设计-数字化校园学工信息系统
  17. 计算机网络中NTFS概念及功能,什么是NTFS
  18. 创建邮件服务器(メールサーバー)Radish
  19. 安防大数据在智慧城市建设中的地位与深度应用
  20. 如何使用计算机打印机共享的打印机驱动,怎样设置打印机共享_不同系统电脑设置打印机共享的方法 - 驱动管家...

热门文章

  1. 【§极品网游之我叫mt online电脑版免费中文版§】
  2. 分书问题(book)-C++
  3. 在cad中出现“代理对象不允许复制 因此写块操作被拒绝”的解决办法
  4. 员工的成长管理——改进、教练与辅导
  5. win10计算机管理器丢失,win10 realtek管理器消失怎么办?
  6. Fluent算例4:脏几何处理的歧管
  7. ETL过程及常用工具简介
  8. 网络营销的方式你知道几种?
  9. paddle第十九期3天训练营3.14-day0~1
  10. 文字检测(白底黑字)