2017年12月_zhangphil

12月 11月 10月 09月 08月 07月 06月 05月 04月 03月 02月 01月

原创黎曼积分求解可微曲线的弧线长度

黎曼积分求解可微曲线的弧线长度假设曲线y=f(x)在区间[a,b]内光滑、可微且连续。那么可以根据微积分求解y=f(x)在a如图：从微分的思想入手建立数学函数式，假设s为曲线上(x,f(x))到(x+dx,f(x+dx))两点连线。这两点在水平方向的长度为dx，在垂直方向的y坐标轴长度为dy，根据直角三角形的勾股定理可知：其中，由f’(x)=dy/dx，得到dy =f’(x) dx从而：即ds的长

2017-12-31 00:22:15 5558

原创 MATLAB计算黎曼积分曲线围成的面积

MATLAB计算黎曼积分曲线围成的面积假设一个曲线方程f(x)= x.^3-x.^2-2*x。f(x)与笛卡尔坐标x坐标轴有交点，如图：计算该曲线与x（1设所求面积为S，那么：但是f(x)与x坐标轴相交形成的两块面积，在x区域[-1,0]为正，[0,2]为负，因此要对[0,2]区域的面积分开计算，分别为：与然后取绝对值相加。matlab: syms x f;f=x.^3-x.^2-2*x;li

2017-12-29 20:11:01 14270

原创 Newton-Raphson切线法解高次方程近似根

Newton-Raphson切线法解高次方程近似根对于一般的一次，二次方程来说，求解方程的根比较简单。但是对于四次、五次甚至更高次方程，求解方程的f(x)=0的根变得十分困难甚至不可能完成。为此Newton（牛顿）在1736年 Method of Fluxions 中发表文章提出一种解决方案，事实上，牛顿所提出的这种方案，另一位数学家Joseph Raphson于1690年已经发现。为此，牛顿法

2017-12-27 15:58:54 7238

原创线性化微分数学解释Einstein狭义相对论质能方程E=MC^2

线性化微分数学解释Einstein狭义相对论质能方程E=MC^2要理解爱因斯坦在狭义相对论中的质能方程是如何推导出来的，需要先了解数学中的微分方程及其线性化方程的知识。现在先从最简单的微分方程开始。以简单的曲线方程y=x*x和它的切线方程为例。假设取y=x*x上一点(1,1)，过(1,1)点的切线方程很容易求得，根据一般的过曲线上点（a,f(a)）的切线方程公式：f(x)=f’(a)(x-a)+f

2017-12-25 20:46:28 6783 14

原创经济金融领域简单数学建模和分析：MATLAB成本曲线方程和销售收入直线方程

经济金融领域简单数学建模和分析：MATLAB成本曲线方程和销售收入直线方程MATLAB代码：x=[0:0.1:5]; y=9*x; plot(x,y,'r','LineWidth',0.5)hold on;y=x.^3-6*x.^2+15*x;plot(x,y,'b','LineWidth',0.5)hold on;grid on;结果如图1：根据数学图形进行经济现象分析。

2017-12-24 16:28:38 8063

原创 Rolle中值定理的两个数学推论证明

Rolle中值定理的两个数学推论证明中值定理的两个数学推论的证明过程，体现的数学思想比较有趣，我把它备忘记录下来。Rolle中值定理的数学推论1：简单的说吧，就是，假设I区间可微、连续，如果f’(x)=0，那么f(x)=C，C为常数。可以这么理解，比如常见的常数函数f(x)=2这种常数方程，在连续可微的区间I内，永远f’(x)=0，那么f(x)一定是一个常数。证明过程还是要利用中值定理：取x1,x

2017-12-23 14:50:45 4663

原创 Matlab验算拉格朗日中值定理

