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

登录注册

应《网络安全法》要求，自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);

回复此楼

» 猜你喜欢

一志愿北交大材料工程总分358求调剂已经有6人回复
312求调剂已经有7人回复
材料334求调剂已经有19人回复
生物与医药273求调剂已经有6人回复
301求调剂已经有10人回复
专硕0854初试考材科基，求调剂已经有5人回复
0703化学已经有15人回复
0703调剂，一志愿天津大学319分已经有9人回复
085600材料与化工301分求调剂院校已经有10人回复
求调剂已经有8人回复

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

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

返回列表

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

[考研] 化学357分，考研调剂 +11	.Starry. 2026-04-04	12/600	2026-04-06 06:28 by houyaoxu
[考研] 308求调剂 +11	倘若起风了呢 2026-04-05	11/550	2026-04-05 23:21 by 来看流星雨10
[考研] 生物与医药086000调剂一志愿西北农林320分 +3	美美女士 2026-04-03	3/150	2026-04-05 21:55 by 学员8dgXkO
[考研] 求调剂，一志愿厦门大学，生物与医药，总分272，本科211 +3	Electron1cc 2026-04-01	4/200	2026-04-05 20:24 by lys0704
[考研] 277求调剂 +5	考研调剂lxh 2026-04-05	5/250	2026-04-05 19:03 by chy09050039
[考研] 一志愿南京航空航天大学，080500材料科学与工程学硕 +10	@taotao 2026-03-30	10/500	2026-04-05 17:57 by jj987
[考研] 312求调剂 +3	Say Never 2026-04-04	3/150	2026-04-05 11:19 by guoweigw
[考研] 材料与化工363求推荐 +7	zh096 2026-04-04	7/350	2026-04-05 09:11 by 陌秋26
[考研] 材料求调剂 +10	呢呢妮妮 2026-04-01	10/500	2026-04-04 23:12 by 无际的草原
[考研] 本科211 分数293请求调剂 +4	莲菜就是藕吧 2026-04-01	4/200	2026-04-04 22:32 by hemengdong
[考研] 材料工程085601数二英一335求调剂 +6	双马尾痞老板2 2026-03-31	6/300	2026-04-04 22:29 by hemengdong
[考研] 一志愿0817化学工程与技术，求调剂 +24	我不是只因 2026-04-02	28/1400	2026-04-04 15:15 by dongzh2009
[考研] 348分环境工程·调剂 +10	吴彦祖24k 2026-04-03	11/550	2026-04-04 14:19 by 无际的草原
[考研] 282求调剂 +20	ycy1201 2026-04-01	22/1100	2026-04-04 00:42 by userper
[考研] 311求调剂 +20	zchqwer 2026-04-01	22/1100	2026-04-03 22:09 by lglzsd
[考研] 320调剂 +4	农业工程与信息� 2026-04-03	4/200	2026-04-03 21:40 by lbsjt
[考研] 321求调剂 +17	y-yh 2026-04-01	20/1000	2026-04-03 12:57 by y-yh
[考研] 266求调剂 +4	学员97LZgn 2026-04-02	4/200	2026-04-02 13:03 by yulian1987
[考研] 322求调剂 +5	熹僖XX 2026-03-31	6/300	2026-04-02 10:08 by 求调剂zz
[考研] 合肥区域性重点一本招收调剂 +4	6266jl 2026-03-30	8/400	2026-03-31 18:43 by 6266jl