مهندسی مکانیک و هوافضا
هر آنچه که یک مهندس بخواهد اینجا هست!
 
آخرین مطالب
 
محبوبترین ها ی سایت
آدرس جدید سایت WWW.AEROS.IR


1-اگر خواهان تبادل لینک با ما هستید پس از گذاشتن لینک ما در وبلاگ خود ، آدرس وبلاگ خود را برای ما در زیر همین پست کامنت کنید . 

2-استفاده بدون مبنع از مطالب وبلاگ بلامانع است .




برچسب‌ها: WWW.AEROS.IR ,
نرم کردن نمودار در متلب ( smoothing plot )

smooth

Smooth response data

Syntax

yy = smooth(y)
gpuarrayYY = smooth(gpuarrayY)
yy = smooth(y,span)
yy = smooth(y,method)
yy = smooth(y,span,method)
yy = smooth(y,'sgolay',degree)
yy = smooth(y,span,'sgolay',degree)
yy = smooth(x,y,...)

Description

yy = smooth(y) smooths the data in the column vector y using a moving average filter. Results are returned in the column vector yy. The default span for the moving average is 5.

The first few elements of yy are given by

yy(1) = y(1)
yy(2) = (y(1) + y(2) + y(3))/3
yy(3) = (y(1) + y(2) + y(3) + y(4) + y(5))/5
yy(4) = (y(2) + y(3) + y(4) + y(5) + y(6))/5
...

Because of the way endpoints are handled, the result differs from the result returned by the filter function.

gpuarrayYY = smooth(gpuarrayY) performs the operation on a GPU. The input gpuarrayY is a gpuArray column vector. The output gpuarrayYY is a gpuArray column vector. This syntax requires the Parallel Computing Toolbox™.

    Note:   You can use gpuArray x and y inputs with the smooth function, but this is only recommended with the default 'method', 'moving'. Using GPU data with other methods does not offer any performance advantage.

yy = smooth(y,span) sets the span of the moving average to spanspan must be odd.

yy = smooth(y,method) smooths the data in y using the method method and the default span. Supported values for method are listed in the table below.

method

Description

'moving'

Moving average (default). A lowpass filter with filter coefficients equal to the reciprocal of the span.

'lowess'

Local regression using weighted linear least squares and a 1st degree polynomial model

'loess'

Local regression using weighted linear least squares and a 2nd degree polynomial model

'sgolay'

Savitzky-Golay filter. A generalized moving average with filter coefficients determined by an unweighted linear least-squares regression and a polynomial model of specified degree (default is 2). The method can accept nonuniform predictor data.

'rlowess'

A robust version of 'lowess' that assigns lower weight to outliers in the regression. The method assigns zero weight to data outside six mean absolute deviations.

'rloess'

A robust version of 'loess' that assigns lower weight to outliers in the regression. The method assigns zero weight to data outside six mean absolute deviations.

yy = smooth(y,span,method) sets the span of method to span. For the loess and lowess methods, span is a percentage of the total number of data points, less than or equal to 1. For the moving average and Savitzky-Golay methods, span must be odd (an even span is automatically reduced by 1).

yy = smooth(y,'sgolay',degree) uses the Savitzky-Golay method with polynomial degree specified by degree.

yy = smooth(y,span,'sgolay',degree) uses the number of data points specified by span in the Savitzky-Golay calculation. span must be odd and degree must be less than span.

yy = smooth(x,y,...) additionally specifies x data. If x is not provided, methods that require x data assume x = 1:length(y). You should specify x data when it is not uniformly spaced or sorted. If x is not uniform and you do not specify methodlowess is used. If the smoothing method requires x to be sorted, the sorting occurs automatically.

Examples

                                                                           read more..

ادامه مطلب

مرتبط با: متلب ( MATLAB ) ,
مراحل ایجاد کلید نیم‌فاصله در نرم‌افزار WORD

تدا نرم‌افزار Word2007 یا Word2010 را باز می‌کنیم و سپس مراحل زیر را یک به یک اجرا می‌نماییم.

