| 查看: 1401 | 回复: 14 | ||
[求助]
急求各位大神，看下我这个编程有什么问题，万分感谢 已有1人参与
|
||
这段程序作出了两个三维图形，是k和alpha关于beta和delta的。。。 clc;clear c=1;q=10;p=5.5; delta=0.10:0.01:1; beta=0.55:0.005:1; for i=1:length(delta) for j=1:length(beta) t=(1-delta(i)+beta(j)*delta(i))/beta(j)/delta(i)/q; m=beta(j)*delta(i)*q; u=m*c*t-(1-delta(i))*(p-c); K(i,j)=1/3/m*3^(1/3)*(u^2*(9*m+3^(1/2)*(u^2+27*m^2)^(1/2)))^(1/3)-1/3*u^2/m*3^(2/3)/(u^2*(9*m+3^(1/2)*(u^2+27*m^2)^(1/2)))^(1/3) alpha(i,j)=((K(i,j)*m+u)/2/K(i,j)/t/m; test(i,j)=alpha(i,j)-beta(j)*delta(i); end end mesh(delta,beta,K) xlabel('\delta') ylabel('\beta') zlabel('K') figure mesh(delta,beta,alpha) xlabel('\delta') ylabel('\beta') zlabel('\alpha') |
回复此楼» 猜你喜欢
291求调剂
已经有5人回复
-
药学专硕调剂
已经有10人回复
-
一志愿211，0703化学305分求调剂
已经有32人回复
-
求调剂
已经有9人回复
-
211本科材料化工求调剂
已经有13人回复
-
286求调剂
已经有23人回复
-
085410 273求调剂
已经有7人回复
-
291分调剂
已经有6人回复
-
310求调剂
已经有8人回复
-
药学305求调剂
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求哪位大神帮忙画一个图，用chemdraw画一个球棒图，万分感谢！
已经有3人回复
- 分离到一个化合物，打了谱，解不出，请大神帮忙，万分感谢！
已经有9人回复
- 炭石墨复合材料浸渍制备工艺，求大神指点！！万分感谢！
已经有11人回复
- 求大神赐教模拟陶瓷透波率的方法
已经有3人回复
- 求大神帮你拟合一组数据~万分感激。。。
已经有6人回复
- 哪位大神能告诉给我发一张天然气集输站场平面布置图，万分感谢！！
已经有6人回复
- 数学的大神们，求解一个数学问题，大神们帮帮忙，小女子感激万分
已经有2人回复
- 哪位大神帮忙翻译一下啊，关于混凝土自修复的，万分感激........
已经有5人回复
- 跪求大神帮看matlab小程序
已经有3人回复
- 中科院教材竟出现低级错误，求数学大神解答，万分感谢
已经有22人回复
- 【求大神帮忙】功能陶瓷类论文投稿 求大神推荐审稿超级快的期刊 万分感谢！
已经有9人回复
- 各位大神们有没有 PVC线缆喷霜试验的测试方法阿 ？？万分感谢
已经有5人回复
- 有中今年的BMVC的大神么，帮忙带一下poster，万分感谢
已经有15人回复
- excel单变量求解自动运算问题
已经有6人回复
- 紧急求助，利用Matlab对实验数据进行拟合求解参数。
已经有27人回复
charmantang
禁虫 (正式写手)
★ ★ ★ ★ ★
感谢参与，应助指数 +1
燕尾蝶12: 金币+5, ★★★★★最佳答案, 谢谢，很热心的回答 2014-09-04 15:54:32
感谢参与，应助指数 +1
燕尾蝶12: 金币+5, ★★★★★最佳答案, 谢谢，很热心的回答 2014-09-04 15:54:32
|
本帖内容被屏蔽 |
3楼2014-09-02 20:02:51
charmantang
禁虫 (正式写手)
|
本帖内容被屏蔽 |
5楼2014-09-02 21:21:59
多写了一半括号，但是如何使k和alpha图像只显示它们在[0,1]区间的部分呢？希望大神能解答。。。 
2楼2014-09-02 16:25:56
4楼2014-09-02 21:13:21
6楼2014-09-02 21:39:53
|
不知道我表达的是否清楚，其实我最终的目的是除了做出三维图，还有就是判断test(i,j)的正负，您能再帮我看下下面的程序吗？ 在这个程序中test(i,j)正负不定，所以才想着要使k和alpha在[0,1]范围内，这样可以筛掉不满足条件的delta(i)和beta(j)。不知道这种思路对不对？不知道您有更好的办法吗？万分感谢。。。。 clc;clear c=1;q=10;p=5.5; delta=0.10:0.01:1; beta=0.55:0.005:1; for i=1:length(delta) for j=1:length(beta) t=(1-delta(i)+beta(j)*delta(i))/beta(j)/delta(i)/q; m=beta(j)*delta(i)*q; u=m*c*t-(1-delta(i))*(p-c); K(i,j)=1/3/m*3^(1/3)*(u^2*(9*m+3^(1/2)*(u^2+27*m^2)^(1/2)))^(1/3)-1/3*u^2/m*3^(2/3)/(u^2*(9*m+3^(1/2)*(u^2+27*m^2)^(1/2)))^(1/3); alpha(i,j)=(K(i,j)*m+u)/2/K(i,j)/t/m; test(i,j)=alpha(i,j)-delta(i); end end mesh(delta,beta,K) xlabel('\delta') ylabel('\beta') zlabel('K') set(gca,'zlim',[0,1]) figure mesh(delta,beta,alpha) xlabel('\delta') ylabel('\beta') zlabel('\alpha') set(gca,'zlim',[0,1]) 再次感谢。。。。 |
7楼2014-09-02 22:31:20
hytao2012
铁杆木虫 (正式写手)
木头虫子
- 应助: 53 (初中生)
- 金币: 6345.8
- 散金: 115
- 红花: 11
- 帖子: 479
- 在线: 206.6小时
- 虫号: 2050091
- 注册: 2012-10-08
- 性别: GG
- 专业: 固体力学
8楼2014-09-03 11:33:26
9楼2014-09-03 12:40:15
charmantang
禁虫 (正式写手)
|
本帖内容被屏蔽 |
10楼2014-09-03 14:30:01
