24小时热门版块排行榜    

查看: 1208  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

lili870113

金虫 (小有名气)

[求助] 急求

为啥我的程序这样算出来是一个具体的数字而不是矩阵
r0=10.18*10^(-6);
h=17.8*10^(-6);
l = 100;
right = 10*10^(-6);
r = linspace(0,right,l);
zz = -h*r/(r0*sqrt(r0^2-r.^2));
theta = 180/pi.*(atan(-zz));

subplot(1,3,2)
plot(rr,theta(1:col));
xlabel('\it r \rm (μm)','fontname','Times New Roman');
ylabel('\it θ(r) \rm (°)','fontname','Times New Roman');
计算出来的时候zz等于一个具体的数字,theta总提示索引超出矩阵维度
回复此楼
开心快乐每一天
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csp101

金虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lili870113(月只蓝代发): 金币+20, 求助者申请代发金币 2018-08-13 13:38:22
月只蓝: 应助指数+1 2018-08-13 13:38:33
引用回帖:
7楼: Originally posted by lili870113 at 2018-08-11 10:17:17
我定义了,但是还是不对...

right = 10*10^(-6);
l = 100;
r = linspace(0,right,l);
r0=10.18*10^(-6)*ones(1,length(r));
h=17.8*10^(-6)*ones(1,length(r));
zz = -h.*r./(r0.*sqrt(r0.^2-r.^2));
theta = 180/pi.*(atan(-zz));
subplot(1,3,2)
plot(zz,theta);
xlabel('\it r \rm (μm)','fontname','Times New Roman');
ylabel('\it θ(r) \rm (°)','fontname','Times New Roman');
手机回复不方便,特意上了电脑,除了定义矩阵以外,你没有用点乘,所以才导致的计算得到的只是一个数
8楼2018-08-11 10:29:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

csp101

金虫 (著名写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lili870113(月只蓝代发): 金币+20, 求助者申请代发金币 2018-08-13 13:38:06
因为你的根号里面r0是一个数啊,你应该用ones(1,length(r))*r0来定义一个矩阵,其内部元素均为r0,这样才可以和r进行矩阵运算,另外,你作图时候的rr没有定义吧

发自小木虫IOS客户端
2楼2018-08-11 10:02:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lili870113

金虫 (小有名气)

引用回帖:
2楼: Originally posted by csp101 at 2018-08-11 10:02:39
因为你的根号里面r0是一个数啊,你应该用ones(1,length(r))*r0来定义一个矩阵,其内部元素均为r0,这样才可以和r进行矩阵运算,另外,你作图时候的rr没有定义吧

您好,请问一下这个ones(1,length(r))*r0是用来定义r0的吗?可以麻烦您说详细点吗?初学者,望见谅、
开心快乐每一天
3楼2018-08-11 10:11:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csp101

金虫 (著名写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lili870113(月只蓝代发): 金币+10, 求助者申请代发金币 2018-08-13 13:38:45
引用回帖:
3楼: Originally posted by lili870113 at 2018-08-11 10:11:30
您好,请问一下这个ones(1,length(r))*r0是用来定义r0的吗?可以麻烦您说详细点吗?初学者,望见谅、...

是的,定义r0,因为你要用r0和r进行运算,但是r是矩阵,而你原程序中r0是一个数,所以计算的时候MATLAB默认r0只与r中的一个元素(应该是第一个)进行计算了,自然得不到矩阵

发自小木虫IOS客户端
4楼2018-08-11 10:14:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:29 by ZPa0EcMwuECS
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:22 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:18 by ZPa0EcMwuECS
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:03 by ZPa0EcMwuECS
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-30 00:59 by winsaint
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 00:52 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 12/600 2026-08-30 00:31 by ZPa0EcMwuECS
[基金申请] 中青基了要发朋友圈吗? +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[基金申请] 国社科申报系统有变化 +3 kynobel 2026-08-23 3/150 2026-08-29 21:57 by 余韵清
[基金申请] 面上意见出来了 +8 黄鸟于飞Chao 2026-08-29 16/800 2026-08-29 20:45 by 黄鸟于飞Chao
[基金申请] 有没有仍没收到信息的 +5 德尚中行 2026-08-27 6/300 2026-08-29 18:10 by manplx
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见