1.    بر روی نوار Insert کلیک می‌کنیم و سپس در انتهای نوار سراغ گزینه Symbol  می‌رویم و آن را انتخاب می‌کنیم.


 

به ادامه ی مطلب مراجعه کنید

 


ادامه مطلب

مرتبط با: آموزش آفیس ,
آیرودینامیک فعال در طرح مفهومی مرسدس بنز IAA

انسپت افراطی جدید مرسدس‌بنز، نشانه‌هایی از آینده طراحی‌های این کمپانی در خود دارد. این کانسپت سدان مجهز به برخی اِلمان‌های جدید بدنه به صورت هوشمند است که می‌توانند دو حالت متفاوت آئرودینامیک برای خودرو ایجاد کنند.

 

این کانسپت مرسدس‌بنز که با ضریب دراگ 0.19 به مانند ماهی در هوا جریان می‌یابد، اینتلیجنت آئرودینامیک اتومبیل (IAA) نام دارد که قابلیت تبدیل از حالت «دیزاین» به حالت نرم‌تر «آئرودینامیک» را در بدنه خود دارد. در حالت آئرودینامیک، که با فشار یک دکمه توسط راننده و یا به صورت اتوماتیک در سرعت ۸۰ کیلومتر بر ساعت فعال می‌شود، ۸ عضو بدنه در عقب به حالت کام (Kamm-tail) درمی‌آیند و طول کلی خودرو به اندازه ۱۵/۴ اینچ یا ۳۹۱ میلی‌متر افزایش می‌دهند. همچنین باله‌ها روی سپر جلو به اندازه یک اینچ اضافه می‌شوند تا بتوانند جریان هوای بیشتری را به سمت عقب منتقل کنند. همین‌طور چرخ‌ها نیز «مقطع» خود را به اندازه حدود ۲/۲ اینچ (۵۶ میلی‌متر) کاهش می‌دهند و هواکش‌های جلو نیز ۲/۴ اینچ (۶۱ میلی‌متر) عقب می‌نشینند. این تبدیل‌های آئرودینامیک در نهایت منجر به یک حالت بسیار بهینه در بدنه می‌شود که کمترین میزان دراگ ممکن را ایجاد خواهد کرد.

به ادامه مطلب مراجعه کنید

ادامه مطلب

مرتبط با: تونل باد , هوافضا ,
ترسیم نمودار در متلب - 4

LineSpec (Line Specification)

Line specification string syntax

Description

Plotting functions accept string specifiers as arguments and modify the graph generated accordingly. Three components can be specified in the string specifiers along with the plotting command. They are:

  • Line style

  • Marker symbol

  • Color

For example:

plot(x,y,'-.or')

plots y versus x using a dash-dot line (-.), places circular markers (o) at the data points, and colors both line and marker red (r). Specify the components (in any order) as a quoted string after the data arguments. Note that linespecs are single strings, not property-value pairs.

Plotting Data Points with No Line

If you specify a marker, but not a line style, only the markers are plotted. For example:

plot(x,y,'d')

read more..

به ادامه مطلب مراجعه کنید

ادامه مطلب

مرتبط با: متلب ( MATLAB ) ,
همه ی دستورات متلب

عملگرهای ریاضی و منطقی

+ , - , / , * , ^ , > , < , >= , <= , == , ~= , ~ , || , &&

 

مدیریت m فایل

% , … , %%

 

رسم نمودارهای دو بعدی

plot , fplot , comet , bar , barh, stairs , stem , pie , semilogx , semilogy , loglog , hist , polar

 

رسم نمودارهای سه بعدی

plot3 , meshgrid , mesh , surf , contour , comet3 , contour3 , meshc , clabel , surfl , …

 

متعلقات مربوط به نمودارهای دو بعدی و سه بعدی

xlabel , ylabel , zlabel , title , legend , grid , text , gtext , hold , figure , subplot , zoom , axis , xlim , ylim , clf , cla , close , colorbar , colormap , axes , ginput ,  …

 

انتگرالگیری عددی

trapz , quad , dblquad , triplequad , …

 

حل معادله دیفرانسیل عددی

