منبع : ایران سازه
دانلود راهنمای استفاده از متلب (MATLAB)
برای دانلود از ایران سازه اینجا کلیک کنید...
منبع : وبلاگ مهندسی آب
سوال) در هر ماتریس با سطر و ستون بزرگتر یا مساوی 3، عناصر داخلی وجود دارد که توسط 8 عنصر احاطه شده است. مثلا در یک ماتریس 3×3 عنصر 2,2 چنین خاصیتی دارد. برنامهای کلی بنویسید که برای یک ماتریس N×N برای هر عنصر داخلی ماتریس، عنصر ماکزیمم از عناصر محصورکننده را پیدا نموده و آن را چاپ نماید. مثال:
11 7 -2
-7 14 3
2 -3 5
عنصر داخلی این ماتریس برابر 14 و عنصر ماکزیمم محصورکننده برابر 11 است.
Integer , Dimension (:,:), Allocatable :: Mat
Integer N
Integer Maximum
Print*, "Please Enter N (Matrix is N*N):"
Read*, N
Allocate (Mat (N,N))
Print*, "Please Enter Matrix Elements (Mat (1,1), Mat(1,2), ...):"
Do I = 1, N
Do J = 1, N
Read*, Mat (I,J)
End Do
End Do
Do I = 1, N
Do J = 1, N
If (I > 1 .And. I < N .And. J > 1 .And. J < N) Then
Maximum = Max (Mat(I-1,J), Mat(I+1,J), Mat(I,J-1), Mat(I,J+1), &
Mat(I-1,J+1), Mat(I+1,J-1), Mat(I+1,J+1), Mat(I-1,J-1))
Print*, Maximum
End If
End Do
End Do
End
منبع : ایران سازه
دانلود پروژه تحلیل خرپای دو بعدی با استفاده از Matlab به همراه اصل برنامه
پروژه : تحلیل خرپای دو بعدی با استفاده از Matlab
استاد : دکتر کوهستانی
دانشجو : مصطفی فرجیان
دانشگاه تبریز
نقل قول :
برنامه نوشته شده برای تحلیل دو بعدی خرپا با استفاده از Matlab در محیط ویندوز 64bit قابل اجرا می باشد.
دانلود برنامه از ایران سازه اینجا کلیک کنید...
دانلود راهنما از ایران سازه اینجا کلیک کنید...
معرفی سایت تخصصی نرم افزار متلب (mathworks.ir)
نکته : امیدوارم در راستای مهندسی عمران فعالیت خوبی داشته باشید...
منبع : mathworks.ir
حل معادله آونگ ساده را با استفاده از روش رانگ کوتا
% simple pendelum - rung-kutta method
% saeed babanezhad
clc
clear all
%-----------------
g=9.8;
L=1;
N=200;
w=sqrt(g/L);
t=linspace(0,5,N); % time between 0-5s
h=t(2)-t(1); % time step
f=@(alpha) alpha;
g=@(teta) -w^2*sin(teta);
teta=zeros(1,N);
alpha=zeros(1,N);
teta(1)=10*pi/180; % initial angle(radian)
alpha(1)=0; % initial angle speed
for i=1:N-1
k1= f(alpha(i));
kp1= g(teta(i));
k2=f(alpha(i)+h/2*kp1);
kp2=g(teta(i)+h/2*k1);
k3=f(alpha(i)+h/2*kp2);
kp3=g(teta(i)+h/2*k2);
k4=f(alpha(i)+h*kp2);
kp4=g(teta(i)+h*k2);
teta(i+1)=teta(i)+h/6*(k1+2*k2+2*k3+k4);
alpha(i+1)=alpha(i)+h/6*(kp1+2*kp2+2*kp3+kp4);
end
plot(t,teta)
مشخصات مدیر وبسایت
عناوین یادداشتهای وبلاگ
بایگانی آرشیو ماهانه وبسایت
کلمات کلیدی وبسایت