这部分重新从定积分学了

1,lnx 的导数就是x^(-1) = 1/x

那么求∫(1/x)dx = ln|x|+C

 2,初值问题.就是求∫f(x)dx = F(x) + C 求C

.

3,Houdini体积映射:

#define PI 3.1415926

vector gmin ;
vector gmax ;
getbbox(,gmin,gmax); // get sphere radius
float R = gmax.x - gmin.x;
//printf("gmax.x / gmaxmin : %f,%f \n", gmax.x, gmin.x);
float r = R / 2.000f; vector center = (gmax + gmin)/;
int segments = ;
float dx = R / segments; function float cylinder_volume(float r ; float h)
{
return PI * r *r * h;
} float v = 0.0f;
for(int i=;i<segments;i++)
{
float ci = lerp(-r, r ,1.0f/segments * i);
float cylinder_r = sqrt(r*r - ci*ci);
//printf("r/ci/cylinder_r %f/%f/%f \n" ,r, ci , cylinder_r);
float cl_volume = cylinder_volume(cylinder_r, dx);
v += cl_volume;
}
//printf("volume is %f \n" , v);
setdetailattrib(geoself(),"volume",v);

下图可以看到体积小的是黑色.虽然方法积的是圆柱形,最终为球,为了学习积分,没用球的方程

两条曲线之间的面积:

A = ∫ [f(x)-g(x) ] dx  积分上限b,下限a

数值积分: 

有时候不知道反导数怎么求?梯形法,simpson法

曲线积分:求曲线长度:

求普通曲线的长度:

L = ∫ [sqrt(1+dx/dy)^2 ] dy 或者 L=∫ [sqrt(1+dy/dx)^2] dx

求参数方程如图:

x= cos^3t

y = sin^3t

0<= t <= 2π

解方法就是:

√ [(dx/dt)^2 + (dy/dt)^2]

float pi = 3.141596;
float t = fit(@ptnum,,@numpt-,,) * * pi;
@P.x = cos(t) * cos(t) * cos(t);
@P.y = sin(t) * sin(t) * sin(t);

HoudiniCode

蒙特卡洛求PI

曲线下的面积/矩形面积 = 曲线下的总点数/随机点的总数

<1>求PI

import random,math

import numpy as np
import matplotlib.pyplot as plt N = 100000
distance = lambda x, y: math.sqrt(x * x + y * y) np.random.seed(100)
x = np.random.uniform(-1,1,N)
np.random.seed(200)
y = np.random.uniform(-1,1,N) def Monte_carlo_PI():
hitNum = 0
for index in xrange(0,N,1):
if distance(x[index],y[index]) <=1 :
hitNum += 1
proportion_seed = float(hitNum)/float(N)
return proportion_seed * 4 if __name__ == "__main__":
print Monte_carlo_PI() fig = plt.figure(figsize=(10,10)) ax = fig.add_subplot(111)
ax.spines['top'].set_color('none')
ax.spines['right'].set_color('none')
ax.xaxis.set_ticks_position('bottom')
ax.spines['bottom'].set_position(('data', 0))
ax.yaxis.set_ticks_position('left')
ax.spines['left'].set_position(('data', 0)) #plt.plot(xs,ys,'r')
#plt.scatter(0, 0, s=100000)
circle2 = plt.Circle((0, 0), 1, color='r',alpha = 0.5)
ax.add_patch(circle2) plt.scatter(x, y,s=1,c="black")
plt.grid(True)
plt.show()

python_mente_claro

功:

弹簧做功:弹簧的积分区域一定是变化了多少,从0开始积分变化了多少。如下

1,弹簧系数: 把弹簧从自然长度2米拉到5米做功为1800焦耳,求弹簧的弹力系数.

弹簧力是:F = kx (系数*变化的距离)

弹簧做功:∫F(x)dx  积分阈为[0,变化长度]

所以W ∫F(x)dx = k/2 * x^2    [0,3]

W(3) - W(0) = 1800 ,k = 400牛顿/米

抽水做功:

第一步(薄片体积)        :       ΔV = PI * r^2  Δy                  这里薄片体积表示为圆柱

第二步(薄片力)            :      F(y) = ΔV * 比重

第三步(薄片需要的功)   :    ΔW = 力 X 力的作用距离 = F(y) *   (距离-y)

第四步(黎曼和总功)      :    w =Σ ΔW

矩和质心

....