ode (ode45 , ode23 , …) bvp4c , bvp5c , odeset , …

به ادامه ی مطلب مراجعه کنید

 

ادامه مطلب

مرتبط با: متلب ( MATLAB ) ,
مقایسه دو رشته (string) با یکدیگر، برای چک کردن مثل هم بودن آنها، به صورت حساس به بزرگی و کوچکی حروف (case sensitive)، با دستور strcmp در متلب

با استفاده از دستور strcmp در متلب، می توانیم دو رشته (string) را با یکدیگر مقایسه کنیم و ببینیم که آیا دقیقا مثل هم هستند یا خیر. شیوه مقایسه کردن، به صورت حساس بودن به بزرگی و کوچکی حروف (case sensitive) نیز خواهد بود، یعنی برای مثل هم بودن، تنها معیار مثل هم بودن حروف آنها نیست، بلکه از نظر بزرگ بودن یا کوچک بودن حروف نیز باید مثل هم باشند.

اگر دو رشته (string) دقیقا مثل هم باشند، دستور strcmp ، مقدار 1 را برمی گرداند، اما اگر دو رشته متفاوت باشند، دستور strcmp ، مقدار 0 را برمی گرداند.

به مثال زیر توجه کنید :

مثال

clear all
close all
clc

A = strcmp('Yes''No')
B = strcmp('Yes''Yes')
C = strcmp('YES''yes')

سه خط اول کدها، برای عدم تداخل برنامه فعلی با برنامه قبلی اجرا شده در متلب می باشد.

نتیجه :


A =

     0


B =

     1


C =

     0

منبع : کلیدستان


مرتبط با: متلب ( MATLAB ) ,
ترسیم نمودار در متلب - 3

قبل از اینکه بخواهم نحوه ایجاد نمودارهای 2 بعدی رو به شما عزیزان بگم باید یادآوری کنم فرض من بر این هست که شما دوستان مباحث ابتدایی متلب مثل ایجاد بردارها و ماتریس ها و کار کردن با m فایل ها را بلدید. اما ان شاء الله در آینده من ابتدایی ترین مباحث رو هم خدمتتون ارائه خواهم کرد. اما فعلاً نه!

فرض کنید که می خواهیم نمودار سینوس  را در فاصله بین صفر تا 2*pi رسم کنیم:

t=0:0.01:2*pi;

u=sin(t);

گام متغیر t  را برابر با 0.01 قرار دادم تا وقتیکه نمودار سینوس کشیده میشه حالت شکسته نداشته باشه.

دستور رسم نمودارهای دو بعدی در متلب دستور plot هستش:

plot(t,u)

با اجرای دستور بالا میبینید که نمودار سینوس کشیده میشه.

دستوراتی که از الان به بعد می گویم حتما باید زیر دستور plot نوشته بشه.

grid on    یا    grid

دستور grid on یا grid باعث شبکه بندی صفحه مختصات میشه.

xlabel('t')

ylabel('u')

title('graph u vs t ')

legend('sin')

دستور xlabel باعث گذاشتن برچسب بر روی محور x ها میشه. بین دو علامت ' '  هرچی بنویسید در روی محور x ها نشون داده میشه.

دستور ylabel باعث گذاشتن برچسب بر روی محور y ها، دستور title باعث گذاشتن عنوان بر روی نمودار و دستور legend باعث گذاشتن راهنما در گوشه نمودار میشه. بین دو علامت ' ' هرچی بنویسید همون نشون داده خواهد شد.

فرض کنید که می خواهید در روی بر چسب محور x علامت ریاضی آلفا یا بتا یا هر حرف یونانی دیگه قرار داده بشه. برای این کار در بین دو علامت ' ' کلمه alpha یا beta را تایپ میکنید و قبلش علامت \ را می گذارید:

xlabel('\alpha')

title('\beta')

حال فرض کنید که می خواهید اندیس بالا یا پایین در روی نوشته هاتون داشته باشید. اندیس پایین با علامت _ ایجاد میشه و اندیس بالا با علامت ^ ایجاد میشه:

