質問編集履歴
2
コード説明書き足し
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -14,7 +14,7 @@
|
|
|
14
14
|
|
|
15
15
|
該当箇所
|
|
16
16
|
|
|
17
|
-
hp = plot3(d
|
|
17
|
+
hp = plot3(ActiveAIdx,t,conv_X(ActiveAIdx,:));
|
|
18
18
|
|
|
19
19
|
|
|
20
20
|
|
|
@@ -134,13 +134,13 @@
|
|
|
134
134
|
|
|
135
135
|
|
|
136
136
|
|
|
137
|
-
dof = [1 2];
|
|
138
137
|
|
|
139
138
|
|
|
139
|
+
ActiveADispIdx = ActiveADispIdx(1:2:end);
|
|
140
140
|
|
|
141
141
|
figure(101); hold on
|
|
142
142
|
|
|
143
|
-
hp = plot3(d
|
|
143
|
+
hp = plot3(ActiveADispIdx,t,conv_X(ActiveADispIdx,:)); %問題箇所
|
|
144
144
|
|
|
145
145
|
|
|
146
146
|
|
|
@@ -154,7 +154,9 @@
|
|
|
154
154
|
|
|
155
155
|
|
|
156
156
|
|
|
157
|
-
plot
|
|
157
|
+
同一面上で plot(t,conv_X(ActiveADisdIdx,:));
|
|
158
|
+
|
|
159
|
+
の描画はできました．
|
|
158
160
|
|
|
159
161
|
|
|
160
162
|
|
1
コード説明書き足し
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -32,7 +32,7 @@
|
|
|
32
32
|
|
|
33
33
|
|
|
34
34
|
|
|
35
|
-
cKeyName = '2dof_pch_A';
|
|
35
|
+
cKeyName = '2dof_pch_A'; %モデルデータ名
|
|
36
36
|
|
|
37
37
|
|
|
38
38
|
|
|
@@ -50,7 +50,7 @@
|
|
|
50
50
|
|
|
51
51
|
|
|
52
52
|
|
|
53
|
-
src = load([cKeyName, '.mat']);
|
|
53
|
+
src = load([cKeyName, '.mat']); %モデルデータ
|
|
54
54
|
|
|
55
55
|
|
|
56
56
|
|
|
@@ -68,7 +68,7 @@
|
|
|
68
68
|
|
|
69
69
|
|
|
70
70
|
|
|
71
|
-
dof = 2;
|
|
71
|
+
dof = 2; %自由度
|
|
72
72
|
|
|
73
73
|
|
|
74
74
|
|
|
@@ -100,7 +100,7 @@
|
|
|
100
100
|
|
|
101
101
|
|
|
102
102
|
|
|
103
|
-
% solve
|
|
103
|
+
% solve %インパルス応答計算
|
|
104
104
|
|
|
105
105
|
A0 = A(ActiveAIdx, ActiveAIdx);
|
|
106
106
|
|
|
@@ -112,13 +112,11 @@
|
|
|
112
112
|
|
|
113
113
|
|
|
114
114
|
|
|
115
|
-
%% convolution
|
|
115
|
+
%% convolution %畳み込み積分
|
|
116
116
|
|
|
117
|
-
|
|
117
|
+
F = 5*cos(7*t); %調和外力
|
|
118
118
|
|
|
119
|
-
F = 5*cos(7*t);
|
|
120
|
-
|
|
121
|
-
conv_X = zeros(2*dof, iNumStep);
|
|
119
|
+
conv_X = zeros(2*dof, iNumStep);
|
|
122
120
|
|
|
123
121
|
|
|
124
122
|
|