房地产价格与住房保障规模 数学建模(5)

y=x0(1,:); m1=size(y,1); z=x0(2:n,:);m2=size(z,1); for i=1:m1

for j=1:m2 t(j,:)=z(j,:)-y(i,:); end

m=min(min(abs(t'))); M=max(max(abs(t'))); p=0.5;

k=(m+p*M)./(abs(t)+p*M); rt=sum(k')/size(k,2); r(i,:)=rt; end r

程序二：

#include void main() {

double

a1[10]={98000.45,206068.65,325164.35,461338.35,620925.05,806733.65,1024256.35,1292020.05,1608248.85,1951713.55}; double

a2[10]={99214.55,208869.75,329202.45,465025.25,624903.55,809840.95,102655.35,1291965.65,1606011.05,1946517.95}; double a3[10]={1197,2133,3722,6190,7912,10475,13516,16819,20663,24690};

double 10]={100.4,201.1,300.3,401.5,505.4,607.2,708.7,813.5,919.4,1018.7};

double

a5[10]={12582,27883,45519,65536,89702,118481,153285,198907,252131,311653};

double

a6[10]={36.22,73.88,112.97,153.5,195.26,238.25,282.15,327.09,372.77,419.36};

double

a7[10]={542.4,1142,1731,2353,2959.4,3478.6,4175.4,5009.51,5992.15,7126.25};

double a8[10]={2063,4113,6225,8395,10645,13004,15782,18950,22415,26143};

for(int i1=0;i1<9;i1++) { long double temp1; temp1=(a1[i1]+a1[i1+1])/2; cout<

19

cout<

cout<

for(int i3=0;i3<9;i3++) {

long double temp3;

temp3=(a3[i3]+a3[i3+1])/2; cout<

cout<

for(int i4=0;i4<9;i4++) {

long double temp4;

temp4=(a4[i4]+a4[i4+1])/2; cout<

cout<

for(int i5=0;i5<9;i5++) {

long double temp5;

temp5=(a5[i5]+a5[i5+1])/2; cout<

cout<

for(int i6=0;i6<9;i6++) {

long double temp6;

temp6=(a6[i6]+a6[i6+1])/2; cout<

cout<

for(int i7=0;i7<9;i7++) {

long double temp7;

temp7=(a7[i7]+a7[i7+1])/2; cout<

cout<

for(int i8=0;i8<9;i8++) {

long double temp8;

temp8=(a8[i8]+a8[i8+1])/2; cout<

20

cout<

程序三：

#include void main() {

double a1[10]={98000.45, 108068.2, 119095.7, 135174.0, 159586.7, 185808.6, 217522.7, 267763.7, 316228.8, 343464.7};

double a2[10]={99214.55, 109655.2, 120332.7, 135822.8, 159878.3, 184937.4, 216314.4, 265810.3, 314045.4, 340506.9};

double a3[10]={1197, 1036, 1589, 2468, 1722, 2563, 3041, 3303, 3844, 4027}; double a4[10]={100.4, 100.7, 99.2, 101.2,103.9,101.8, 101.5, 104.8, 105.9, 99.3}; double a5[10]={12582, 15301, 17636, 20017, 24166, 28779, 34804, 45622, 54224, 59522};

double a6[10]={36.22, 37.66, 39.09, 40.53, 41.76, 42.99, 43.90, 44.94, 45.68, 46.59};

double a7[10]={542.4, 599.6, 589.0, 622.0, 606.4, 519.2, 696.8, 834.11, 982.64, 1134.1};

double a8[10]={2063, 2053, 2112, 2170, 2250, 2359, 2778, 3168, 3465, 3728};

for(int i1=0;i1<9;i1++) { long double temp1; temp1=a1[i1]/a1[i1+1]; cout<

21

}

{

long double temp4; temp4=a4[i4]/a4[i4+1]; cout<

cout<

for(int i5=0;i5<9;i5++) {

long double temp5; temp5=a5[i5]/a5[i5+1]; cout<

cout<

for(int i6=0;i6<9;i6++) {

long double temp6; temp6=a6[i6]/a6[i6+1]; cout<

cout<

for(int i7=0;i7<9;i7++) {

long double temp7; temp7=a7[i7]/a7[i7+1]; cout<

cout<

for(int i8=0;i8<9;i8++) {

long double temp8; temp8=a8[i8]/a8[i8+1]; cout<

cout<

程序四：

a=[-152.035 -265.617 -393.251 -541.132 -713.829 -915.495 -1158.14 -1450.13 -1779.98;1 1 1 1 1 1 1 1 1];

b=[-152.035 1;-265.617 1;-393.251 1;-541.132 1;-713.829 1;-915.495 1;-1158.14 1;-1450.13 1;-1779.982 1];

22

c=a*b; d=inv(c); e=d*a;

f=[108.0682; 119.0957; 135.1740; 159.5867; 185.8086; 217.5227; 267.7637; 316.2288; 343.4647]; g=e*f

程序五：

a=[-154.042 -269.036 -397.114 -544.964 -717.372 -456.248 -697.311 -1448.990 -1776.260;1 1 1 1 1 1 1 1 1];

b=[-154.042 1;-269.036 1;-397.114 1;-544.964 1;-717.372 1;-456.248 1;-697.311 1;-1448.990 1;-1776.260 1]; c=a*b; d=inv(c); e=d*a;

f=[109.6552; 120.3327; 135.8228; 159.8783; 184.9374; 216.3144; 265.8103; 314.0454; 340.5069]; g=e*f

程序六：

a=[-1665 -2927.5 -4956 -7051 -9193.5 -11995.5 -15167.5 -18741 -22676.5;1 1 1 1 1 1 1 1 1]

b=[-1665 1;-2927.5 1;-4956 1;-7051 1;-9193.5 1;-11995.5 1;-15167.5 1;-18741 1;-22676.5 1] c=a*b f=inv(c) g=f*a

h=[1036;1589;2468;1722;2563;3041;3303;3844;4027] i=g*h

程序七：

a=[-150.75 -250.7 -350.9 -453.45 -556.3 -657.95 -761.1 -866.45 -969.05;1 1 1 1 1 1 1 1 1]

b=[-150.75 1;-250.7 1;-350.9 1;-453.45 1;-556.3 1;-657.95 1;-761.1 1;-866.45 1;-969.05 1]

c=a*b

23