版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

>论坛更新日志 (3301)
>考研 (229)
>导师招生 (214)
>虫友互识 (183)
>文献求助 (106)
>休闲灌水 (86)
>考博 (58)
>硕博家园 (57)
>公派出国 (43)
>论文道贺祈福 (39)
>论文投稿 (38)
>博后之家 (27)
>教师之家 (23)
>基金申请 (20)
>找工作 (19)
>SciFinder/Reaxys (12)

返回列表

zzucbx

新虫 (初入文坛)

应助: 2 (幼儿园)
金币: 165
散金: 60
帖子: 23
在线: 9.4小时
虫号: 1814752
注册: 2012-05-13
专业: 光学

[求助] MATLAB 程序优化

这里一个程序，算了两天了，还是没画出来图形，高手能不能优化一下，让我得到图形，万分感谢！
format long
warning off
f=zeros(9,9);
m=-4:1:4;
n=-4:1:4;
a=2.5;
b=2.5;
c=0.5;
d=20;
l=60;
e=(d/c)^2;
o=c/d;
for p=1:9;
for q=1:9;
f(p,q)=(-(a^3*b)/4-(a*b^3)/4+2*a*b*l^2-3*a^3*b*m(p)^2-3*a*b^3*n(q)^2)/(2*l^4);
end
end;
syms x y s k;
z=e*f(1,1)*(heaviside(x+4/4+o)-heaviside(x+4/4-o)).*(heaviside(y+4/4+o)-heaviside(y+4/4-o))+...
e*f(1,2)*(heaviside(x+4/4+o)-heaviside(x+4/4-o)).*(heaviside(y+3/4+o)-heaviside(y+3/4-o))+...
e*f(1,3)*(heaviside(x+4/4+o)-heaviside(x+4/4-o)).*(heaviside(y+2/4+o)-heaviside(y+2/4-o))+...
e*f(1,4)*(heaviside(x+4/4+o)-heaviside(x+4/4-o)).*(heaviside(y+1/4+o)-heaviside(y+1/4-o))+...
e*f(1,5)*(heaviside(x+4/4+o)-heaviside(x+4/4-o)).*(heaviside(y+0/4+o)-heaviside(y+0/4-o))+...
e*f(1,6)*(heaviside(x+4/4+o)-heaviside(x+4/4-o)).*(heaviside(y-1/4+o)-heaviside(y-1/4-o))+...
e*f(1,7)*(heaviside(x+4/4+o)-heaviside(x+4/4-o)).*(heaviside(y-2/4+o)-heaviside(y-2/4-o))+...
e*f(1,8)*(heaviside(x+4/4+o)-heaviside(x+4/4-o)).*(heaviside(y-3/4+o)-heaviside(y-3/4-o))+...
e*f(1,9)*(heaviside(x+4/4+o)-heaviside(x+4/4-o)).*(heaviside(y-4/4+o)-heaviside(y-4/4-o))+...
e*f(2,1)*(heaviside(x+3/4+o)-heaviside(x+3/4-o)).*(heaviside(y+4/4+o)-heaviside(y+4/4-o))+...
e*f(2,2)*(heaviside(x+3/4+o)-heaviside(x+3/4-o)).*(heaviside(y+3/4+o)-heaviside(y+3/4-o))+...
e*f(2,3)*(heaviside(x+3/4+o)-heaviside(x+3/4-o)).*(heaviside(y+2/4+o)-heaviside(y+2/4-o))+...
e*f(2,4)*(heaviside(x+3/4+o)-heaviside(x+3/4-o)).*(heaviside(y+1/4+o)-heaviside(y+1/4-o))+...
e*f(2,5)*(heaviside(x+3/4+o)-heaviside(x+3/4-o)).*(heaviside(y+0/4+o)-heaviside(y+0/4-o))+...
e*f(2,6)*(heaviside(x+3/4+o)-heaviside(x+3/4-o)).*(heaviside(y-1/4+o)-heaviside(y-1/4-o))+...
e*f(2,7)*(heaviside(x+3/4+o)-heaviside(x+3/4-o)).*(heaviside(y-2/4+o)-heaviside(y-2/4-o))+...
e*f(2,8)*(heaviside(x+3/4+o)-heaviside(x+3/4-o)).*(heaviside(y-3/4+o)-heaviside(y-3/4-o))+...
e*f(2,9)*(heaviside(x+3/4+o)-heaviside(x+3/4-o)).*(heaviside(y-4/4+o)-heaviside(y-4/4-o))+...
e*f(3,1)*(heaviside(x+2/4+o)-heaviside(x+2/4-o)).*(heaviside(y+4/4+o)-heaviside(y+4/4-o))+...
e*f(3,2)*(heaviside(x+2/4+o)-heaviside(x+2/4-o)).*(heaviside(y+3/4+o)-heaviside(y+3/4-o))+...
e*f(3,3)*(heaviside(x+2/4+o)-heaviside(x+2/4-o)).*(heaviside(y+2/4+o)-heaviside(y+2/4-o))+...
e*f(3,4)*(heaviside(x+2/4+o)-heaviside(x+2/4-o)).*(heaviside(y+1/4+o)-heaviside(y+1/4-o))+...
e*f(3,5)*(heaviside(x+2/4+o)-heaviside(x+2/4-o)).*(heaviside(y+0/4+o)-heaviside(y+0/4-o))+...
e*f(3,6)*(heaviside(x+2/4+o)-heaviside(x+2/4-o)).*(heaviside(y-1/4+o)-heaviside(y-1/4-o))+...
e*f(3,7)*(heaviside(x+2/4+o)-heaviside(x+2/4-o)).*(heaviside(y-2/4+o)-heaviside(y-2/4-o))+...
e*f(3,8)*(heaviside(x+2/4+o)-heaviside(x+2/4-o)).*(heaviside(y-3/4+o)-heaviside(y-3/4-o))+...
e*f(3,9)*(heaviside(x+2/4+o)-heaviside(x+2/4-o)).*(heaviside(y-4/4+o)-heaviside(y-4/4-o))+...
e*f(4,1)*(heaviside(x+1/4+o)-heaviside(x+1/4-o)).*(heaviside(y+4/4+o)-heaviside(y+4/4-o))+...
e*f(4,2)*(heaviside(x+1/4+o)-heaviside(x+1/4-o)).*(heaviside(y+3/4+o)-heaviside(y+3/4-o))+...
e*f(4,3)*(heaviside(x+1/4+o)-heaviside(x+1/4-o)).*(heaviside(y+2/4+o)-heaviside(y+2/4-o))+...
e*f(4,4)*(heaviside(x+1/4+o)-heaviside(x+1/4-o)).*(heaviside(y+1/4+o)-heaviside(y+1/4-o))+...
e*f(4,5)*(heaviside(x+1/4+o)-heaviside(x+1/4-o)).*(heaviside(y+0/4+o)-heaviside(y+0/4-o))+...
e*f(4,6)*(heaviside(x+1/4+o)-heaviside(x+1/4-o)).*(heaviside(y-1/4+o)-heaviside(y-1/4-o))+...
e*f(4,7)*(heaviside(x+1/4+o)-heaviside(x+1/4-o)).*(heaviside(y-2/4+o)-heaviside(y-2/4-o))+...
e*f(4,8)*(heaviside(x+1/4+o)-heaviside(x+1/4-o)).*(heaviside(y-3/4+o)-heaviside(y-3/4-o))+...
e*f(4,9)*(heaviside(x+1/4+o)-heaviside(x+1/4-o)).*(heaviside(y-4/4+o)-heaviside(y-4/4-o))+...
e*f(5,1)*(heaviside(x+0/4+o)-heaviside(x+0/4-o)).*(heaviside(y+4/4+o)-heaviside(y+4/4-o))+...
e*f(5,2)*(heaviside(x+0/4+o)-heaviside(x+0/4-o)).*(heaviside(y+3/4+o)-heaviside(y+3/4-o))+...
e*f(5,3)*(heaviside(x+0/4+o)-heaviside(x+0/4-o)).*(heaviside(y+2/4+o)-heaviside(y+2/4-o))+...
e*f(5,4)*(heaviside(x+0/4+o)-heaviside(x+0/4-o)).*(heaviside(y+1/4+o)-heaviside(y+1/4-o))+...
e*f(5,5)*(heaviside(x+0/4+o)-heaviside(x+0/4-o)).*(heaviside(y+0/4+o)-heaviside(y+0/4-o))+...
e*f(5,6)*(heaviside(x+0/4+o)-heaviside(x+0/4-o)).*(heaviside(y-1/4+o)-heaviside(y-1/4-o))+...
e*f(5,7)*(heaviside(x+0/4+o)-heaviside(x+0/4-o)).*(heaviside(y-2/4+o)-heaviside(y-2/4-o))+...
e*f(5,8)*(heaviside(x+0/4+o)-heaviside(x+0/4-o)).*(heaviside(y-3/4+o)-heaviside(y-3/4-o))+...
e*f(5,9)*(heaviside(x+0/4+o)-heaviside(x+0/4-o)).*(heaviside(y-4/4+o)-heaviside(y-4/4-o))+...
e*f(6,1)*(heaviside(x-1/4+o)-heaviside(x-1/4-o)).*(heaviside(y+4/4+o)-heaviside(y+4/4-o))+...
e*f(6,2)*(heaviside(x-1/4+o)-heaviside(x-1/4-o)).*(heaviside(y+3/4+o)-heaviside(y+3/4-o))+...
e*f(6,3)*(heaviside(x-1/4+o)-heaviside(x-1/4-o)).*(heaviside(y+2/4+o)-heaviside(y+2/4-o))+...
e*f(6,4)*(heaviside(x-1/4+o)-heaviside(x-1/4-o)).*(heaviside(y+1/4+o)-heaviside(y+1/4-o))+...
e*f(6,5)*(heaviside(x-1/4+o)-heaviside(x-1/4-o)).*(heaviside(y+0/4+o)-heaviside(y+0/4-o))+...
e*f(6,6)*(heaviside(x-1/4+o)-heaviside(x-1/4-o)).*(heaviside(y-1/4+o)-heaviside(y-1/4-o))+...
e*f(6,7)*(heaviside(x-1/4+o)-heaviside(x-1/4-o)).*(heaviside(y-2/4+o)-heaviside(y-2/4-o))+...
e*f(6,8)*(heaviside(x-1/4+o)-heaviside(x-1/4-o)).*(heaviside(y-3/4+o)-heaviside(y-3/4-o))+...
e*f(6,9)*(heaviside(x-1/4+o)-heaviside(x-1/4-o)).*(heaviside(y-4/4+o)-heaviside(y-4/4-o))+...
e*f(7,1)*(heaviside(x-2/4+o)-heaviside(x-2/4-o)).*(heaviside(y+4/4+o)-heaviside(y+4/4-o))+...
e*f(7,2)*(heaviside(x-2/4+o)-heaviside(x-2/4-o)).*(heaviside(y+3/4+o)-heaviside(y+3/4-o))+...
e*f(7,3)*(heaviside(x-2/4+o)-heaviside(x-2/4-o)).*(heaviside(y+2/4+o)-heaviside(y+2/4-o))+...
e*f(7,4)*(heaviside(x-2/4+o)-heaviside(x-2/4-o)).*(heaviside(y+1/4+o)-heaviside(y+1/4-o))+...
e*f(7,5)*(heaviside(x-2/4+o)-heaviside(x-2/4-o)).*(heaviside(y+0/4+o)-heaviside(y+0/4-o))+...
e*f(7,6)*(heaviside(x-2/4+o)-heaviside(x-2/4-o)).*(heaviside(y-1/4+o)-heaviside(y-1/4-o))+...
e*f(7,7)*(heaviside(x-2/4+o)-heaviside(x-2/4-o)).*(heaviside(y-2/4+o)-heaviside(y-2/4-o))+...
e*f(7,8)*(heaviside(x-2/4+o)-heaviside(x-2/4-o)).*(heaviside(y-3/4+o)-heaviside(y-3/4-o))+...
e*f(7,9)*(heaviside(x-2/4+o)-heaviside(x-2/4-o)).*(heaviside(y-4/4+o)-heaviside(y-4/4-o))+...
e*f(8,1)*(heaviside(x-3/4+o)-heaviside(x-3/4-o)).*(heaviside(y+4/4+o)-heaviside(y+4/4-o))+...
e*f(8,2)*(heaviside(x-3/4+o)-heaviside(x-3/4-o)).*(heaviside(y+3/4+o)-heaviside(y+3/4-o))+...
e*f(8,3)*(heaviside(x-3/4+o)-heaviside(x-3/4-o)).*(heaviside(y+2/4+o)-heaviside(y+2/4-o))+...
e*f(8,4)*(heaviside(x-3/4+o)-heaviside(x-3/4-o)).*(heaviside(y+1/4+o)-heaviside(y+1/4-o))+...
e*f(8,5)*(heaviside(x-3/4+o)-heaviside(x-3/4-o)).*(heaviside(y+0/4+o)-heaviside(y+0/4-o))+...
e*f(8,6)*(heaviside(x-3/4+o)-heaviside(x-3/4-o)).*(heaviside(y-1/4+o)-heaviside(y-1/4-o))+...
e*f(8,7)*(heaviside(x-3/4+o)-heaviside(x-3/4-o)).*(heaviside(y-2/4+o)-heaviside(y-2/4-o))+...
e*f(8,8)*(heaviside(x-3/4+o)-heaviside(x-3/4-o)).*(heaviside(y-3/4+o)-heaviside(y-3/4-o))+...
e*f(8,9)*(heaviside(x-3/4+o)-heaviside(x-3/4-o)).*(heaviside(y-4/4+o)-heaviside(y-4/4-o))+...
e*f(9,1)*(heaviside(x-4/4+o)-heaviside(x-4/4-o)).*(heaviside(y+4/4+o)-heaviside(y+4/4-o))+...
e*f(9,2)*(heaviside(x-4/4+o)-heaviside(x-4/4-o)).*(heaviside(y+3/4+o)-heaviside(y+3/4-o))+...
e*f(9,3)*(heaviside(x-4/4+o)-heaviside(x-4/4-o)).*(heaviside(y+2/4+o)-heaviside(y+2/4-o))+...
e*f(9,4)*(heaviside(x-4/4+o)-heaviside(x-4/4-o)).*(heaviside(y+1/4+o)-heaviside(y+1/4-o))+...
e*f(9,5)*(heaviside(x-4/4+o)-heaviside(x-4/4-o)).*(heaviside(y+0/4+o)-heaviside(y+0/4-o))+...
e*f(9,6)*(heaviside(x-4/4+o)-heaviside(x-4/4-o)).*(heaviside(y-1/4+o)-heaviside(y-1/4-o))+...
e*f(9,7)*(heaviside(x-4/4+o)-heaviside(x-4/4-o)).*(heaviside(y-2/4+o)-heaviside(y-2/4-o))+...
e*f(9,8)*(heaviside(x-4/4+o)-heaviside(x-4/4-o)).*(heaviside(y-3/4+o)-heaviside(y-3/4-o))+...
e*f(9,9)*(heaviside(x-4/4+o)-heaviside(x-4/4-o)).*(heaviside(y-4/4+o)-heaviside(y-4/4-o));
h=int(z^2*cos(1i*k*x),x,-sqrt(1-y^2)+s/2,sqrt(1-y^2)-s/2);
j=subs(h,y,-sqrt(1-(s/2)^2)/100:sqrt(1-(s/2)^2)/10000:sqrt(1-(s/2)^2)/100);
g=sum(j(

)*(sqrt(1-(s/2)^2)/100);%得出是g(s,0)
t=0:0.01:2;
g=subs(g,s,t)/0.867349799305449;
plot(t,g);

回复此楼

» 猜你喜欢

283求调剂已经有11人回复
深圳大学硕士招生（2026秋，传感器方向，仅录取第一志愿）已经有8人回复
0703化学调剂已经有7人回复
344求调剂已经有3人回复
304求调剂已经有4人回复
0703化学调剂 290分有科研经历，论文在投已经有6人回复
材料专硕306英一数二已经有3人回复
环境工程调剂已经有3人回复
290求调剂已经有10人回复
321求调剂已经有4人回复

» 本主题相关价值贴推荐，对您同样有帮助:

使用拉格朗日优化，但Matlab 运行不出结果来已经有7人回复
优化 matlab编程已经有4人回复
matlab 遗传算法工具箱gaot gatbx下载已经有183人回复
求教matlab线性优化求最小值的问题已经有7人回复
求助：遗传算法的matlab编程求解资源优化问题已经有7人回复
MATLAB RBF神经网络已经有14人回复
原创：用MATLAB做的GUI做数据拟合，通过全局优化算法（含原文件，code，教程）已经有55人回复
matlab程序运行时快时慢已经有8人回复
Matlab中save问题已经有5人回复
求粒子群优化算法进行变量优选的matlab程序已经有20人回复
急求一个matlab编码，各位如果知道一定要帮我一把啊已经有10人回复
MATLAB经典算法程序已经有531人回复
求回归方程参数估计的最优化算法matlab代码已经有9人回复
求助有关Matlab有约束非线性最优化问题已经有9人回复
Matlab人工神经网络工具箱的位置及使用已经有5人回复
【求助】matlab 二次规划的优化的问题已经有4人回复
【求助】多目标pso算法的MATLAB程序已经有3人回复
【求助】智能优化算法用哪些仿真工具？已经有5人回复
【分享】推荐一本好书《精通MATLAB最优化计算》已经有880人回复

1楼 2013-09-25 21:50:20

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

jerkwin

专家顾问 (正式写手)

专家经验: +14
应助: 454 (硕士)
金币: 20699.1
散金: 148
红花: 81
帖子: 813
在线: 2648.3小时
虫号: 1023452
注册: 2010-05-19
专业: 理论和计算化学
管辖: 分子模拟

【答案】应助回帖

感谢参与，应助指数 +1

你这式子长的吓死人, 能算出来才是怪事.

赞一下

回复此楼

2楼2013-09-25 22:19:29

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

babyoflion

金虫 (著名写手)

应助: 183 (高中生)
金币: 4776.5
散金: 301
红花: 30
帖子: 2507
在线: 671.9小时
虫号: 2297725
注册: 2013-02-25
性别: GG
专业: 光学

都是高手啊

回复此楼

3楼2013-09-25 23:01:37

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zzucbx

新虫 (初入文坛)

应助: 2 (幼儿园)
金币: 165
散金: 60
帖子: 23
在线: 9.4小时
虫号: 1814752
注册: 2012-05-13
专业: 光学

引用回帖:

2楼: Originally posted by jerkwin at 2013-09-25 22:19:29
你这式子长的吓死人, 能算出来才是怪事.

原本只是一个卷积，但是在上面用conv算不出来，只能手写上去了。谢谢帮助

赞一下

回复此楼

4楼2013-09-26 08:17:09

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

jerkwin

专家顾问 (正式写手)

专家经验: +14
应助: 454 (硕士)
金币: 20699.1
散金: 148
红花: 81
帖子: 813
在线: 2648.3小时
虫号: 1023452
注册: 2010-05-19
专业: 理论和计算化学
管辖: 分子模拟

【答案】应助回帖

你不会用矩阵写, 至少也用循环写啊
总不能一个一个都写出来吧

赞一下

回复此楼

5楼2013-09-26 08:58:02

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 zzucbx 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 344求调剂 +3	knight344 2026-03-16	3/150	2026-03-16 09:42 by 无际的草原
[考研] 化学工程321分求调剂 +6	大米饭！ 2026-03-15	6/300	2026-03-16 07:58 by wang_dand
[考博] 东华理工大学化材专业26届硕士博士申请 +6	zlingli 2026-03-13	6/300	2026-03-15 20:00 by ryzcf
[考研] 309求调剂 +4	花与叶@ 2026-03-10	4/200	2026-03-14 21:26 by a不易
[考研] 304求调剂 +5	小熊joy 2026-03-14	5/250	2026-03-14 21:07 by peike
[考研] 255求调剂 +3	李嘉慧， 2026-03-12	4/200	2026-03-14 16:58 by 有只狸奴
[考研] 求调剂，药学 +3	归零lbm 2026-03-09	5/250	2026-03-14 02:21 by JourneyLucky
[考研] 求调剂 +3	清风问长安 2026-03-09	3/150	2026-03-14 02:15 by JourneyLucky
[考研] 一志愿华中农业大学071010，总分三百二,求调剂 +3	困困困困坤坤 2026-03-10	3/150	2026-03-14 00:35 by JourneyLucky
[考研] 招收0805（材料）调剂 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 材料与化工（0856）304求B区调剂 +6	邱gl 2026-03-12	7/350	2026-03-13 23:24 by 邱gl
[考研] 求材料调剂 +5	隔壁陈先生 2026-03-12	5/250	2026-03-13 22:03 by 星空星月
[考研] 四川大学085601材料工程专硕初试294求调剂 +4	祝我们好在冬天 2026-03-11	4/200	2026-03-13 21:39 by peike
[考研] 材料工程调剂 +4	咪咪空空 2026-03-11	4/200	2026-03-13 19:57 by JourneyLucky
[考研] 310求调剂 +3	【上上签】 2026-03-11	3/150	2026-03-13 16:16 by JourneyLucky
[考研] 【0856】化学工程（085602）313 分，本科学科评估A类院校化学工程与工艺，诚求调剂 +7	小刘快快上岸 2026-03-11	7/350	2026-03-13 16:06 by ruiyingmiao
[考研] 268求调剂 +4	好运连绵不绝 2026-03-12	4/200	2026-03-13 10:45 by hyswxzs
[考研] 296求调剂 +3	大口吃饭身体健 2026-03-13	3/150	2026-03-13 10:31 by 学员8dgXkO
[考研] 283求调剂，材料、化工皆可 +8	苏打水7777 2026-03-11	10/500	2026-03-13 09:06 by Linda Hu
[考研] 293求调剂，一志愿陕师大生物学 +3	??????.?.??? 2026-03-09	3/150	2026-03-11 10:02 by 学员8dgXkO