وبلاگ تخصصی متلب و نقشه برداری

برنامه درونیابی محلی به روش second exact fit در متلب

سه شنبه, ۲۶ آبان ۱۳۹۴، ۱۱:۳۰ ق.ظ

برنامه را در mfile کپی کنید

%daroonyabi mahalli-Tin Based-second exact fit%
%zP=?%
clc
clear all
format long g
a=input('x,y,z a ra be soorat [x y z] vared konid=');
b=input('x,y,z b ra be soorat [x y z] vared konid=');
c=input('x,y,z c ra be soorat [x y z] vared konid=');
d=input('x,y,z d ra be soorat [x y z] vared konid=');
e=input('x,y,z e ra be soorat [x y z] vared konid=');
f=input('x,y,z f ra be soorat [x y z] vared konid=');
p=input('x,y P ra be soorat [x y] vared konid=');
ma=[1 a(1,1) a(1,2) a(1,1)^2 a(1,2)^2 a(1,1)*a(1,2)
1 b(1,1) b(1,2) b(1,1)^2 b(1,2)^2 b(1,1)*b(1,2)
1 c(1,1) c(1,2) c(1,1)^2 c(1,2)^2 c(1,1)*c(1,2)
1 d(1,1) d(1,2) d(1,1)^2 d(1,2)^2 d(1,1)*d(1,2)
1 e(1,1) e(1,2) e(1,1)^2 e(1,2)^2 e(1,1)*e(1,2)
1 f(1,1) f(1,2) f(1,1)^2 f(1,2)^2 f(1,1)*f(1,2)];
l=[a(1,3);b(1,3);c(1,3);d(1,3);e(1,3);f(1,3)];
x=(ma'*ma)^-1*ma'*l;
zp=x(1,1)+(x(2,1)*p(1,1))+(x(3,1)*p(1,2))+(x(4,1)*(p(1,1))^2)+(x(5,1)*(p(1,2))^2)+(x(6,1)*p(1,1)*p(1,2))

نظرات  (۰)

هیچ نظری هنوز ثبت نشده است

ارسال نظر

ارسال نظر آزاد است، اما اگر قبلا در بیان ثبت نام کرده اید می توانید ابتدا وارد شوید.
شما میتوانید از این تگهای html استفاده کنید:
<b> یا <strong>، <em> یا <i>، <u>، <strike> یا <s>، <sup>، <sub>، <blockquote>، <code>، <pre>، <hr>، <br>، <p>، <a href="" title="">، <span style="">، <div align="">
تجدید کد امنیتی