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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 请教matlab用数组给参数矩阵赋值的问题
51/1返回列表
查看: 2405  |  回复: 11
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

加油大楠楠

铜虫 (小有名气)

[求助] 请教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数列中只有一个数 ，请帮忙
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2013-04-15 09:14:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

加油大楠楠

铜虫 (小有名气)
送红花一朵
引用回帖:
7楼: Originally posted by 材料廖 at 2013-04-21 18:10:57
开始那个由于你在rad中用了0，导致了出现了NAN，我照你的程序修改了下，可以出部分结果，程序如下：
j=1;z=20;r=32;
for l=32:34   %设定范围
    i=1;  tezheng=0;
    for rad=1:l; %这两个for的步长可以根据 ...

我大概明白了  再厚脸皮问你一下 这个hh矩阵能不能用极坐标画出来 画成三维的图，就是图中x坐标为极轴长，y坐标为角度，hh中对应的数值为z值
赞一下 回复此楼
8楼2013-04-22 20:32:35
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

材料廖

木虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
你这个肯定只有n=1的值呀，都没有循环，你可以将你的这段程序
f(n,1)=e;%将范数记录
n=n+1;
if n==101;%B数组有100行所以到一百终止
end
改成：
for n = 1:100
f(n,1) = e;
end
这样应该就可以了！
赞一下 回复此楼
2楼2013-04-15 12:12:41
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

加油大楠楠

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by 材料廖 at 2013-04-15 12:12:41
你这个肯定只有n=1的值呀，都没有循环，你可以将你的这段程序
f(n,1)=e;%将范数记录
n=n+1;
if n==101;%B数组有100行所以到一百终止
end
改成：
for n = 1:100
f(n,1) = e;
end
这样应该就可以了！

不好意思回复晚了  那你再帮我看看这整段程序   我把上次的程序稍微改了下，就是对于每个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
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

材料廖

木虫 (正式写手)

【答案】应助回帖

我还是不是很明白你要做什么，你这个程序写的我有点看不懂。但有点建议，一：你这个应该没有必要用符号计算，符号计算很慢。二：你这个程序编写的不怎样，基础不是很好，写的有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
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