正文

积分从入门到放弃

AlgebraMaster  AlgebraMaster  2022-10-01  234

关键词：

这部分重新从定积分学了

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(1,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)/2;
int segments = 20;
float dx = R / segments;

function float cylinder_volume(float r ; float h)
{
    return PI * r *r * h;
}

float v = 0.0f;
for(int i=0;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);

View Code

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

两条曲线之间的面积:

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,0,@numpt-1,0,1) * 2 * 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

矩和质心

....

java从入门到放弃

（1）Java入门变量与注释用户输入if语句java中的循环while 查看详情

深度学习:从入门到放弃

https://zhuanlan.zhihu.com/p/22976342 首发于深度学习:从入门到放弃写文章登录 FCN学习:SemanticSegmentation余俊1年前感谢@huangh12 @郑途 @麦田守望者对标签图像生成的研究和讨论，这几天研究了一下，补充如下。-------------------... 查看详情

[资源]深度学习从入门到放弃

Relationship:　　MachineLearning---->DeepLearning　　　　　　　　　　　　　　　　　　　　　　　　　　---->DeepReinforcementLearning[LearningRoadMap]　　　　　　　　　　　　　　ReinforcementLearningPapers:　　DeepLearningPapersReadin 查看详情

vue从入门到放弃(代码片段)

----------------------------------------------------点击这里《专栏目录》查看更多--------------------------------------------------------------------------------------------------------点击这里《专栏目录》查看更多---------------- 查看详情

vue从入门到放弃(代码片段)

凸优化从入门到放弃（目录）

...课程时候的笔记和心得，其中应用篇没有学习。凸优化从入门到放弃00-凸优化引言01-凸集02-凸函数（暂无）03-凸优化问题（暂无）04-对偶（暂无）05-无约束优化算法（暂无）06-等式约束优化算法（暂无）07-内点法(不等式约束优... 查看详情

visualstudioformac从入门到放弃1

MAC 第一步：从微软官网下载：https://www.visualstudio.com/vs/visual-studio-mac/第二步：安装软件过程出现Itwasnotpossibletocompleteanautomaticinstallation.Thismightbeduetoaproblemwithyournetwork,proxyserversoranunsolvablei 查看详情

kubelet从入门到放弃：拓扑管理(上)

<Kubelet从入门到放弃>系列将对Kubelet组件由基础知识到源码进行深入梳理。上一篇zouyee带各位看了CPU管理的相关内容，其中提及拓扑管理，本文将对此进行详细剖析,拓扑管理在Kubernetes1.18时提升为Beta。TopologyManager功能可实现C... 查看详情

观察者模式（从放弃到入门）

观察者模式（从放弃到入门）观察者模式从放弃到入门需求1错误示范需求2观察者模式认识观察者模式需求3Java中自带的观察者模式思考源码ObservableObserver解析一些体会今天分享第二个模式，观察者模式。相信作Android... 查看详情

webpack从入门到放弃之路

公司的中流砥柱要走啦!!!!我要接手这些摊子啦!!!!!硬着头皮上吧!/(ㄒoㄒ)/~~第一部分:webpack使用部分第二部分:自动化部署部分第三部分:前端视频部分tobecontinue… 查看详情

初识pytorch：从安装到入门，从入门到放弃(代码片段)

目录PyTorch安装配置安装验证PyTorchPyTorch是Facebook团队于2017年1月发布的一个深度学习框架，虽然晚于TensorFlow，也没有TensorFlow火，但目前已经与TensorFlow奇虎相当。而且PyTorch采用了Python语言的接口，可以说它才是Pytho... 查看详情

[补档][从入门到放弃]——网络流学习索引

...啊喂）最长的博文2017-7-29大佬讲课笔记网络流——从入门到放弃没错，就是它==各种题解的索引搭配飞行员士兵占领王者之剑happiness切糕餐巾血帆海盗查看详情

vue从入门到放弃(代码片段)

vue2.0创建项目命令：(base)goldwater@hujinshuivue%cnpminstall-gvue-cli(base)goldwater@hujinshuivue%vueinitwebpackvue2.0-demo?Projectnamevue2.0-demo#项目名称?ProjectdescriptionAVue.jsproject#项目描述?Au 查看详情

深度学习---从入门到放弃pytorch基础(代码片段)

深度学习—从入门到放弃(一)pytorchTensor类似于numpy的array，pandas的dataframe；在pytorch里的数据结构是tensor，即张量tensor简单操作1.Flattenandreshape###Originalz:tensor([[0,1],[2,3],[4,5],[6,7],[8,9],[10,11]])Flatte 查看详情

转-vue.js2.0从入门到放弃---入门实例

weex从入门到放弃

一、安装(1)npminstall-gweex-toolkit(2)weexinitawesome-project二、运行项目npmrundev:监测文件变化，执行package.json中的script下的dev，实际执行watch指定npmrunserve:启动服务器两者必须要一起运行，不然每次更改文件都得执行npmrunbuild，先编译，... 查看详情

volatile从入门到放弃

1.引言如果你对java的volatile有着疑虑请阅读本文,如果你想对volatile想有一个更深的了解,请阅读本文.本文主要讲的是volatile的写happen-before在volatile读之前所涉及相关的原理,以及在Hotspot中相关代码的实现. ... 查看详情