[]
MatLab
merhaba,
matlab'da output'u
u1=(8131307261700661*x)/72057594037927936
şeklinde aliyorum.
0.1128*x şeklinde nasil alabilirim ?
sadeleştirmeyi nasil yaparim?
matlab'da output'u
u1=(8131307261700661*x)/72057594037927936
şeklinde aliyorum.
0.1128*x şeklinde nasil alabilirim ?
sadeleştirmeyi nasil yaparim?
sembolik matematik mi çalışıyorsun ? bu çıktıyı veren kod ne ?
- kimlanbu
(30.05.11 15:17:29)
evet sembolik çalişiyorum. yeni yeni isindiriyorum kendimi matlaba. script de şu şekilde:
clear all
clc
syms x
k=(2+3*x)*cos(pi*x/3);
f=sin(pi*x/3);
%shape functions
fi = [x;x^3;x^5];
%stiffness matrices
for i=1:3
for j=1:3
K(i,j)=0;
K(i,j)=int(k*diff(fi(i,1),x)*diff(fi(j,1),x),x,0,1);
end
end
%force vectors
for i=1:3
F(i,1)=int(f*fi(i,1),x,0,1);
end
%%Ritz Coefficients and Displacements
%for N=1
a1=K(1,1)^-1*F(1,1);
a1 = double(a1)
u1=a1*fi(1,1)
%for N=2
a2=K(1:2,1:2)^-1*F(1:2,1);
a2=double(a2)
u2=a2(1,1)*fi(1,1)+a2(2,1)*fi(2,1)
%for N=3
a3=K(1:3,1:3)^-1*F(1:3,1);
a3=double(a3)
u3=a3(1,1)*fi(1,1)+a3(2,1)*fi(2,1)+a3(3,1)*fi(3,1)
clear all
clc
syms x
k=(2+3*x)*cos(pi*x/3);
f=sin(pi*x/3);
%shape functions
fi = [x;x^3;x^5];
%stiffness matrices
for i=1:3
for j=1:3
K(i,j)=0;
K(i,j)=int(k*diff(fi(i,1),x)*diff(fi(j,1),x),x,0,1);
end
end
%force vectors
for i=1:3
F(i,1)=int(f*fi(i,1),x,0,1);
end
%%Ritz Coefficients and Displacements
%for N=1
a1=K(1,1)^-1*F(1,1);
a1 = double(a1)
u1=a1*fi(1,1)
%for N=2
a2=K(1:2,1:2)^-1*F(1:2,1);
a2=double(a2)
u2=a2(1,1)*fi(1,1)+a2(2,1)*fi(2,1)
%for N=3
a3=K(1:3,1:3)^-1*F(1:3,1);
a3=double(a3)
u3=a3(1,1)*fi(1,1)+a3(2,1)*fi(2,1)+a3(3,1)*fi(3,1)
- okkusapik
(30.05.11 15:21:39 ~ 15:36:54)
evde değilim kodu çalıştıramıyorum, hata vermesi muhtemel. u1 hesaplandıktan sonra aşağıdaki satırı dene bir...
u1= sym(u1, 'd')
u1= sym(u1, 'd')
- kimlanbu (30.05.11 16:50:01 ~ 16:52:23)
hata verdi maalesef
??? Error using ==> sym.sym>assumptions at 2196
Second argument d not recognized.
Error in ==> sym.sym>tomupad at 2172
assumptions(S,x.s,a);
Error in ==> sym.sym>sym.sym at 123
S.s = tomupad(x,a);
Error in ==> hw3new at 32
u1== sym(u1, 'd')
??? Error using ==> sym.sym>assumptions at 2196
Second argument d not recognized.
Error in ==> sym.sym>tomupad at 2172
assumptions(S,x.s,a);
Error in ==> sym.sym>sym.sym at 123
S.s = tomupad(x,a);
Error in ==> hw3new at 32
u1== sym(u1, 'd')
- okkusapik (30.05.11 16:52:17 ~ 16:53:16)
ya o pi orada durmak zorunda mı? pi bi sabit olarak şey edip en sonda hesaba sokmayı bi dener misin?
- bl (31.05.11 15:36:19)
öyle submit ettim artik ödevi eyvallah denemeye gerek kalmadi
- okkusapik (04.06.11 20:09:45)
1