积分从入门到放弃<2>的更多相关文章

  1. CYQ.Data 从入门到放弃ORM系列:开篇:自动化框架编程思维

    前言: 随着CYQ.Data 开始回归免费使用之后,发现用户的情绪越来越激动,为了保持这持续的激动性,让我有了开源的念头. 同时,由于框架经过这5-6年来的不断演进,以前发的早期教程已经太落后了,包括 ...

  2. [精品书单] C#/.NET 学习之路——从入门到放弃

    C#/.NET 学习之路--从入门到放弃 此系列只包含 C#/CLR 学习,不包含应用框架(ASP.NET , WPF , WCF 等)及架构设计学习书籍和资料. C# 入门 <C# 本质论&g ...

  3. OpenStack从入门到放弃

    OpenStack从入门到放弃 目录: 为何选择云计算/云计算之前遇到的问题 什么是云计算 云服务模式 云应用形式 传统应用与云感知应用 openstack及其相关组件介绍 flat/vlan/gre ...

  4. 绕过校园网的共享限制 win10搭建VPN服务器实现--从入门到放弃

    一.开篇立论= =.. 上次说到博主在电脑上搭建了代理服务器来绕过天翼客户端的共享限制,然而经过实际测试还不够完美,所以本着生命不息,折腾不止的精神,我又开始研究搭建vpn服务器= =... (上次的 ...

  5. 《区块链:从入门到放弃》之obc安装步骤

    obc安装步骤 朋友们可能会好奇,厨师不研究菜谱怎么改研究兵法了,哈哈,我原本是app出身,最近被安排去预研区块链和比特币技术,2个月下来,颇有斩获.期间得到IBM的CC同学指导我一步一步安装obc的 ...

  6. win10搭建代理服务器实现绕过校园网的共享限制--从入门到放弃

    博主所在学校特别坑爹,校园网被电信一家垄断了,而且最恶心的还是电信要求一条网线只能供一台电脑上网,不许接路由器共享网络= =- (还有电信2M价格是380+每年,20m是500每年,而且网速都很慢= ...

  7. WPF从入门到放弃系列第二章 XAML

    本文是作者学习WPF从入门到放弃过程中的一些总结,主要内容都是对学习过程中拜读的文章的整理归纳. 参考资料 XAML 概述 (WPF):https://msdn.microsoft.com/zh-cn ...

  8. Android -- 带你从源码角度领悟Dagger2入门到放弃

    1,以前的博客也写了两篇关于Dagger2,但是感觉自己使用的时候还是云里雾里的,更不谈各位来看博客的同学了,所以今天打算和大家再一次的入坑试试,最后一次了,保证最后一次了. 2,接入项目 在项目的G ...

  9. Android -- 带你从源码角度领悟Dagger2入门到放弃(二)

    1,接着我们上一篇继续介绍,在上一篇我们介绍了简单的@Inject和@Component的结合使用,现在我们继续以老师和学生的例子,我们知道学生上课的时候都会有书籍来辅助听课,先来看看我们之前的Stu ...

随机推荐

  1. vim: 基本知识;

    1. 函数: function!   funcName(para.) content; endfunction 如果添加!,将覆盖已存在的重名函数: 注: 该博文为扩展型: 2.调用外部命令: exe ...

  2. jmeter计算身份证校验位

    idcard_no='111111198101017000' ; idcard_no_array = idcard_no.substring(0,17).toCharArray(); int[] c ...

  3. 剑指Offer_编程题_3

    题目描述 输入一个链表,从尾到头打印链表每个节点的值. /** * struct ListNode { * int val; * struct ListNode *next; * ListNode(i ...

  4. SQL语法基础之CREATE语句

    SQL语法基础之CREATE语句 作者:尹正杰 版权声明:原创作品,谢绝转载!否则将追究法律责任. 一.查看帮助信息 1>.使用“?”来查看MySQL命令的帮助信息 mysql> ? CR ...

  5. Kafka吞吐量测试案例

    Kafka吞吐量测试案例 作者:尹正杰 版权声明:原创作品,谢绝转载!否则将追究法律责任. 领英公司参考连接:https://www.slideshare.net/JiangjieQin/produc ...

  6. Centos6.6搭建Maven私服

    操作系统:Centos6.6 私服Ip:10.0.210.112 JDK:1.7 (已安装并配置好了环境变量) 1:上 传 nexus-2.11.2-03-bundle.tar.gz到/root/ne ...

  7. Hadoop记录-技术网站

    Grafan监控:http://docs.grafana.org/ Tez:http://tez.apache.org/install.html 阿里巴巴镜像:https://opsx.alibaba ...

  8. Hbuilder开发app时生成ios要的mobileprovision和p12文件步骤.

    1.在MAC电脑.钥匙串串访问->证书助理->从证书颁发机构请求证书,创建一个证书为certSigningRequest文件 2.在Apple Developer中的Certificate ...

  9. 鼠标右键Table的td弹出多级菜单,双击td编辑

    <html xmlns="http://www.w3.org/1999/xhtml"> <head> <meta http-equiv="C ...

  10. Emacs 快速指南(中文翻译)

      Emacs 快速指南 目录 1. 小结(SUMMARY) 2. 基本的光标控制(BASIC CURSOR CONTROL) 3. 如果 EMACS 失去响应(IF EMACS STOPS RESP ...