应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 急求
91/1返回列表
查看: 1209  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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总提示索引超出矩阵维度
回复此楼

» 猜你喜欢
开心快乐每一天
1楼 2018-08-11 09:56:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

lili870113

金虫 (小有名气)
引用回帖:
4楼: Originally posted by csp101 at 2018-08-11 10:14:18
是的，定义r0,因为你要用r0和r进行运算，但是r是矩阵，而你原程序中r0是一个数，所以计算的时候MATLAB默认r0只与r中的一个元素（应该是第一个）进行计算了，自然得不到矩阵
...

那能麻烦你告诉我一下具体怎么加吗？非常感谢了
赞一下 回复此楼
开心快乐每一天
5楼2018-08-11 10:15:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csp101

金虫 (著名写手)
引用回帖:
5楼: Originally posted by lili870113 at 2018-08-11 10:15:34
那能麻烦你告诉我一下具体怎么加吗？非常感谢了...

？不是告诉你了吗……把r0定义成一个矩阵

发自小木虫IOS客户端
赞一下 回复此楼
6楼2018-08-11 10:16:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lili870113

金虫 (小有名气)
引用回帖:
6楼: Originally posted by csp101 at 2018-08-11 10:16:46
？不是告诉你了吗……把r0定义成一个矩阵
...

我定义了，但是还是不对
赞一下 回复此楼
开心快乐每一天
7楼2018-08-11 10:17:17
已阅   回复此楼   关注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');
手机回复不方便，特意上了电脑，除了定义矩阵以外，你没有用点乘，所以才导致的计算得到的只是一个数
赞一下(1人) 回复此楼
8楼2018-08-11 10:29:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lili870113

金虫 (小有名气)
引用回帖:
8楼: Originally posted by csp101 at 2018-08-11 10:29:26
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));
su ...

非常感谢，
回复此楼
开心快乐每一天
9楼2018-08-11 10:39:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lili870113 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:33 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.55.1.O.54,科目全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:36 by ZPa0EcMwuECS
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:36 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:34 by ZPa0EcMwuECS
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:25 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 01:03 by ZPa0EcMwuECS
[基金申请] 中青基了要发朋友圈吗？ +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 09:26 by G6APbkg8SA6w
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 我不理解！ +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) +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
信息提示
关闭
请填处理意见
关闭