| 查看: 1219 | 回复: 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总提示索引超出矩阵维度 |
回复此楼» 猜你喜欢
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
梦想
已经有7人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
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
5