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

登录注册

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

24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

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);

回复此楼

» 猜你喜欢

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

使用拉格朗日优化，但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 的主题更新

返回列表

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

[考研] 081700，311，求调剂 +14	冬十三 2026-04-04	15/750	2026-04-06 01:57 by BruceLiu320
[考研] 085600，320分求调剂 +7	大馋小子 2026-04-01	8/400	2026-04-05 21:19 by 学员8dgXkO
[考研] 工科277分求调剂材料 +8	上了上了上哦 2026-04-05	9/450	2026-04-05 13:05 by wwytracy
[考研] 318求调剂 +11	ykyhsa 2026-04-05	13/650	2026-04-05 12:44 by aidndnjck
[考研] 考研调剂 +3	mcbbc 2026-04-04	3/150	2026-04-05 10:03 by barlinike
[考研] 290求调剂 +7	luoziheng 2026-04-04	7/350	2026-04-04 23:17 by lqwchd
[考研] 333求调剂 +9	阿科逸 2026-03-31	9/450	2026-04-04 18:25 by macy2011
[考研] 22408求调剂 354分可跨专业 +3	hannnnnnn 2026-04-04	3/150	2026-04-04 14:35 by 土木硕士招生
[考研] 292求调剂 +21	是妍子也是研子 2026-03-30	22/1100	2026-04-03 21:44 by qlm5820
[硕博家园] 求老师收留 +9	lllq123 2026-04-03	9/450	2026-04-03 13:48 by 呼吸都是减肥
[考研] 材料340分调剂 +7	夏夜晚风_long 2026-04-02	9/450	2026-04-02 21:20 by dongzh2009
[考研] 考研调剂 +12	Amber00 2026-03-31	12/600	2026-04-02 09:04 by sanrepian
[考研] 11408 321分求调剂 +3	huchun12138 2026-03-30	4/200	2026-04-01 22:48 by guanxin1001
[考研] 292求调剂 +17	木虫er12138 2026-04-01	17/850	2026-04-01 21:37 by 七度不信任
[考研] 0817化工学硕调剂 +11	努力上岸中！ 2026-03-31	11/550	2026-04-01 20:30 by 赖春艳
[考研] 省双一流重点一本大学招收调剂 +4	wwwwffffff 2026-03-31	7/350	2026-04-01 15:23 by wwwwffffff
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +5	双马尾痞老板2 2026-03-31	5/250	2026-04-01 09:04 by oooqiao
[考研] 340求调剂 +4	希望如此i 2026-03-31	4/200	2026-03-31 16:40 by 690616278
[考研] 物理学调剂 +4	小羊36 2026-03-30	4/200	2026-03-31 16:16 by lishahe
[考研] 求调剂 +8	11ggg 2026-03-30	8/400	2026-03-31 13:56 by nanaliuyun