拉格朗日中值定理的定义：如果函数f(x)在[a,b]上连续，且在(a,b)可导，则函数f(x)上必有一点p，使得：(f(b)-f(a) )/(b-a)=f'(p)。该定理可以认为如果函数满足拉格朗日中值定理所有条件，那么f(x)上两点连线构成的直线，与过f(p)点的直线平行。MATLAB代码：syms x y; y=x.^2; %绘出y=x*x曲线，黑色线。e1=ezplot(y,[0

2017-12-20 16:02:25 5432

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

Matlab计算微分方程曲线求导及过曲线上点的切线方程求解f(x)=x^2一元二次方程上某点的切线方程并绘制出方程的切线图。点(4,f(4))是曲线方程f(x)上的一个点，求出该点的切线并绘制出来。画出f(x)= x^2方程曲线。对f(x)进行求导得到f’(x)=2*x。根据一般的过点(a,b)的斜切线方程求出切线方程：m为导数值。变形得到过(4,f(4))的切线方程：matlab绘出的图：红色线

2017-12-19 16:33:49 13831

原创 Machine Learning：最小二乘法数学原理及简单推导

Machine Learning：最小二乘法数学原理及简单推导假设给定一系列散列值（数据集）记为D={（x1,y1），（x2，y2），（x3，y3），，，（xn，yn）}，找到一个函数y=ax+b（也可记得f(x)=ax+b）使得f(x)函数尽可能拟合D。求解函数f(x)的方法很多种。最小二乘法寻找拟合函数f(x)的原理和思想关键：平方差之和最小，即使得Q最小。即求解最小值。因为(x1,y1)，（

2017-12-18 15:21:44 9227 1

原创机器学习之线性回归的最小二乘法求解

机器学习之线性回归的最小二乘法求解假设现在一个普通的一阶线性方程，y=2*x+2*t。t是随机噪音，生成的散列点(x,y)会沿直线y=2*x上下摆动。利用最小二乘法做一次简单的一阶“曲线”拟合。用matlab做数据实验：t=randn(1,101); x=[-10:0.2:10]; y=2*x+t*2; s=scatter(x,y); s.LineWidth = 0.6;

2017-12-18 09:38:09 1167

原创 Machine Learning 之Logistic回归算法中最小二乘法的Matlab曲线拟合

Machine Learning 之Logistic回归算法中最小二乘法的Matlab曲线拟合逻辑回归是机器学习（Machine Learning）中常见的机器学习算法，在处理逻辑回归（Logistic Regression）离散数据点集时，最常用的算法是最小二乘法。古代欧洲没有“平方”的叫法，“二乘”其实就是平方。逻辑回归是相对于线性回归而言，线性回归可以较好拟合连续值。但是现实世界中的数据样本

2017-12-15 19:14:54 7162 2

原创 matlab三维山峰/山脉/山地曲面数据图

matlab三维山峰/山脉/山地曲面数据图可以使用peaks函数。比如直接peaks(80)：peaks(80) z = 3*(1-x).^2.*exp(-(x.^2) - (y+1).^2) ... - 10*(x/5 - x.^3 - y.^5).*exp(-x.^2-y.^2) ... - 1/3*exp(-(x+1).^2 - y.^2) peaks(80)的三维数据

2017-12-15 14:22:32 21818 1

原创 AI神经网络激活函数sigmoid及matlab的sigmf

AI神经网络激活函数sigmoid及matlab的sigmf神经网络中引入激活函数sigmoid作用是逻辑回归（logistic regression），引入非线性化。数学中的标准sigmoid输出范围是(0,1)。sigmoid的数学定义：在matlab中，对于sigmoid的定义实现是sigmf，但是sigmf包含多个参数：用MATLAB跑出不同的sigmoid函数曲线：x1=-10:0.1:

2017-12-15 10:52:22 9491

原创 Matlab数值分析K-means聚类算法

