| 查看: 1706 | 回复: 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') |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求哪位大神帮忙画一个图，用chemdraw画一个球棒图，万分感谢！
已经有3人回复
- 分离到一个化合物，打了谱，解不出，请大神帮忙，万分感谢！
已经有9人回复
- 炭石墨复合材料浸渍制备工艺，求大神指点！！万分感谢！
已经有11人回复
- 求大神赐教模拟陶瓷透波率的方法
已经有3人回复
- 求大神帮你拟合一组数据~万分感激。。。
已经有6人回复
- 哪位大神能告诉给我发一张天然气集输站场平面布置图，万分感谢！！
已经有6人回复
- 数学的大神们，求解一个数学问题，大神们帮帮忙，小女子感激万分
已经有2人回复
- 哪位大神帮忙翻译一下啊，关于混凝土自修复的，万分感激........
已经有5人回复
- 跪求大神帮看matlab小程序
已经有3人回复
- 中科院教材竟出现低级错误，求数学大神解答，万分感谢
已经有22人回复
- 【求大神帮忙】功能陶瓷类论文投稿 求大神推荐审稿超级快的期刊 万分感谢！
已经有9人回复
- 各位大神们有没有 PVC线缆喷霜试验的测试方法阿 ？？万分感谢
已经有5人回复
- 有中今年的BMVC的大神么，帮忙带一下poster，万分感谢
已经有15人回复
- excel单变量求解自动运算问题
已经有6人回复
- 紧急求助，利用Matlab对实验数据进行拟合求解参数。
已经有27人回复
|
不知道我表达的是否清楚，其实我最终的目的是除了做出三维图，还有就是判断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
多写了一半括号，但是如何使k和alpha图像只显示它们在[0,1]区间的部分呢？希望大神能解答。。。 2楼2014-09-02 16:25:56
4楼2014-09-02 21:13:21
charmantang
禁虫 (正式写手)
|
本帖内容被屏蔽 |
5楼2014-09-02 21:21:59
20