24小时热门版块排行榜    

查看: 537  |  回复: 2

54681811

铜虫 (初入文坛)

[求助] matlab 修改问题

??? Error using ==> rdivide
Matrix dimensions must agree.
以下是程序的代码>>
n=1.84;
l=0.8;
d=4.43e-6;
a=3.78;
r=0.15;
K=0.0523;
dnt=4.7e-6;
wp=0.02;
np=1-808/1064;
nt=0.5;
p=2:1:30;
f1=3.14*K*wp*wp/np/dnt/nt;
f2=1/(1-exp(-a*l));
ff1=f1.*f2./p;

f3=3.14*K*wp*wp/2/np/nt/d/r/(n-1);
ff2=f3.*f2./p;

f=1./(7e-2)./(1./ff1+1./ff2);
L1=2;
L2=10;
L3=15;
L4=15;
L5=50;
lambda=1.064*10^-3;
R=200;
n=1;


m=length(f);
w=zeros(1,m);
for n=1:1:m
fs=f(n);
m0=[1,2.17;0,1];
m1=[1,L1;0,1];
m2=[1,0;0,1];
m3=[1,L1;0,1];
m4=[1,2.17;0,1];
m5=[1,0;-1/82.5,1];
m6=[1,2.17;0,1];
m7=[1,L2;0,1];
m8=[1,0;-1/90,1];
m9=[1,L3;0,1];
m10=[1,3/1.6;0,1];
m11=[1,L4;0,1];
m12=[1,50/1.6;0,1];
m13=[1,L5;0,1];
m14=[1,0;-2/R,1];
m15=[1,L5;0,1];
m16=[1,50/1.6;0,1];
m17=[1,L4;0,1];
m18=[1,3/1.6;0,1];
m19=[1,L3;0,1];
m20=[1,0;-1/90,1];
m21=[1,L2;0,1];
m22=[1,2.17;0,1];
m23=[1,0;-1/82.5,1];

M=m0*m1*m2*m3*m4*m5*m6*m7*m8*m9*m10*m11*m12*m13*m14*m15*m16*m17*m18*m19*m20*m21*m22*m23;
A=M(1,1);
B=M(1,2);
C=M(2,1);
D=M(2,2);
if(abs((A+D)/2)<=1)
    w(n)=sqrt(lambda*abs(B)/(pi*sqrt(1-((A+D)/2)^2)));
    n=n+1;
end
end

z=3:1:200;
ff=3.1415926.*w.*w./1064e-6;
w1=w.*(1+(z./ff).^2).^0.5;
plot(z,w1,'r')
xlabel('谐振腔内位置(mm)')
ylabel('光斑半径(mm)')
回复此楼

» 猜你喜欢

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

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

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
54681811: 金币+5 2013-01-12 09:43:55
z跟ff不是一般大啊
showmethemoney
2楼2013-01-10 13:30:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ygjianqiang

铜虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
dbb627: 金币+1, 谢谢应助! 2013-01-10 15:37:35
54681811: 金币+6 2013-01-12 09:44:06
因为没有给出数学模型,所以也不太清楚你的应该怎么改。
你的z大小为1*198
而ff大小为1*29
因此两者是大小不一样的,造成无法在同一式中进行计算
3楼2013-01-10 15:19:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 54681811 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +16 Lanmanbaby 2026-08-09 25/1250 2026-08-11 15:51 by 医学老男孩
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +7 Tide man 2026-08-10 8/400 2026-08-11 15:49 by 铁帽子农民
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +3 工厂打螺丝 2026-08-08 3/150 2026-08-11 13:42 by 会议编辑
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] FileCode能看出啥? +7 要乐观耀哥 2026-08-10 17/850 2026-08-11 10:50 by 要乐观耀哥
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +5 布布和一二 2026-08-10 6/300 2026-08-11 10:15 by 医学老男孩
[基金申请] filecode +9 documentary 2026-08-10 9/450 2026-08-11 10:08 by wk7465
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
信息提示
请填处理意见