| 查看: 2877 | 回复: 11 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
请教matlab用数组给参数矩阵赋值的问题
|
||
|
syms x y z d1 d2 d3 r l;n=1;%定义矩阵参数 a=1/l*[z+r-d1 0 0;0 -0.866*y-0.5*z+r-d2 0;0 0 0.866*y-0.5*z+r-d3]; b=[x y z+r-d1;x y-0.866*r+0.866*d2 z-0.5*r+0.5*d2;x y+0.866*r-0.866*d3 z-0.5*r+0.5*d3]; c=inv(a); d=c*b;%d为我需要的比矩阵 到这一步矩阵的计算都不成问题 d=subs(d,{'r','l','x','y','z','d1','d2','d3'},{31.8,40.2,B(n,1),B(n,2),B(n,3),B(n,4),B(n,5),B(n,6)}); %给参数赋值，其中x y z d1 d2 d3是用B矩阵的第n行对应的列来赋值的 e = norm(d,'fro'); %求范数 f(n,1)=e;%将范数记录 n=n+1; if n==101;%B数组有100行所以到一百终止 end 最后结果f数列中只有一个数 ，请帮忙 |
回复此楼» 猜你喜欢
只有每年这种时候来逛逛小木虫
已经有25人回复
-
什么时候开奖？
已经有7人回复
-
看来今天不会放榜了？
已经有3人回复
-
时间戳又变了
已经有18人回复
-
我面上完蛋了
已经有7人回复
-
时间戳又变了8-15
已经有27人回复
-
今天基金会出结果吗？20260819
已经有16人回复
-
科研孤儿太难了
已经有16人回复
-
应该是下周三26日公布了吧？
已经有4人回复
-
投稿咨询
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab稀疏矩阵问题求教
已经有3人回复
- 如何用matlab求解矩阵系数的二阶微分方程
已经有13人回复
- matlab函数调用与数组的问题
已经有19人回复
- 用matlab来实现分类汇总
已经有3人回复
- 如何用Matlab画矩阵的散点图？
已经有3人回复
- 菜鸟求助~matlab怎么样让一个1xn的向量中的每个元素是个矩阵？
已经有11人回复
- 请教向量Q满足高斯分布，用matlab如何表示？
已经有11人回复
- 请教关于matlab矩阵转换问题
已经有6人回复
- matlab线性矩阵不等式求解做优化问题
已经有10人回复
- 如何用c++来求解矩阵的相关问题。
已经有12人回复
- 【分享】Matlab经典教程——从入门到精通【已搜索无重复】
已经有797人回复
- VS2005下matlab调用
已经有8人回复
- 【求助】请教在在matlab里面怎么输入含参数矩阵指数形式-----------【已完结】
已经有5人回复
- 【求助】关于matlab的高维数组
已经有6人回复
- 【求助】初学matlab
已经有11人回复
- 【求助】柱形统计图matlab怎么画？
已经有4人回复
- 【求助】MATLAB的矩阵扩展问题
已经有2人回复
材料廖
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 4889.5
- 散金: 200
- 红花: 11
- 帖子: 891
- 在线: 295.1小时
- 虫号: 1794189
- 注册: 2012-05-04
- 性别: GG
- 专业: 金属材料的合金相、相变及
9楼2013-04-22 22:01:14
材料廖
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 4889.5
- 散金: 200
- 红花: 11
- 帖子: 891
- 在线: 295.1小时
- 虫号: 1794189
- 注册: 2012-05-04
- 性别: GG
- 专业: 金属材料的合金相、相变及
2楼2013-04-15 12:12:41
|
不好意思回复晚了 那你再帮我看看这整段程序 我把上次的程序稍微改了下，就是对于每个l，我要求他xyx为不同值的时候的范数，然后相加，最后把每个l的范数和提取出来，十分感谢了 我运行总是报错 j=1;z=20;r=32; for l=32:34 %设定范围 i=1; tezheng=0; for rad=0:l; %这两个for的步长可以根据运算效率适当放大 for thea=0:pi/10:2*pi; m=rad*cos(thea); n=rad*sin(thea); d11=1/3*(3*r+3*m+sqrt(9*l^2-9*n^2-9*z^2)); d12=1/3*(6*r-3*m+3*n*sqrt(3)+sqrt(36*l^2-27*m^2-18*m*n*sqrt(3)-9*n^2-36*z^2)); d13=1/3*(6*r-3*m-3*n*sqrt(3)+sqrt(36*l^2-27*m^2+18*m*n*sqrt(3)-9*n^2-36*z^2)); B(i,  =[m,n,z,d11,d12,d13];syms x1 y1 z1 d1 d2 d3 l1; a=(1/l1)*[z1+r-d1 0 0;0 -0.866*y1-0.5*z1+r-d2 0;0 0 0.866*y1-0.5*z1+r-d3]; b=[x1 y1 z1+r-d1;x1 y1-0.866*r+0.866*d2 z1-0.5*r+0.5*d2;x1 y1+0.866*r-0.866*d3 z1-0.5*r+0.5*d3]; c=inv(a); d=c*b;%d为雅克比矩阵 d=subs(d,{'l1','x1','y1','z1','d1','d2','d3'},{l,B(i,1),B(i,2),B(i,3),B(i,4),B(i,5),B(i,6)}); e = norm(d,'fro'); %2范数 hh(i,1)=1/e;%把这个范数记录在hh中 tezheng=hh(i,1)+tezheng;%相加对于相同l值的对应的范数 i=i+1; end end tezheng1(j, =tezheng/(i-1);%把不同l值对应的范数求和再除以相应的l值 存在数组tezheng1中j=j+1; end 错误信息是 clear ??? Error using ==> mupadmex Error in MuPAD command: Out of memory Error in ==> sym.horzcat at 25 y = mupadmex('mllib::horzcat',strs{:}); Error in ==> jacobi at 16 a=(1/l1)*[z1+r-d1 0 0;0 -0.866*y1-0.5*z1+r-d2 0;0 0 0.866*y1-0.5*z1+r-d3]; >> |
3楼2013-04-18 09:38:21
材料廖
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 4889.5
- 散金: 200
- 红花: 11
- 帖子: 891
- 在线: 295.1小时
- 虫号: 1794189
- 注册: 2012-05-04
- 性别: GG
- 专业: 金属材料的合金相、相变及
【答案】应助回帖
|
我还是不是很明白你要做什么，你这个程序写的我有点看不懂。但有点建议，一：你这个应该没有必要用符号计算，符号计算很慢。二：你这个程序编写的不怎样，基础不是很好，写的有C语言的感觉，可能你的想法你没用程序表达出来，要不你说说你要干啥吧，我看能不能帮你。话说直了点，望不怪。 下面是我没用符号计算程序，速度快了很多，但是还是存在错误，出不了结果，因为我还是没有理解你要干啥。 j=1;z=20;r=32; for l=32:34 %设定范围 i=1; tezheng=0; for rad=0:l; %这两个for的步长可以根据运算效率适当放大 for thea=0:pi/10:2*pi; m=rad*cos(thea); n=rad*sin(thea); d11=1/3*(3*r+3*m+sqrt(9*l^2-9*n^2-9*z^2)); d12=1/3*(6*r-3*m+3*n*sqrt(3)+sqrt(36*l^2-27*m^2-18*m*n*sqrt(3)-9*n^2-36*z^2)); d13=1/3*(6*r-3*m-3*n*sqrt(3)+sqrt(36*l^2-27*m^2+18*m*n*sqrt(3)-9*n^2-36*z^2)); B(i, =[m,n,z,d11,d12,d13];a=(1/B(i,1))*[B(i,3)+r-B(i,4) 0 0;0 -0.866*B(i,2)-0.5*B(i,3)+r-B(i,5) 0;0 0 0.866*B(i,2)-0.5*B(i,3)+r-B(i,6)]; b=[B(i,1) B(i,2) B(i,3)+r-B(i,4);B(i,1) B(i,2)-0.866*r+0.866*B(i,5) B(i,3)-0.5*r+0.5*B(i,5);B(i,1) B(i,2)+0.866*r-0.866*B(i,6) B(i,3)-0.5*r+0.5*B(i,6)]; % c=inv(a); d=a\b;%d为雅克比矩阵 e = norm(d,'fro'); %2范数 hh(i,1)=1/e;%把这个范数记录在hh中 tezheng=hh(i,1)+tezheng;%相加对于相同l值的对应的范数 i=i+1; end end tezheng1(j, =tezheng/(i-1+eps);%把不同l值对应的范数求和再除以相应的l值 存在数组tezheng1中j=j+1; end |
4楼2013-04-19 18:36:16
