基于有限差分法求解的二维谐振子的MATLAB程序如下。哪位大神能帮我做个注明啊,完全看不懂啊,,急0

____丿呆呆丶2017.04.15浏览20次分享举报

tic clc clear L=20; W=20; N=20; M=20; hx=L/(2*N); hy=W/(2*M); S=zeros((2*M-1)*(2*N-1)); for m=1:2*M-1 D(m,m)=-1/(hy.^2); end for m=1:2*N-2 mx=(2*M-1)*(m-1)+1; my=(2*M-1)*(m... tic

clc

clear

L=20;

W=20;

N=20;

M=20;

hx=L/(2*N);

hy=W/(2*M);

S=zeros((2*M-1)*(2*N-1));

for m=1:2*M-1

D(m,m)=-1/(hy.^2);

end

for m=1:2*N-2

mx=(2*M-1)*(m-1)+1;

my=(2*M-1)*(m-1)+2*M-1;

nx=(2*M-1)*(m-1)+(2*M-1)+1;

ny=(2*M-1)*(m-1)+2*(2*M-1);

S(mx:my,nx:ny)=D;

S(nx:ny,mx:my)=D;

end

for n=1:2*N-1

C=zeros(2*M-1);

for m=1:2*M-1

C(m,m)=2*(1/hx^2+1/hy^2)+((n-M)*hx)^2+((m-N)*hy)^2;

end

for m=1:2*M-2

C(m,m+1)=-1/hx.^2;

C(m+1,m)=-1/hx.^2;

end

cx=(2*M-1)*(n-1)+1;

cy=(2*M-1)*(n-1)+2*M-1;

S(cx:cy,cx:cy)=C;

end

[V,E]=eig(S); Y=V(:,2); k=1;

for j=1:2*M-1;

for i=1:2*N-1;

ZZ(i,j)=Y(k);

k=k+1;

end

end

y=linspace(-W/2,W/2,2*N-1);

x=linspace(-L/2,L/2,2*M-1);

[x1,y1]=meshgrid(x,y);

surf(x1,y1,-ZZ);

shading interp;

zlabel('¦•');

toc 展开

matlab二维谐振子,基于有限差分法求解的二维谐振子的MATLAB程序如下。哪位大神能帮我做个注明啊,完全看不懂啊,,急...相关推荐

  1. matlab if语句多个执行举例,初学Matlab,有两个语句,if语句和switch语句,有两个例子哪位大神能帮我讲讲...

    问题描述: 初学Matlab,有两个语句,if语句和switch语句,有两个例子哪位大神能帮我讲讲 if logical_expression statements elseif logical_ex ...

  2. matlab结果输出的代码,哪位大神能帮我看一下下列代码输出的结果是啥!

    该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 N=6000; %信息序列的长度 K=11; %11个抽头系数 %实数信道 actual_isi=[0.05 -0.063 0.088 -0.126 -0 ...

  3. matlab联立两个方程组求解,实验二 Matlab求解数学问题(终稿)2

    实验二 MATLAB求解数学问题 2.1实验目的 掌握MATLAB在大学数学问题中的基本应用,会使用MATLAB软件求解高等数学.线性代数和概率统计中的常见问题. 2.2实验要求 掌握MATLAB简单 ...

  4. matlab梯形模糊数,基于梯形模糊数的OWA方法与matlab应用

    基于梯形模糊数的OWA 方法与matlab 应用 一.梯形模糊数定义 1.定义记(,,,),=-∞<<<< c a b d c a b d ,称为%A 为梯形模糊数,当0> ...

  5. 基于matlab的电机故障诊断,基于BP神经网络的电机故障诊断的研究(MATLAB,程序)

    基于BP神经网络的电机故障诊断的研究(MATLAB,程序)(课题申报表,任务书,开题报告,中期检查表,外文翻译,论文25600字,程序,答辩PPT) 摘 要 随着我国经济发展的不断提高,电气化的应用越 ...

  6. matlab 2010无法运行程序,matalb r2010a安装后打开出现一系列警告,无法运行,哪位大神帮...

    该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 Matlab R2010a安装完成后的错误:Warning: D:\Program Files\MATLAB\R2010b\toolbox\local\p ...

  7. 【优化求解】基于遗传算法求解电动汽车充电管理优化问题Matlab代码

    1 简介 在现有的电池技术和充电条件下,快换式充电站成为国内纯电动公交最主要的能量补给摸式.针对电池组充电电费过高和增加配电网峰谷差的问题,以及换电站内快换工位和备用电池空闲的情况,在保证车辆正常运营 ...

  8. 【路径规划】基于遗传算法求解灾情巡视路径问题附matlab代码

    1 内容介绍 灾情巡视属于旅行商问题,具有广泛的应用价值.假定有若干巡视组,分工协作对所辖区域内的各灾民聚集地进行巡视,需要对各巡视组的巡视任务,巡视路线进行合理的分配和设计.在现实生活中,各被巡视地 ...

  9. 【优化求解】基于遗传算法求解电力系统最优潮流问题matlab代码

    1 简介 遗传算法是一种通过在整个解空间多渠道同时搜索以找到全局最优解的寻优方法,针对电力系统无功优化问题,提出基于遗传算法的电力系统无功优化的求解方法,用实例验证了遗传算法对电力系统无功优化的效果, ...

最新文章

  1. php写js跳转页面,php实现页面的三种跳转
  2. Asp.Net 学习资源列表(转)
  3. v-for和v-if的问题
  4. 文东工作室开通微信公众号了!欢迎订阅!~
  5. APScheduler——定时任务框架
  6. 手把手教你玩转CSS3 3D技术
  7. python中shuffle是什么意思_选择vs.Shuffle,Python
  8. [nsis]安装包反编译
  9. CDOJ 485 UESTC 485 Game (八数码变形,映射,逆cantor展开)
  10. cf的服务器型号,常用CF法兰的型号和尺寸表
  11. css的div垂直居中的方法,CSS使DIV内容垂直居中的方法
  12. Photoshop使用十大技巧
  13. 电脑蓝牙模式接收手机文件
  14. python面向对象编程实例:烤地瓜
  15. java二进制保存图片_Java中如何把图片转换成二进制流
  16. java 封装类_java的基本类型和对应的封装类
  17. ARM架构Generic Interrupt Controller(GIC)详解之术语介绍
  18. Python---Excel文件xls格式转为xlsx格式
  19. ubc 文学院 计算机,UBC大学电气与计算机教授获加拿大国家最高科学奖
  20. 通过加速计和地磁传感器实现方位角

热门文章

  1. wex5部署教程到数据库
  2. 正则基本知识和常用正则
  3. 修改el-popover和el-select样式
  4. @ImportResource()注解的使用
  5. Java面试——Spring系列总结
  6. Module build failed: Error: Node Sass version 5.0.0 is incompatible with ^4.0.0.
  7. 后端开发 java_Java后端开发三年,你不得不了解的JVM
  8. c++如何快速写出get set_如何快速写出一篇合格的毕业论文
  9. java 反射 内存_Java 反射
  10. python中msg是什么意思_python下载.msg文件的附件