符号运算太慢!

把表达式写出来有太占地方!不过运行快啊!

>> zhouch

ans =

76.4040

>>

function zhouch

quad(@zhchjf,0,2*pi/5)

function y=zhchjf(theta)

z1=6;

z2=z1-1;

R=48.78;

r=8.13;

Q=R-r;

e=7.05;

k=12.6;

f=r/e;

a=sin(z1*theta);

b=f+cos(z1*theta);

m=atan(a./b);

p=m-theta;       %theta

c=sin(theta+p);

n=asin(f*c);

o=n-p;          %phi

g=(z1/z2)*o;    %psi

x2t=Q*sin(theta)+e*sin(z2*theta)-k*sin(p);       %原始齿形

y2t=Q*cos(theta)-e*cos(z2*theta)+k*cos(p);

x3t=x2t.*cos(o-g)-y2t.*sin(o-g)-e*sin(g);        %共轭齿形

y3t=x2t.*sin(o-g)+y2t.*cos(o-g)-e*cos(g);

ds2=(((813./20.*cos(theta)+141./4.*cos(5.*theta)-63./5.*cos(atan(sin(6.*theta)./(271./235+cos(6.*theta)))-theta).*((6.*cos(6.*theta)./(271./235+cos(6.*theta))+6.*sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)-1)).*cos(-1./5.*asin(271./235.*sin(6.*theta)./(271./235+cos(6.*theta))./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(1./2))+1./5.*atan(sin(6.*theta)./(271./235+cos(6.*theta)))-1./5.*theta)-(813./20.*sin(theta)+141./20.*sin(5.*theta)-63./5.*sin(atan(sin(6.*theta)./(271./235+cos(6.*theta)))-theta)).*sin(-1./5.*asin(271./235.*sin(6.*theta)./(271./235+cos(6.*theta))./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(1./2))+1./5.*atan(sin(6.*theta)./(271./235+cos(6.*theta)))-1./5.*theta).*(-47.*(1626./235.*cos(6.*theta)./(271./235+cos(6.*theta))./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(1./2)+1626./235.*sin(6.*theta).^2./(271./235+cos(6.*theta)).^2./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(1./2)-271./470.*sin(6.*theta)./(271./235+cos(6.*theta))./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(3./2).*(12.*sin(6.*theta)./(271./235+cos(6.*theta)).^2.*cos(6.*theta)+12.*sin(6.*theta).^3./(271./235+cos(6.*theta)).^3))./(55225-73441.*sin(6.*theta).^2./(271./235+cos(6.*theta)).^2./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)).^(1./2)+1./5.*(6.*cos(6.*theta)./(271./235+cos(6.*theta))+6.*sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)-1./5)-(-813./20.*sin(theta)+141./4.*sin(5.*theta)-63./5.*sin(atan(sin(6.*theta)./(271./235+cos(6.*theta)))-theta).*((6.*cos(6.*theta)./(271./235+cos(6.*theta))+6.*sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)-1)).*sin(-1./5.*asin(271./235.*sin(6.*theta)./(271./235+cos(6.*theta))./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(1./2))+1./5.*atan(sin(6.*theta)./(271./235+cos(6.*theta)))-1./5.*theta)-(813./20.*cos(theta)-141./20.*cos(5.*theta)+63./5.*cos(atan(sin(6.*theta)./(271./235+cos(6.*theta)))-theta)).*cos(-1./5.*asin(271./235.*sin(6.*theta)./(271./235+cos(6.*theta))./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(1./2))+1./5.*atan(sin(6.*theta)./(271./235+cos(6.*theta)))-1./5.*theta).*(-47.*(1626./235.*cos(6.*theta)./(271./235+cos(6.*theta))./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(1./2)+1626./235.*sin(6.*theta).^2./(271./235+cos(6.*theta)).^2./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(1./2)-271./470.*sin(6.*theta)./(271./235+cos(6.*theta))./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(3./2).*(12.*sin(6.*theta)./(271./235+cos(6.*theta)).^2.*cos(6.*theta)+12.*sin(6.*theta).^3./(271./235+cos(6.*theta)).^3))./(55225-73441.*sin(6.*theta).^2./(271./235+cos(6.*theta)).^2./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)).^(1./2)+1./5.*(6.*cos(6.*theta)./(271./235+cos(6.*theta))+6.*sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)-1./5)-141./20.*cos(6./5.*asin(271./235.*sin(6.*theta)./(271./235+cos(6.*theta))./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(1./2))-6./5.*atan(sin(6.*theta)./(271./235+cos(6.*theta)))+6./5.*theta).*(282.*(1626./235.*cos(6.*theta)./(271./235+cos(6.*theta))./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(1./2)+1626./235.*sin(6.*theta).^2./(271./235+cos(6.*theta)).^2./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(1./2)-271./470.*sin(6.*theta)./(271./235+cos(6.*theta))./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(3./2).*(12.*sin(6.*theta)./(271./235+cos(6.*theta)).^2.*cos(6.*theta)+12.*sin(6.*theta).^3./(271./235+cos(6.*theta)).^3))./(55225-73441.*sin(6.*theta).^2./(271./235+cos(6.*theta)).^2./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)).^(1./2)-6./5.*(6.*cos(6.*theta)./(271./235+cos(6.*theta))+6.*sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)+6./5)).^2+((813./20.*cos(theta)+141./4.*cos(5.*theta)-63./5.*cos(atan(sin(6.*theta)./(271./235+cos(6.*theta)))-theta).*((6.*cos(6.*theta)./(271./235+cos(6.*theta))+6.*sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)-1)).*sin(-1./5.*asin(271./235.*sin(6.*theta)./(271./235+cos(6.*theta))./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(1./2))+1./5.*atan(sin(6.*theta)./(271./235+cos(6.*theta)))-1./5.*theta)+(813./20.*sin(theta)+141./20.*sin(5.*theta)-63./5.*sin(atan(sin(6.*theta)./(271./235+cos(6.*theta)))-theta)).*cos(-1./5.*asin(271./235.*sin(6.*theta)./(271./235+cos(6.*theta))./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(1./2))+1./5.*atan(sin(6.*theta)./(271./235+cos(6.*theta)))-1./5.*theta).*(-47.*(1626./235.*cos(6.*theta)./(271./235+cos(6.*theta))./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(1./2)+1626./235.*sin(6.*theta).^2./(271./235+cos(6.*theta)).^2./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(1./2)-271./470.*sin(6.*theta)./(271./235+cos(6.*theta))./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(3./2).*(12.*sin(6.*theta)./(271./235+cos(6.*theta)).^2.*cos(6.*theta)+12.*sin(6.*theta).^3./(271./235+cos(6.*theta)).^3))./(55225-73441.*sin(6.*theta).^2./(271./235+cos(6.*theta)).^2./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)).^(1./2)+1./5.*(6.*cos(6.*theta)./(271./235+cos(6.*theta))+6.*sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)-1./5)+(-813./20.*sin(theta)+141./4.*sin(5.*theta)-63./5.*sin(atan(sin(6.*theta)./(271./235+cos(6.*theta)))-theta).*((6.*cos(6.*theta)./(271./235+cos(6.*theta))+6.*sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)-1)).*cos(-1./5.*asin(271./235.*sin(6.*theta)./(271./235+cos(6.*theta))./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(1./2))+1./5.*atan(sin(6.*theta)./(271./235+cos(6.*theta)))-1./5.*theta)-(813./20.*cos(theta)-141./20.*cos(5.*theta)+63./5.*cos(atan(sin(6.*theta)./(271./235+cos(6.*theta)))-theta)).*sin(-1./5.*asin(271./235.*sin(6.*theta)./(271./235+cos(6.*theta))./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(1./2))+1./5.*atan(sin(6.*theta)./(271./235+cos(6.*theta)))-1./5.*theta).*(-47.*(1626./235.*cos(6.*theta)./(271./235+cos(6.*theta))./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(1./2)+1626./235.*sin(6.*theta).^2./(271./235+cos(6.*theta)).^2./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(1./2)-271./470.*sin(6.*theta)./(271./235+cos(6.*theta))./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(3./2).*(12.*sin(6.*theta)./(271./235+cos(6.*theta)).^2.*cos(6.*theta)+12.*sin(6.*theta).^3./(271./235+cos(6.*theta)).^3))./(55225-73441.*sin(6.*theta).^2./(271./235+cos(6.*theta)).^2./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)).^(1./2)+1./5.*(6.*cos(6.*theta)./(271./235+cos(6.*theta))+6.*sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)-1./5)+141./20.*sin(6./5.*asin(271./235.*sin(6.*theta)./(271./235+cos(6.*theta))./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(1./2))-6./5.*atan(sin(6.*theta)./(271./235+cos(6.*theta)))+6./5.*theta).*(282.*(1626./235.*cos(6.*theta)./(271./235+cos(6.*theta))./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(1./2)+1626./235.*sin(6.*theta).^2./(271./235+cos(6.*theta)).^2./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(1./2)-271./470.*sin(6.*theta)./(271./235+cos(6.*theta))./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2).^(3./2).*(12.*sin(6.*theta)./(271./235+cos(6.*theta)).^2.*cos(6.*theta)+12.*sin(6.*theta).^3./(271./235+cos(6.*theta)).^3))./(55225-73441.*sin(6.*theta).^2./(271./235+cos(6.*theta)).^2./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)).^(1./2)-6./5.*(6.*cos(6.*theta)./(271./235+cos(6.*theta))+6.*sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)./(1+sin(6.*theta).^2./(271./235+cos(6.*theta)).^2)+6./5)).^2).^(1./2);