xlabel('x_i')

legend('R^2')

اما فرض کنید که می خواهید بر روی برچسب محور x (یا محور y یا در عنوان و یا در راهنما) عبارت R را بنویسید که در اندیس پایین اون عبارت index باشه. شاید فکر کنید باید به صورت زیر بنویسید:

xlabel(R_index)

اما با این کار فقط حرف i به صورت اندیس پایین نشون داده خواهد شد. و یا مثلا اگر که بخواهید index اندیس بالا نوشته بشه دستور زیر فقط حرف i را در اندیس بالا نشان میده.

legend('R^index')

برای اینکه این دستور درست کار بکنه 2 راه وجود داره که راه اول کمی غیر منطقی به نظر میرسه:

راه اول:

xlabel(R_i_n_d_e_x)

legend(R^i^n^d^e^x)

این کار برای اندیسهایی که تعداد حروفشون کمه شاید به نظر غیر منطقی نیاد اما اگرکه کلمه اندیس بالا یا پایین 50 حرف باشه غیر منطقیه که بین 50 حرف از _ یا ^ استفاده کنیم. برای همین بهتره از راه دوم استفاده کنیم.

راه دوم: در این راه عبارت اندیس را بین دو آکولاد قرار می دهیم:

xlabel(R_{index})

ylabel(R^{index})

فرض کنید می خواهید در راهنمای نمودار خود عبارت آلفا به توان سیگما را نشان دهید:

legend(\alpha^{\sigma})

سایر علائم و حروف یونانی رو می تونید در قسمت help با جستجوی عبارت text properties در یک جدول کاملی پیدا کنید.

 

اما متوجه شدید که نمودار وقتی کشیده شد با رنگ پیش فرض یعنی آبی کشیده شد. اگر که بخواهید رنگ نمودار را تغییر دهید باید در همان دستور plot این کار را انجام دهید:

راه اول: تغییر رنگ به یکی از رنگهای سیاه، سفید، فیروزه ای، قرمز، سبز، آبی، زرد و بنفش.

plot(t,u,' ')

بین دو علامت ' ' هر کدام از حروف زیر را که بنویسید رنگ متناسب با حرف نوشته شده تغییر خواهد کرد:

c : فیروزه ای

m : بنفش

b: آبی (که اگر نوشته نشود هم رنگ آبی پیش فرض می باشد)

r: قرمز

g: سبز

k: سیاه

w: سفید

y: زرد

برای مثال می خواهید رنگ نمودار قرمز شود:

plot(t,u,'r')

راه دوم: می خواهید از رنگهای ترکیبی دیگر استفاده کنید:

plot(t,u,'color',[r g b])

به جای سه عدد r و g و b  باید عددی بین 0 و 1 وارد کنید. عدد مربوط به r معرف رنگ قرمز، عدد مربوط به g معرف رنگ سبز و عدد مربوط به b معرف رنگ آبی میباشد. مسلما اگر سه عدد را به صورت [0 0 1] قرار دهید رنگ به صورت قرمز خواهد شد و اگر هر سه را برابر با 1 قرار دهید مشکی و اگر هر سه را برابر با 0 قرار دهید سفید خواهد شد. پش با انتخاب سه عدد بین 0 و 1 میلیونها رنگ می توانید داشته باشید:

plot(t,u,'color',[1 0.2 0.06])

بار دیگر به دستور ابتدایی بر می گردییم. مشاهده کردید که وقتی نمودار را کشیدید با یک ضخامت پیش فرض نمودار برای شما کشیده شد. اگر که بخواهید ضخامت را هم عوض کنید باید از دستور LineWidth در دستور plot به صورت زیر استفاده کنید:

plot(t,u,'LineWidth',a)

یا

plot(t,u,'g','LineWidth',a)

یا

plot(t,u,'color',[1 .03 .5],'LineWidth',a)

منظور از a عددی است که ضخامت خط نمودار را  تعیین می کنه مثلا a=2:




مرتبط با: متلب ( MATLAB ) ,
ترسیم نمودار در متلب - 2

