24小时热门版块排行榜

>论坛更新日志 (3766)
>虫友互识 (516)
>文献求助 (356)
>基金申请 (350)
>导师招生 (266)
>休闲灌水 (196)
>论文投稿 (173)
>硕博家园 (140)
>考博 (126)
>招聘信息布告栏 (96)
>博后之家 (69)
>攻关文献(高奖励) (43)
>找工作 (37)
>考研 (36)
>教师之家 (24)
>医学 (21)

返回列表

饶玲raoling

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 68
帖子: 32
在线: 8.1小时
虫号: 3499385
注册: 2014-10-26
专业: 原子和分子物理

[交流] MATLAB程序已有2人参与

clear;
d=2;
a=1.0115;
b=0.876;
h=b*sqrt(1-d^2/(4*a^2));
for x=0:pi/20:pi;
s=(-h*a^2*sin(x)+a*sqrt(h^2*a^2*sin(x).^2-(b^2*cos(x).^2+a^2*sin(x).^2)*(h^2-b^2)))/(b^2*cos(x).^2+a^2*sin(x).^2);
plot(s.*cos(x),s.*sin(x))
end
这个代码为什么不能输出结果，求帮助？

回复此楼

» 猜你喜欢

【急招】合肥工大核聚变材料计算方向2026级工程博士生已经有4人回复
大豆异黄酮分离已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有155人回复
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标已经有10人回复
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名已经有1人回复
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生已经有4人回复
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究已经有12人回复
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有1人回复

1楼 2015-11-05 10:39:06

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

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68145.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.8小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

CODE:

clear;

d=2;

a=1.0115;

b=0.876;

h=b*sqrt(1-d^2/(4*a^2));

x=0:pi/20:pi;

s=(-h*a^2*sin(x)+a*sqrt(h^2*a^2*sin(x).^2-(b^2*cos(x).^2+a^2*sin(x).^2)*(h^2-b^2)))/(b^2*cos(x).^2+a^2*sin(x).^2);  

plot(s.*cos(x),s.*sin(x))

附图1.png

赞一下

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

2楼2015-11-05 12:22:09

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

饶玲raoling

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 68
帖子: 32
在线: 8.1小时
虫号: 3499385
注册: 2014-10-26
专业: 原子和分子物理

引用回帖:

2楼: Originally posted by 月只蓝 at 2015-11-05 12:22:09
clear;
d=2;
a=1.0115;
b=0.876;
h=b*sqrt(1-d^2/(4*a^2));
x=0:pi/20:pi;
s=(-h*a^2*sin(x)+a*sqrt(h^2*a^2*sin(x).^2-(b^2*cos(x).^2+a^2*sin(x).^2)*(h^2-b^2)))/(b^2*cos(x).^2+a^2*sin(x).^2);
plot( ...

这个结果是不正确的，理论上应该是底边的总长度是2的

赞一下

回复此楼

3楼2015-11-05 15:11:55

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