IFS为迭代函数系统，是一种构造分形的方法。 方法的核心是以概率p的方式对初始点进行仿射平移变换。 变换公式为： x(k+1)=a*x(k) + b*y(k) + ey(k+1)=c*x(k) + d*y(k) + f 关键是确定IFS码表，只要码表确定了，后面的编程套路基本一样。 下面 ...

这里用到的还是最小二乘方法，和上一次这篇文章原理差不多。 就是首先构造最小二乘函数，然后对每一个系数计算偏导，构造矩阵乘法形式，最后解方程组。 比如有一个二次曲面：z=ax^2+by^2+cxy+dx+ey+f 首先构造最小二乘函数，然后计算系数偏导（我直接手写了）： 解方程组（下图 ...

　　Julia集是一个在复平面上形成分形的点的集合，它最早由法国数学家Gaston Julia发现。 　　Julia集合可以由下式进行反复迭代得到：f(z) = z2 + c, 其中z是复平面某一点，c是一个复常数。把这个公式反复迭代，最终会得到一个复数C,然后根据C的模的大小，把这个点映射 ...

对于函数f(x)=λsin(πx)，λ∈(0,1]，使用matlab计算随着λ逐渐增大，迭代x=f(x)的值，代码如下： function y=diedai(f,a,x1) N=32; y=zeros(N,1); for i=1:1e4 x2=f(a,x1); x1 ...

可以参考之前的B样条曲线插值的方法。 要想利用控制点进行B样条曲面插值，前提是控制点不能是类似点云一样的无序点，一定要是分布在UV坐标空间中的栅格点。 方法就是首先对U方向进行B样条曲线插值，然后利用插值的结果再在V方向上进行B样条曲线插值即可。 下面用matlab自带的Peaks函数 ...

本文整理了一些三维绘图的内容，代码都比较简单 1、surf、surfc、surfl之后可以试试shading interp 2、绘制隐函数 3、参数方程绘制曲面 4、三维到平面：等高线+引力线 5、快速绘制球面(椭球面) 6、绘制旋转面 7、绕轴旋转 8、分段曲面 1、surf、surfc ...

近一段时间一直在研究分形，写了几个分形相关的程序，这是其中一个。程序中里面包含近20种分形图形的生成算法。 （1） 科赫(Koch)雪花 （2） 列维(levy)曲线 （3） 龙形曲线（Dragon ...

插值和拟合都是数据优化的一种方法，当实验数据不够多时常常须要用到这样的方法来绘图。 在matlab中都有特定的函数来完毕这些功能。 这两种方法的确别在于： 当測量值是准确的，没有误差时，一般用插值； 当測量值与真实值有误差时。一般用数据拟合。 插值 ...