简单的线性回归问题-TensorFlow+MATLAB·

import tensorflow as tf
import numpy as np
import matplotlib.pyplot as plt
W = tf.Variable(tf.zeros([2, 1], name="weight_age"))
b = tf.Variable(0.0, name="bias")def inference(X):
return tf.matmul(X, W) + bdef loss(X, Y):
Y_predicted = inference(X)
return tf.reduce_sum(tf.square(Y-Y_predicted))def inputs():
weight_age = [[84,46],[93,20],[65,52],[70,30],[76,57],[69,25],[63,28],[72,36],[79,57],[75,44]
,[27,24],[89,31],[65,52],[57,23],[59,60],[69,48],[60,34],[79,51],[75,50],[82,34]
,[59,46],[67,23],[85,37],[55,40],[63,30]]
blood_fat_content = [354,190,405,263,451,302,288,385,402,365,209,290,346
,254,395,434,220,374,308,220,311,181,274,303,244]
return tf.to_float(weight_age), tf.to_float(blood_fat_content)def train(total_loss):
learning_rate = 0.00000001
return tf.train.GradientDescentOptimizer(learning_rate).minimize(total_loss)def evaluate(sess, X, Y):
print sess.run(inference([[80., 25.]]))
print sess.run(inference([[65., 25.]]))init = tf.global_variables_initializer()with tf.Session() as sess:
sess.run(init)
X, Y = inputs()
total_loss = loss(X, Y)
print "The total_loss:",total_loss
train_op = train(total_loss)
training_steps = 10000
for step in range(training_steps):
sess.run([train_op])
if step % 1000 == 0:
print "loss:", sess.run([total_loss])
evaluate(sess, X, Y)
print 'W is:',W.eval()
print 'b is:',b.eval()
sess.close()

% TensorFlow Data Plot Script For Matlab%
weight_age = [[84,46];[93,20];[65,52];[70,30];[76,57];[69,25];[63,28];[72,36];
[79,57];[75,44];[27,24];[89,31];[65,52];[57,23];[59,60];[69,48];[60,34];
[79,51];[75,50];[82,34];[59,46];[67,23];[85,37];[55,40];[63,30]]; % 25*2
weight = weight_age(:,1); %25*1
age = weight_age(:,2); %25*1
blood_fat_content = [354,190,405,263,451,302,288,385,402,365,209,290,346,254,395,434,220,374,308,220,311,181,274,303,244]; % 1*25
Fat = zeros(max(weight),max(age));
for i=1:25
Fat(weight(i),age(i)) = blood_fat_content(i);
end
figure,plot3(weight,age,blood_fat_content','.')
hold on% 平面拟合
A = 3.23492765;
B = 1.94148386;
C = 1.15458;
[x,y] = meshgrid(0:1:90);
z = A*x+B*y+C;mesh(z)