کدهای متلب :


clear all
close all
clc

x = linspace(-2,2,20);
y = 1./(x+(1-x).^2);
p = polyfit(x,y,3);
xx = linspace(-2,2,100);
plot(x,y,'.',xx,polyval(p,xx),'-','MarkerSize',30,'LineWidth',2)

نتیجه :

مثال شماره 16

منبع کدها :

کتاب MATLAB Guide - D.J. Higham - second edition - صفحه 161 از فایل pdf کتاب



مرتبط با: متلب ( MATLAB ) ,
ترسیم نمودار در متلب - 1

کدهای متلب :


clear all
close all
clc

x = [0 pi/4 3*pi/8 3*pi/4 pi]; y = sin(x);
xi = linspace(0,pi,40)';
yn = interp1(x,y,xi,'nearest');
yl = interp1(x,y,xi,'linear');
ys = interp1(x,y,xi,'spline');
yp = interp1(x,y,xi,'pchip');
xx = linspace(0,pi,50);
plot(xi,yn,'*', xi,yl,'+', xi,ys,'v', xi,yp,'o')
legend('nearest','linear','spline','pchip')
hold on
plot(xx,sin(xx),'-',x,y,'.k','MarkerSize',30)
set(gca,'XTick',x)set(gca,'XTickLabel','0|pi/4l3pi/8|3pi/4|pi')
set(gca,'XGrid','on')
axis([-0.25 3.5 -0.1 1.1])
hold off

نتیجه :

مثال شماره 17

منبع کدها :

کتاب MATLAB Guide - D.J. Higham - second edition - صفحه 163 از فایل pdf کتاب



مرتبط با: متلب ( MATLAB ) ,
طیف عدد موج

طیف عدد موج [۱] بیانگر سهم انرژی ادیهایی با اندازه های مختلف (در مقابل ادیهای با فرکانسهای مختلف برای طیف فرکانسی) در انرژی جنبشی اغتشاش می باشد. طیف عدد موج را نمی توان براحتی مستقیماً با استفاده از وسایل موجود اندازه گیری بدست آورد. در عمل آنرا با استفاده از تبدیل فوریه تابع ارتباز مکانی و با استفاده از فرضیه ی تیلور و طیف فرکانسی بدست می آورند.

طیف عدد موج یک و چند بعدی


در صورتی که اندازه گیری های لازم جهت محاسبه ی طیف مورد نظر در یک بعد انجام گیرند به آن طیف، طیف یک بعدی و در صورتیکه در سه بعد انجام گیرد به آن طیف سه بعدی می گویند. طیف فرکانسی یک طیف یک بعدی می باشد.

درمورد طیف یک بعدی عدد موج متاسفانه بایستی اذعان داشت که استفاده از آن مناسب نبوده و می تواند تا حدی گمراه کننده باشد. این امر بدلیل سه بعدی بودن جریان مغشوش در فضای مکانی (یا عدد موج) است.

به عنوان مثال طیف های یک بعدی عدد موج اندازه گیری شده معمولاً دارای مقادیر محدودی در مبدا مختصات (k = 0) می باشند که طبیعتاً نمی تواند به سهم انرژی ادیهای با عدد موج صفر تلقی گردد.این امر بدلیل پدیده ای به نام aliasing اتفاق می افتد.


شکل

توابع ارتباط مکانی اندازه گیری شده در یک جهت x نه تنها متاثر از اثر ادیهای با سایز k که در جهت x حرکت می کنند می باشند بلکه از ادیهایی با k بیشتر که در مایل نسبت به x حرکت می کنند نیز می باشند و لذا طیف انرژی در k در برگیرنده سهم انرژی ناشی از تمامی ادیهای با k های بزرگتر که در جهات دیگر حرکت می کنند نیز می باشد. مشکل aliasing در اعداد موج بالا جدی نبوده . چه اینکه ادیهای کوچک تقریباً دارای اندازه مساوی در تمامی جهات می باشند.(isotropy)

به ادامه ی مطلب مراجعه کنید

ادامه مطلب

تعداد کل صفحات: 26