How do I manually insert variables I have declared symbolically into a zero matrix?

2018 8 23
1 个回答
回答已采纳
更新时间：2018 8 24
2 次查看（30 天）

0 个投票

I have been trying to insert symbolic variable into a zero matrix.
clear all
clc
syms k1 k2 k3 k4
global_k1 = zeros(4,4);
global_k2 = zeros(4,4);
global_k3 = zeros(4,4);
global_k4 = zeros(4,4);
global_k1(1,1) = k1;
global_k1(1,2) = -k1;
global_k1(2,1) = -k1;
global_k1(2,2) = k1;
global_k2(2,2) = k2;
global_k2(2,3) = -k2;
global_k2(3,2) = -k2;
global_k2(3,3) = k2;
global_k3(2,2) = k3;
global_k3(2,3) = -k3;
global_k3(3,2) = -k3;
global_k3(3,3) = k3;
global_k4(2,2) = k4;
global_k4(2,4) = -k4;
global_k4(4,2) = -k4;
global_k4(4,4) = k4;
k_stiff = global_k1+global_k2+global_k3+global_k4;
disp('The stiffness matrix is: ');
disp(k_stiff);
It displays the following error when I run the program:
The following error occurred converting from sym to double:
Unable to convert expression into double array.
Error in HW1_1 (line 9)
global_k1(1,1) = k1;
How should I resolve this error?

请先登录，再回答此问题。

 采纳的回答

madhan ravi
madhan ravi 2018-8-24
编辑：madhan ravi 2018-8-24

0 个投票

>> syms k1 k2 k3 k4
global_k1 = sym(zeros(4,4));
global_k2 = sym(zeros(4,4));
global_k3 = sym(zeros(4,4));
global_k4 = sym(zeros(4,4));
global_k1(1,1) = k1;
global_k1(1,2) = -k1;
global_k1(2,1) = -k1;
global_k1(2,2) = k1;
global_k2(2,2) = k2;
global_k2(2,3) = -k2;
global_k2(3,2) = -k2;
global_k2(3,3) = k2;
global_k3(2,2) = k3;
global_k3(2,3) = -k3;
global_k3(3,2) = -k3;
global_k3(3,3) = k3;
global_k4(2,2) = k4;
global_k4(2,4) = -k4;
global_k4(4,2) = -k4;
global_k4(4,4) = k4;
global_k1
global_k2
global_k3
global_k4
k_stiff = global_k1+global_k2+global_k3+global_k4;
disp('The stiffness matrix is: ')
disp(k_stiff);
global_k1 =
[ k1, -k1, 0, 0]
[ -k1, k1, 0, 0]
[ 0, 0, 0, 0]
[ 0, 0, 0, 0]
global_k2 =
[ 0, 0, 0, 0]
[ 0, k2, -k2, 0]
[ 0, -k2, k2, 0]
[ 0, 0, 0, 0]
global_k3 =
[ 0, 0, 0, 0]
[ 0, k3, -k3, 0]
[ 0, -k3, k3, 0]
[ 0, 0, 0, 0]
global_k4 =
[ 0, 0, 0, 0]
[ 0, k4, 0, -k4]
[ 0, 0, 0, 0]
[ 0, -k4, 0, k4]
The stiffness matrix is:
[ k1, -k1, 0, 0]
[ -k1, k1 + k2 + k3 + k4, - k2 - k3, -k4]
[ 0, - k2 - k3, k2 + k3, 0]
[ 0, -k4, 0, k4]

13 个评论
显示 11更早的评论 隐藏 11更早的评论

Salil Chawla
Salil Chawla 2018-8-24
I wanted to insert the indexed values in the zero matrix, not just add them. Thanks anyway
Salil Chawla
Salil Chawla 2018-8-24
Not the result I intended. The addition is based on the indexed position. The way you suggested would just add the matrices
madhan ravi
madhan ravi 2018-8-24
Can you show the result you intend? Just to understand what you actually want.
madhan ravi
madhan ravi 2018-8-24
Plus your actual question was How to resolve the error?
Salil Chawla
Salil Chawla 2018-8-24
This is basically to find the elemental and global stiffness matrix for finite element analysis.
madhan ravi
madhan ravi 2018-8-24
编辑：madhan ravi 2018-8-24
Can you show an example of the indexed addition ? If you define the process what has to be found ? It would be narrowed instead of generalisation because it’s totally vague for now.
Salil Chawla
Salil Chawla 2018-8-24
Something of the effect of this:
clear all
clc
syms k1 k2 k3 k4
global_k1 = [k1 -k1 0 0; -k1 k1 0 0; 0 0 0 0; 0 0 0 0];
global_k2 = [0 0 0 0; 0 k2 -k2 0; 0 -k2 k2 0; 0 0 0 0];
global_k3 = [0 0 0 0; 0 k3 -k3 0; 0 -k3 k3 0; 0 0 0 0];
global_k4 = [0 0 0 0; 0 k4 0 -k4 ; 0 0 0 0; 0 -k4 0 k4];
k_stiff = global_k1 + global_k2 + global_k3 + global_k4;
disp('The stiffness matrix is: ');
disp(k_stiff);
The difference is that here I had to manually input the symbolic values into a zero matrix and then add them. What I intended to do was to declare zero matrices for variable global_k1, k2 etc and then put variables into the matrices by calling variables into different positions by calling the function
global_k1(1,1) = k1;
global_k1(1,2) = -k1;
global_k1(2,1) = -k1;
global_k1(2,2) = k1;
And so on. Hope that clears the picture.
madhan ravi
madhan ravi 2018-8-24
编辑：madhan ravi 2018-8-24
Code edited - first zero matrices were made and then symbolic variables inserted according to the index.If still you didn’t get the required result, just type the desired result in the next comment. I’ll give you the apt code.
Salil Chawla
Salil Chawla 2018-8-24
The stiffness matrix is:
[ k1, -k1, 0, 0]
[ -k1, k1 + k2 + k3 + k4, - k2 - k3, -k4]
[ 0, - k2 - k3, k2 + k3, 0]
[ 0, -k4, 0, k4]
This is the result I am looking for
madhan ravi
madhan ravi 2018-8-24
see the code I just posted the results you wanted and it works fine for me for your information , I am using Matlab 2017a release .

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Assumptions 的更多信息

产品

版本

R2017b

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by