积分从入门到放弃

Posted 2020-10-15 AlgebraMaster

tags:

篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了积分从入门到放弃相关的知识，希望对你有一定的参考价值。

这部分重新从定积分学了

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()