| 查看: 1102 | 回复: 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总提示索引超出矩阵维度 |
回复此楼» 猜你喜欢
材料考研调剂
已经有3人回复
-
材料调剂
已经有12人回复
-
英一数一408，总分284，二战真诚求调剂
已经有14人回复
-
085410 一志愿211 22408分数359求调剂
已经有4人回复
-
271求调剂
已经有19人回复
-
385分 生物学（071000）求调剂
已经有3人回复
-
一志愿安徽大学计算机科学与技术学硕，331分求调剂
已经有3人回复
-
318求调剂，计算材料方向
已经有8人回复
-
291求调剂
已经有25人回复
-
一志愿北京科技大学085601材料工程英一数二初试总分335求调剂
已经有6人回复
csp101
金虫 (著名写手)
- 应助: 9 (幼儿园)
- 金币: 9688.1
- 散金: 120
- 红花: 7
- 帖子: 1444
- 在线: 81.8小时
- 虫号: 3076622
- 注册: 2014-03-21
- 性别: GG
- 专业: 零件成形制造
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lili870113(月只蓝代发): 金币+20, 求助者申请代发金币 2018-08-13 13:38:06
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
lili870113
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1593.7
- 散金: 50
- 帖子: 205
- 在线: 41.5小时
- 虫号: 1815521
- 注册: 2012-05-14
- 性别: MM
- 专业: 生物力学
3楼2018-08-11 10:11:30
csp101
金虫 (著名写手)
- 应助: 9 (幼儿园)
- 金币: 9688.1
- 散金: 120
- 红花: 7
- 帖子: 1444
- 在线: 81.8小时
- 虫号: 3076622
- 注册: 2014-03-21
- 性别: GG
- 专业: 零件成形制造
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lili870113(月只蓝代发): 金币+10, 求助者申请代发金币 2018-08-13 13:38:45
lili870113(月只蓝代发): 金币+10, 求助者申请代发金币 2018-08-13 13:38:45
|
是的，定义r0,因为你要用r0和r进行运算，但是r是矩阵，而你原程序中r0是一个数，所以计算的时候MATLAB默认r0只与r中的一个元素（应该是第一个）进行计算了，自然得不到矩阵 发自小木虫IOS客户端 |
4楼2018-08-11 10:14:18
lili870113
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1593.7
- 散金: 50
- 帖子: 205
- 在线: 41.5小时
- 虫号: 1815521
- 注册: 2012-05-14
- 性别: MM
- 专业: 生物力学
5楼2018-08-11 10:15:34
csp101
金虫 (著名写手)
- 应助: 9 (幼儿园)
- 金币: 9688.1
- 散金: 120
- 红花: 7
- 帖子: 1444
- 在线: 81.8小时
- 虫号: 3076622
- 注册: 2014-03-21
- 性别: GG
- 专业: 零件成形制造
6楼2018-08-11 10:16:46
lili870113
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1593.7
- 散金: 50
- 帖子: 205
- 在线: 41.5小时
- 虫号: 1815521
- 注册: 2012-05-14
- 性别: MM
- 专业: 生物力学
7楼2018-08-11 10:17:17
csp101
金虫 (著名写手)
- 应助: 9 (幼儿园)
- 金币: 9688.1
- 散金: 120
- 红花: 7
- 帖子: 1444
- 在线: 81.8小时
- 虫号: 3076622
- 注册: 2014-03-21
- 性别: GG
- 专业: 零件成形制造
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lili870113(月只蓝代发): 金币+20, 求助者申请代发金币 2018-08-13 13:38:22
月只蓝: 应助指数+1 2018-08-13 13:38:33
lili870113(月只蓝代发): 金币+20, 求助者申请代发金币 2018-08-13 13:38:22
月只蓝: 应助指数+1 2018-08-13 13:38:33
|
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
lili870113
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1593.7
- 散金: 50
- 帖子: 205
- 在线: 41.5小时
- 虫号: 1815521
- 注册: 2012-05-14
- 性别: MM
- 专业: 生物力学
9楼2018-08-11 10:39:59
