计算物理基础（彭芳麟）

第二章课后习题3、习题4

习题3 科赫曲线

clc
clear all
close all

new = [0,1i];
subplot(2,2,1);
plot(new)
axis equal
for k = 1:3
old = new;
n = length(old)-1; %线段总数比节点数少1
diff = (old(2:n+1)-old(1:n))/3; %将每条线段三等分
new(1:4:4*n-3) = old(1:n); %起点坐标
new(2:4:4*n-2) = old(1:n)+diff; %插入点P1坐标
new(3:4:4*n-1) = new(2:4:4*n-2)+diff*((1+sqrt(3)*1i)/2); %插入点P2坐标
new(4:4:4*n) = old(1:n)+2*diff; %插入点P3坐标
new(4*n+1) = [1i]; %终点坐标
subplot(2,2,k+1);
plot(new)
axis equal
end

习题4 Minkowski曲线

clc
clear all
close all
new=[0,1];
subplot(2,2,1);
axis([0,1,-0.5,0.5]);
axis equal
for k = 1:4
old = new;
n = length(old)-1; %线段总数比节点数少1
diff = (old(2:n+1)-old(1:n))/4; %将每条线段四等分
p1 = old(1:n)+diff;
p2 = p1+diff*(i);
p3 = p2+diff;
p4 = p3-diff*2*(i);
p = p3-diff*(i);
p5 = p4+diff;
p6 = p5+diff*(i);
new(1:15:15*n-14) = old(1:n);
new(2:15:15*n-13) = p1;
new(3:15:15*n-12) = p1;
new(4:15:15*n-11) = p2;
new(5:15:15*n-10) = p2;
new(6:15:15*n-9) = p3;
new(7:15:15*n-8) = p;
new(8:15:15*n-7) = p;
new(9:15:15*n-6) = p4;
new(10:15:15*n-5) = p4;
new(11:15:15*n-4) = p5;
new(12:15:15*n-3) = p5;
new(13:15:15*n-2) = p6;
new(14:15:15*n-1) = p6;
new(15:15:15*n) = old(2:end);
subplot(2,2,k)
plot(new)
axis([0,1,-0.5,0.5]);
end

计算物理基础（彭芳麟）相关推荐

快速傅里叶变换_计算物理基础：第八章-快速傅里叶变换（FFT）
参考北京师范大学的<计算物理基础> 第八章-快速傅里叶变换计算物理基础_中国大学MOOC(慕课)www.icourse163.org 1.快速傅里叶变换 1.1 离散傅里叶变换及其变换 ...
计算物理基础期中考试2020.5.9
一.填空题 1.ones(3,4)*nan 2.edit 3.第j列的第2个到第5个数组成的列向量. 4.sum(A) 5.x=1:1:3 y=2:4 [x,y]=meshgrid(x,y) 二.计算 ...
matlab 仿真光学实验报告,光学实验数值仿真的三种方法及MATLAB实现
光学实验数值仿真的三种方法及 MATLAB实现 5 结论 (1)数值模拟结果表明三种方法都能对光学实验现象进行正确地仿真,因此在课堂教学中适当应用这种仿真模拟 ,将光学实验中复杂的数学 ...
基于MATLAB的多项式数据拟合方法研究-毕业论文
摘要:本论文先介绍了多项式数据拟合的相关背景,以及对整个课题做了一个完整的认识.接下来对拟合模型,多项式数学原理进行了详细的讲解,通过对文献的阅读以及自己的知识积累对原理有了一个系统的认识.介绍多项式 ...
齐次方程 matlab,齐次弦振动方程的MATLAB解法.docx
齐次弦振动方程的MATLAB解法齐次弦振动方程的MATLAB解法[摘要]弦振动问题是一个典型的波动方程的建立与求解问题.本文通过利用MATLAB特有的方程求解与画图功能,有效地构造和求解了齐次弦振动 ...
matlab 开普勒方程,理论力学中非线性问题的 MATLAB数值解
第20卷第10期大学物理Vol.20No.10 2001年10月COLL EGE PHYSICS Oct.2001 基础物理教学现代化问题理论力学中非线性问题的MAT LAB 数值解收稿日期 ...
弹簧压缩时间 matlab,用matlab解决弹簧振子摆动与时间的关系
用matlab解决弹簧振子摆动与时间的关系用 matlab 解决弹簧振子摆动与时间的关系学院:光电信息班级:应用物理(111160102) 姓名:王梅学号:11116010224201 ...
非线性可视化（3）混沌系统
承接上一篇二维相图. 如果二维相平面中出现了交叉的轨线,则说明这个系统的维度很可能大于二维. 下面就以几个经典的系统作为示范.本章不涉及太多知识点,以展示为主.主要介绍三个经典的非线性混沌系统. 1 ...
matlab标量场作图
mooc计算物理基础学习笔记二维图立体等值线 z=peaks; %下载数据 contour3(peaks) %立体等值线 matlab 中的 peaks 函数是一个典型的多元函数.它本质上是一个二 ...

计算物理基础（彭芳麟）

计算物理基础（彭芳麟）相关推荐

最新文章

热门文章