我写的附录文章1简单介绍了K-means聚类算法。Matlab提供了专用函数kmeans用于聚类的质心。假设随机生成两维样本数据，然后用kmeans算出聚类并标记出质心：r1=randn(5,2)-2r2=randn(5,2)+2X=[r1;r2];opts = statset('Display','final');[idx,C] = kmeans(X,2,'Replicates',2,

2017-12-14 20:43:42 2397

原创数值分析Matlab二维正态(高斯)分布以及协方差矩阵

数值分析Matlab二维正态(高斯)分布以及协方差矩阵主要是使用了matlab的mvnrnd产生随机的正态（高斯）分布二维矩阵，然后绘制出来。代码运行结果生成的正态分布实验数据如图：MATLAB代码：mu1 = [0 0];sigma1 = [4 2 ; 2 4];r1 = mvnrnd(mu1,sigma1,100);scatter(r1(:,1),r1(:,2),'r.');hold

2017-12-14 15:56:30 27172

原创数值分析Matlab绘制三维数据曲面图

主要使用了matlab提供的meshgrid和surf。一个细节是matlab中的 ./ 如果是单数，是一般的数学除法，如果是多维矩阵，则是多维矩阵对应元素的除法。matlab代码：xi=-20:0.3:20;yi=-10:0.3:10;[x,y]=meshgrid(xi,yi);z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2);surf(x,y,z);代码运

2017-12-14 09:57:22 10008 1

原创数值分析Matlab三维数据网格图

主要使用了matlab的meshgrid和mesh网格绘制函数。matlab代码：xi=-20:0.5:20;yi=-20:0.5:20;[x,y]=meshgrid(xi,yi);z=(-x.^2-y.^2);mesh(x,y,z);代码运行结果如图：

2017-12-13 20:42:50 7342

原创数值分析基础工具使用Matlab绘制双曲线

绘制结果如图所示：具体的绘制要求和实现的细节，matlab代码注释中已有：x=-3:0.05:3; %刻度为0.05单位。显示区间为[-3,3]。y1=x.^2+5; %即y1=-2x*x+5的曲线。plot(x,y1,'g','LineWidth',0.5); %绿色的曲线。线宽0.5。hold on;y2=-x.^2-5; %即y2=-2x*x-5的曲线。plot(x,y2,'b'

2017-12-13 19:09:23 11558 3

原创人工智能AI常见的经典K-means聚类算法原理和工作过程

人工智能AI常见的经典K-means聚类算法原理和工作过程K-means聚类算法亦称K聚类均值算法，K-means算法是硬聚类算法中的一种。聚类算法是一类无监督机器学习。K-means算法是计算数据聚集的算法，主要通过不断地取离种子点最近均值的算法。K-means算法是典型的基于距离的聚类算法，以距离作为相似性的评价标准，认为两个对象的距离越近，其相似度就越大。物以类聚，人以群分。K均值算法认为簇

2017-12-12 16:55:12 7438 3

原创人工智能中卷积神经网络基本原理综述

人工智能Artificial Intelligence中卷积神经网络Convolutional Neural Network基本原理综述人工智能（Artificial Intelligence，简称AI）的Deep Learning（深度学习）通过机器学习，把某一层的输出output当做下一层的输入input。在人工智能中，认为output是机器通过深度学习获得的某种“智慧”。深度学习（Deep

2017-12-11 16:10:00 32670 11

原创 simple Android EventBus code

import android.os.Bundle;import android.support.v7.app.AppCompatActivity;import android.util.Log;import org.greenrobot.eventbus.EventBus;import org.greenrobot.eventbus.Subscribe;import org.green

2017-12-08 15:39:04 571

原创 Android piles of RuntimePermissions requests, code improving

Android piles of RuntimePermissions requests, code improvingfor example:AndroidManifest.xml: java:import android.Manifest;import android.content.pm.PackageManager;import android.os.Build;im

2017-12-07 19:32:31 464

原创 Improving Android ORMLite foreign key mapping , based on new version 5.0

Improving Android ORMLite foreign key mapping , based on new version 5.0AClass.java：import com.j256.ormlite.dao.ForeignCollection;import com.j256.ormlite.field.DataType;import com.j256.ormlite.field

2017-12-01 10:10:31 626