24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

小鱼118

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 5779
散金: 51
红花: 4
帖子: 206
在线: 161.5小时
虫号: 1134404
注册: 2010-10-29
性别: MM
专业: 光学

[求助] 求助MATLAB画图，二维曲线，总是出错，求助高手

i=1;
lamda=0.0019;
l=14;
n=1.44;
l1=5;
h=l/2/n;
l2=10;
l3=45;
d1=l1+h;
d2=l2+h;
R2=300;
for f=1:0.1:1000 %热透镜焦距取值范围
M=[1,d1;0,1]*[1,0;0,1]*[1,d1;0,1]*[1-h/f,2*h-h*h/f;-1/f,1-h/f]*[1,d2;0,1]*...
[1,0;0,1]*[1,l3;0,1]*[1,0;-2/R2,1]*[1,l3;0,1]*[1,0;0,1]*[1,d2;0,1]*...
[1-h/f,2*h-h*h/f;-1/f,1-h/f];
A=M(1,1);
B=M(1,2);
C=M(2,1);
D=M(2,2);
N=[1,d2;0,1]*[1,0;0,1]*[1,l3;0,1]*[1,0;-2/R2,1]*[1,l3;0,1]*[1,0;0,1]*...
[1,d2;0,1]*[1-h/f,2*h-h*h/f;-1/f,1-h/f]*[1,d1;0,1]*[1,0;0,1]*...
[1,d1;0,1]*[1-h/f,2*h-h*h/f;-1/f,1-h/f];
a=N(1,1);
b=N(1,2);
c=N(2,1);
d=N(2,2);
if(abs((A+D)/2)<=1)
w01(i)=abs((2*lamda*B)/(pi*sqrt(4-(A+D)^2)))^(1/2);
w02(i)=abs((2*lamda*b)/(pi*sqrt(4-(a+d)^2)))^(1/2);
y(i)=f;
i=i+1;
end
end
wp=0.4;
plot(y,w01,y,w02,y,wp)
MATLAB运行结果为
??? Error using ==> plot
Vectors must be the same lengths.

回复此楼

» 猜你喜欢

求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有174人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有0人回复
英国全奖博士招聘-深度学习与量子物理已经有0人回复
间接带隙半导体有效质量求助已经有0人回复

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

matlab画图求助已经有11人回复
继续求助Matlab绘图！已经有8人回复
Matlab画图求助? (有人说, 标题一定要长------------------------->) 已经有5人回复
求助matlab三维画图中：Inner matrix dimensions must agree. 已经有4人回复
求助matlab画图问题已经有3人回复
求助：谁会用matlab画图？已经有13人回复
【求助】matlab 绘图问题已经有16人回复
【求助】MATLAB多个mat中数据画图问题已经有7人回复
【求助】matlab 遗传算法工具箱绘图命令是什么？已经有6人回复
【求助】请教Matlab高手根据积分画图问题已经有4人回复
【求助】如何采用Matlab绘制电子云密度图已经有8人回复

好好爱自己

1楼 2012-10-09 15:03:04

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

luomia

木虫 (著名写手)

应助: 20 (小学生)
金币: 1779
散金: 4298
红花: 14
沙发: 1
帖子: 1417
在线: 371.2小时
虫号: 1864668
注册: 2012-06-19
性别: GG
专业: 机械摩擦学与表面技术

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小鱼118: 金币+10, ★★★★★最佳答案, 画图成功了，谢谢高手！ 2012-10-10 09:00:32
dbb627: 金币+2, 谢谢应助！ 2012-10-10 17:06:07

我改了一下你原来的程序，w01和w02以及y都给了初始长度10000，循环结束后取出其中的有效长度绘图。你运行一下试试，另外C和c的值在程序中没有使用。
i=1;
lamda=0.0019;
l=14;
n=1.44;
l1=5;
h=l/2/n;
l2=10;
l3=45;
d1=l1+h;
d2=l2+h;
R2=300;
w01=zeros(1,10000);
w02=zeros(1,10000);
y=zeros(1,10000);
for f=1:0.1:1000 %热透镜焦距取值范围
M=[1,d1;0,1]*[1,0;0,1]*[1,d1;0,1]*[1-h/f,2*h-h*h/f;-1/f,1-h/f]*[1,d2;0,1]*...
[1,0;0,1]*[1,l3;0,1]*[1,0;-2/R2,1]*[1,l3;0,1]*[1,0;0,1]*[1,d2;0,1]*...
[1-h/f,2*h-h*h/f;-1/f,1-h/f];
A=M(1,1);
B=M(1,2);

D=M(2,2);
N=[1,d2;0,1]*[1,0;0,1]*[1,l3;0,1]*[1,0;-2/R2,1]*[1,l3;0,1]*[1,0;0,1]*...
[1,d2;0,1]*[1-h/f,2*h-h*h/f;-1/f,1-h/f]*[1,d1;0,1]*[1,0;0,1]*...
[1,d1;0,1]*[1-h/f,2*h-h*h/f;-1/f,1-h/f];
a=N(1,1);
b=N(1,2);

d=N(2,2);
if abs((A+D)/2)<=1
w01(i)=abs((2*lamda*B)/(pi*sqrt(4-(A+D)^2)))^(1/2);
w02(i)=abs((2*lamda*b)/(pi*sqrt(4-(a+d)^2)))^(1/2);
y(i)=f;
i=i+1;
end
end
w01n=w01(1:i-1);
w02n=w02(1:i-1);
yn=y(1:i-1);
wp=0.4*ones(1,i-1);
plot(yn,w01n,yn,w02n,yn,wp)