| 查看: 1706 | 回复: 7 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
matlab循环计算转为矩阵 30金币求助
|
||
|
程序是用来做菲涅尔积分的，由于程序运算时间太长，因此想把最耗时的两重for循环运算部分改为矩阵运算形式，但没找到合适的方法，循环中运行的语句本身就是一个矩阵运算，想请教各位怎么改成完全矩阵运算的程序。 X1=-22e-3; X2=22e-3; Y1=-22e-3; Y2=22e-3; Num0=128; [x,y]=meshgrid(X1  X2-X1)/(Num0-1):X2,Y1Y2-Y1)/(Num0-1):Y2);for i=1:Num0 x1=x(1,i); for j=1:Num0 y1=y(j,1); Fij=Fun2.*exp(sqrt((x1-x).^2+(y1-y).^2)); %Fun2是一个矩阵函数 Fun3(i,j)=sum(sum(Fij)); end end |
回复此楼» 猜你喜欢
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有165人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有18人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab如何实现矩阵多行数据的循环调用
已经有5人回复
- 如何用Matlab画矩阵的散点图？
已经有3人回复
- 菜鸟求助~matlab怎么样让一个1xn的向量中的每个元素是个矩阵？
已经有11人回复
- matlab 有限元计算扩散问题，建立整体矩阵好慢。大家帮忙看看代码
已经有8人回复
- matlab中系数是矩阵的积分
已经有5人回复
- 请教关于matlab矩阵转换问题
已经有6人回复
- 求助用matlab表示一个矩阵
已经有4人回复
- 关于计算机内存和matlab产生矩阵能力的关系的困惑
已经有3人回复
- matlab 不规则矩阵 求平均值
已经有6人回复
- matlab 复数矩阵转为实数矩阵 等值线作图求助 详见帖内
已经有1人回复
- 【求助】求Matlab循环作图
已经有12人回复
- 【求助】matlab 解矩阵方程 重金悬赏解决方案
已经有10人回复
3楼2012-11-24 23:02:47
meatball1982
铜虫 (小有名气)
- 应助: 8 (幼儿园)
- 金币: -3.4
- 红花: 8
- 帖子: 297
- 在线: 108.3小时
- 虫号: 1530793
- 注册: 2011-12-09
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
csgt0: 金币+3, 辛苦啦！ 2012-11-23 09:45:47
csgt0: 金币+30, lz奖励金币 2012-12-17 09:15:17
感谢参与，应助指数 +1
csgt0: 金币+3, 辛苦啦！ 2012-11-23 09:45:47
csgt0: 金币+30, lz奖励金币 2012-12-17 09:15:17
你好。关于两个循环，实在是想不出什么好的方法，主要是求各个点到其它点的距离这一项，试了总是不行。 但如果你确定是用这一组X1，X2，Y1,Y2，就可以偷下懒。 我的程序暂时以Num0=8为例进行说明。只要X关于点对称，Y关于点对称，就好办。 你的F3是一个大的对称的矩阵。我们只要生成这个矩阵的1/8.通过转置，翻转等操作（matlab对矩阵的运算还是挺快的。），可以生成Fun3。 当Num0=8的时候，两层循环比较好。 Elapsed time is 0.000551 seconds. Elapsed time is 0.012936 seconds. 当Num0=32的时候，两种方法差不多。 Elapsed time is 0.063168 seconds. Elapsed time is 0.021089 seconds. 当Num0=128的时候，偷懒的方法要好些，可能是Num0越大，效果越好，但应该是有上限的，这要看你的内存值（用空间换时间就是）。 Elapsed time is 8.402686 seconds. Elapsed time is 0.986670 seconds. clear all clc % X1=-22e-3; % X2=22e-3; % Y1=-22e-3; % Y2=22e-3; % different X,Y X1=-20e-3; X2=20e-3; Y1=-21e-3; Y2=21e-3; Num0=128;% u can change it for 8,32 64 and 128 tm_lu=zeros(Num0/2,Num0/2); tm_lu=zeros(Num0/2,Num0); tmp =zeros(Num0,Num0); Fun2=1; x_l=X1 X2-X1)/(Num0-1):X2;y_l=Y1 Y2-Y1)/(Num0-1):Y2;[x,y]=meshgrid(x_l,y_l); %% your way tic for i=1:Num0 x1=x_l(i); for j=1:Num0 y1=y_l(j); Fij=Fun2.*exp(sqrt((x1-x).^2+(y1-y).^2)); %Fun2是一个矩阵函数 F3(i,j)=sum(Fij(  );end end toc %% tic for i=1:Num0/2 x1=x_l(i); for j=1:i y1=y_l(j); Fij=Fun2.*exp(sqrt((x1-x).^2+(y1-y).^2)); %Fun2是一个矩阵函数 tm(i,j)=sum(Fij( );end end tm_lu=tm+tm'-tril(tm'); tm_u=[tm_lu,flipud(tm_lu)']; F4=[tm_u;flipud(tm_u)]; toc % modified on : 22-Nov-2012 19:11:35 % typed by : m |
2楼2012-11-22 19:23:09
meatball1982
铜虫 (小有名气)
- 应助: 8 (幼儿园)
- 金币: -3.4
- 红花: 8
- 帖子: 297
- 在线: 108.3小时
- 虫号: 1530793
- 注册: 2011-12-09
- 性别: GG
- 专业: 理论和计算化学
4楼2012-11-25 21:58:00
xhan5楼2012-11-26 10:33:53