y=ds2;

matlab 声明变量范围,matlab中syms如何定义变量范围相关推荐

  1. MATLAB及app designer中函数:定义与调用

    本文介绍MATLAB 及其APP designer 中函数定义和调用的相关方法和注意事项: 文章结构如下: 1.脚本文件(.m)中函数调用: A.同一脚本文件下,函数定义在脚本末端 B.不同脚本文件, ...

  2. matlab中数组的声明,如何使用Matlab声明矩阵数组

    我们在Matlab软件中最常用的就是数组了,那么在Matlab中声明矩阵数组有哪些方式呢?下面小编给大家分享总结一下. 工具/材料 Matlab 首先在Matlab中最常用的声明数组的方式就是使用中括 ...

  3. Matlab中的类定义 classdef

    Matlab也可以写面向对象的代码,首先表现在可以定义类,以及可以继承,使用类(class)有很多好处,其中一个重要的好处便是解决变量名冲突和让函数.对象的结构清晰.class的static func ...

  4. matlab中函数的定义和使用

    在matlab的选择框中的的新建中选择新建函数,一般是以.m的文本文件,新建后我们可以在里面写入我们的实现函数,现在我们给出两个实例看一下,分别是简单的求和和Fibonacci sequence. 1 ...

  5. matlab定义字母常数,在matlab中,如何定义函数式子中的未知常数。

    新建一个.m文件,第一行输入下面的格式function [y,m] = abc(x)其中... 定义未知数:syms xx是变量 .syms是定义符号变量.补充:syms是... 这不是解方程,这是曲 ...

  6. matlab hilb,MATLAB在科学计算中的应用:第2章 MATLAB 语言程序设计基础

    <MATLAB在科学计算中的应用:第2章 MATLAB 语言程序设计基础>由会员分享,可在线阅读,更多相关<MATLAB在科学计算中的应用:第2章 MATLAB 语言程序设计基础(7 ...

  7. Matlab在概率统计中的应用问题及解决方案集锦

    前言 关于MATLAB系列的精品专栏大家可参见 MATLAB-30天带你从入门到精通 MATLAB深入理解高级教程(附源码) 喜欢的小伙伴可自行订阅,你的支持就是我不断更新的动力哟! Matlab在概 ...

  8. 热力学分布用matlab,matlab在热物理学中的应用.doc

    matlab在热物理学中的应用.doc 1MATLAB在热物理学中的应用摘要本文阐述了基于MATLAB的数值计算.可视化图形处理.开放式以及可扩充体系结构的特点,并介绍了高性能语言MATLAB在大学物 ...

  9. Matlab基础学习(中)

    1.前记,今天的继续学习知乎专栏:https://zhuanlan.zhihu.com/c_1208096335112843264 所有代码:学习建议,复制粘贴到matlab脚本,保存命名自定,按小结 ...

最新文章

  1. Bzoj4016: [FJOI2014]最短路径树问题
  2. 服务器显示接口类型,查看服务器各接口卡情况的命令是:lspci
  3. 痛并快乐的造轮子之旅:awk访问数据库之旅
  4. 基于深度学习的语义分割代码库
  5. hadoop+hive+flink+hbase交互的版本兼容性
  6. python安装运行时提示不是内部或外部命令怎么办_如何解决cmd运行python提示不是内部命令...
  7. ABP入门系列(11)——编写单元测试
  8. 阿里云 超级码力在线编程大赛初赛 第2场 题目1. 三角魔法
  9. 容器编排技术 -- Kubernetes 为 Namespace 配置Pod配额
  10. 冬天来了,温暖甜品热饮海报设计psd模板,勾住你的胃!
  11. TensorFlow 支持 Unicode,中文 NLP 终于省心了
  12. [LeetCode]Trapping Rain Water
  13. Android获取md5和sha
  14. 【入门】QSS基础入门笔记
  15. easydarwin 安装_win10安装EasyDarwin并作为流媒体服务器的推流和拉流实例
  16. 产品分析 | K12在线教育之猿辅导
  17. 解决小程序直播组件live-player全屏问题
  18. 厦门大学计算机系录取分数线贵州,贵州省多少名可以进厦门大学?附厦门大学近三年录取分数线...
  19. zencart模板修改 (详细)
  20. 场景识别论文阅读感想(初步)

热门文章

  1. 找情感挽回大师拯救情感的究竟是些什么人
  2. MATLAB产生阶跃信号
  3. C语言GDF中是什么作用,深度解密格力GDF-40X21C怎么样呢?质量好不好?真相揭秘必看...
  4. 从零开始学习Memcached之stats sizes命令测试
  5. 概率图模型 —— 串连 NB、LR、MEM、HMM、CRF
  6. 【翻译】使用Modernizr探测HTML5/CSS3新特性
  7. openWRT摄像头选购
  8. 【汇智学堂】-JS菜鸟版过河小游戏(之四---全部角色左侧上船与左侧下船)
  9. qt报错意外标记的class
  10. 联想A798T刷机包 基于百度云V6 集成RE3.1.7美化版 精简冗余文件