今年的五一,没有离开厦门,原因是并不富裕,我渴望去旅游,而这需要经济基础.但也是值得铭记的一次,与舍友沿着山海步道,用脚步丈量了山海线的长度,从2点一直走到6点,近4w的步数,翻山越岭,应了那句话From the mountain to the sea.假期和舍友一起去酒吧喝酒,玩骰子,每个人都醉乎乎的走出酒吧,然后又去唱歌,一直到2点才缓缓返回学校.不出意外的被辅导员问话,好在语重心长,并无大碍.
Apple music在今年上新了独立的application，脱离了iTunes，有了自己的UI设计，在3月份，我并无期待的在PC端安装尝试了一下，只能说是效果一般。用也勉强能用，但最大的缺点是，无法看歌词。
于是，就找到了Lyricease，我愿称其为最好用，最好看的第三方网易云UMP
什么是UMP？
Universal Windows Platform(UWP)直译过来就是通用Windows平台，它是微软基于Win10开发的一个通用平台，微软希望自己所有的设备（包括Surface，Xbox，HoloLens等）都运行于同样的系统，然后让用户像使用苹果设备一样从自己官方的商城来下载安装应用，这样可以让应用市场中心化的同时保证系统的安全，稳定，高效。
clear
n = input('方程个数n=');
A = input('系数矩阵A=');
b = input('常数列向量b=');
if (size(A) ~= [n, n]) | size(B) ~=[n, 1]
disp('输入不正确，要求A是n阶方阵，b是n维列向量')
else if det(A) == 0
disp('系数行列式为0，不能用克拉默法则解此方程组')
else
for i = 1:n
B = A;
B(:,i) = b;
x(i) = det(B)/det(A);
end
end
end
clear
x1 = [1, 0, 2]';x2 = [2, 1, 1]';x3 = [2, 0,1]';x4 = [3, 1, 1]';
x5= [1, 1, 1]';
A = [x1, x2, x3, x4];
[R, jb] = rref(A); len = length(jb);
if len < 5
'The vector group is linearly dependent and serial numbers are'
jb
else
'The vector group is linearly independent'
end
运行结果；
ans =
The vector group is linearly dependent and serial numbers are
jb =
1 2 3
判定二次型的正定值。在Matlab中判定 $ f = 3x_1^2+2x_2^2+7x_3^2-4x_1x_2-4x_2x_3 $是否为正定二次型？
matlab代码：
clear;
A = [3, -2, 0;-2, 2, -2;0, -2, 7]; %输入二次型对应的对称阵
[q, d] = eig(A); %求矩阵A的特征值与特征向量
v = diag(d); %提取A的所有特征值
if all( v > 0) %判断矩阵A的特征值是否为正
disp('二次型为正定')
elseif all(v >= 0)
disp('二次型为半正定') %判断矩阵A的特征值是否均为非负
elseif all(v < 0)
disp('二次型为负定') %判断矩阵A的特征值是否均为负
elseif all(v <= 0)
disp('二次型为半负定') %判断矩阵A的特征值是否均为非正
else
disp('二次型为不定')
end
运行结果：
二次型为正定
6
计算二次型的正（负）惯性指数。在Matlab中计算二次型 $ f = x_1x_2+x_1x_3-3x_2x_3 $的正（负）惯性指数。
Matlab代码：
clear;
A = [0, 0.5, 0.5;0.5, 0, -1.5,;0.5, -1.5, 0];
%输入二次型对应的对称阵
[q, d] = eig(A); %计算矩阵A的特征值与特征向量
n = length(d); %计算矩阵A的阶数
zheng = 0; fu = 0; %将正负惯性指数初始值置为0
for i = 1:n
if d(i) > 0
zheng = zheng + 1;
elseif d(i) < 0
fu = fu + 1;
end
end
zheng, fu
a = int(input("输入的偏移量为:"))
b = input("输入单个大写英文字母为：")
c = (ord(b) - ord("A") + a) % 26
# 视A-Z为26个字母一个循环，取余数
d = chr(ord("A") + c)
print("经过凯撒加密之后，输出的字符为:" + d)
6. 使用 random.randint(a, b) 方法，随机生成三个100以内的自然数，求三个数的和。
year = int(input("输入年份："))
month = int(input("输入月份："))
day = int(input("输入日："))
m = month
if month < 3:
m = month + 12
year = year - 1
c = year // 100
y = year % 100
d = day
w = (y + (y // 4) + (c // 4) - 2 * c + (26 * (m + 1) // 10) + d - 1) % 7
if(w == 1):
message = "星期一"
elif (w == 2):
message = "星期二"
elif (w == 3):
message = "星期三"
elif (w == 4):
message = "星期四"
elif (w == 5):
message = "星期五"
elif (w == 6):
message = "星期六"
else:
message = "星期日"
print("当日为：{}".format(message))
或者通过列表实现。
year = int(input("输入年份："))
month = int(input("输入月份："))
day = int(input("输入日："))
m = month
if month < 3:
m = month + 12
year = year - 1
c = year // 100
y = year % 100
d = day
w = (y + (y // 4) + (c // 4) - 2 * c + (26 * (m + 1) // 10) + d - 1) % 7
w = int(w)
list = ['星期日', '星期一', '星期二', '星期三', '星期四', '星期五', '星期六']
print("当日为：", list[w])
import random
list1 = ['剪刀。', '石头。', '布。']
a = eval(input('剪刀（0），石头（1），布（2）：'))
b = random.randint(0, 2)
if a == b:
print('计算机是', list1[b], '你是', list1[a], '平局。')
elif a - b == 1 or a - b == -2:
print('计算机是', list1[b], '你是', list1[a], '你获胜了。')
else:
print('计算机是', list1[b], '你是', list1[a], '你输了。')
position = input("输入职位，经理或普通员工:")
order = eval(input("输入订单总额(单位：万元)："))
if position == "经理":
if order > 500:
bonus = order * 0.01
else:
bonus = order * 0.005
else:
if order > 50:
bonus = order * 0.02
else:
bonus = order * 0.01
print("该员工的职位为{}，绩效奖金为{}万元".format(position, bonus))
# 逻辑关系不太会处理，while也不熟悉，总之我太菜了，将就着看吧
b = 0
c = 0
sum = 0
count = 0
a = eval(input('Enter an integer, the input ends if it is 0:'))
if a != 0 :
sum = sum + a
if a > 0:
b += 1
elif a < 0:
c += 1
while a != 0 :
a = eval(input('Enter an integer, the input ends if it is 0:'))
count += 1
if a > 0:
b += 1
elif a < 0:
c += 1
sum = sum + a
average = sum / count
print('The number of positives is', b)
print('The number of negatives is', c)
print('The total is', sum)
print('The average is', average)
else:
print('You didn\'t enter anything')
17.输出100以内所有7的倍数。
outNum = 7
i = 1
while outNum < 100:
print(outNum, end=' ')
i = i + 1
outNum = 7 * i
#个人觉得这样处理逻辑关系更加清晰
print('模式A')
for i in range(1, 7):
for j in range(1, i + 1):
if i >= j:
print(j, end=' ')
print()
print('模式B')
for i in range(6, 0, -1):
for j in range(1, i + 1):
if i >= j:
print(j, end=' ')
print()
print('模式C')
for i in range(1, 7):
message = " " * (6 - i)
print(message, end='')
for j in range(6, 0, -1):
if i >= j:
print(j, end=' ')
print()
print('模式D')
for i in range(6, 0, -1):
message = " " * (6 - i)
print(message, end='')
for j in range(1, i + 1):
if i >= j:
print(j, end=' ')
print()
for i in range(1, 10000):
sum = 0
for j in range(1, i // 2 + 1):
if i % j == 0:
sum = sum + j
if sum == i:
print(i, end=" ")
或者你可以这样实现，考虑到运行效率，更建议用第一种逻辑关系。
perfectNum = []
getInput = 10000
for num in range(1, getInput + 1):
divisor = []
for i in range(1, num):
if (num / i) % 1 == 0:
divisor.append(i)
addDivisor = 0
for j in divisor:
addDivisor = addDivisor + j
if addDivisor == num:
perfectNum.append(num)
for k in perfectNum:
print(k, end=" ")
def isLeapYear(year):
if year % 4 == 0 and year % 100 != 0 or year % 400 == 0:
return True
else:
return False
year = int(input("输入一个年份："))
print(isLeapYear(year))
def isPrime(n):
if a <= 0:
return '请输入一个正整数'
else:
for i in range(2, a // 2 + 1):
if a % i == 0:
result = False
break
elif a % i != 0:
result = True
return result
a = int(input('请输入一个正整数：'))
print(isPrime(a))
def isPrime(number):
result = True
for factor in range(2, number // 2 + 1):
if number % factor == 0:
result = False
return result
getInput = int(input('请输入一个正整数：'))
print(isPrime(getInput))
def primeNumbers(a):
count = 0
# 设置从2到输入数字的循环
for num in range(2, a):
# 先默认判断其为素数
b = True
# 开始验证是否为素数
for i in range(2, num // 2 + 1):
# 经历循环，若能被其他数整除，则判断为False赋值给b
if num % i == 0:
b = False
# 当num为素数时候：
if b:
count += 1
print(num, end='\t')
# 执行换行
if count % 10 == 0:
print()
a = int(input('请输入一个正整数：'))
primeNumbers(a)
朋友对我说，他有点搞不懂我。嗯，不止一个人给我这样说，我是什么样的人？我也很难定义，Boldstroom中歌词这样写：“I am so complicated, but I really fucking meant it."挺合适的吧。很多时候，你觉得我并不是你觉得。就像你肯定不知道我在奔溃的时候的泪水。理性的时候极端理性，感性的时候极端感性。我自己也很难给自己下定义。这就是闫立晔，一个？的人。
今年的五一,没有离开厦门,原因是并不富裕,我渴望去旅游,而这需要经济基础.但也是值得铭记的一次,与舍友沿着山海步道,用脚步丈量了山海线的长度,从2点一直走到6点,近4w的步数,翻山越岭,应了那句话From the mountain to the sea.假期和舍友一起去酒吧喝酒,玩骰子,每个人都醉乎乎的走出酒吧,然后又去唱歌,一直到2点才缓缓返回学校.不出意外的被辅导员问话,好在语重心长,并无大碍.
Apple music在今年上新了独立的application，脱离了iTunes，有了自己的UI设计，在3月份，我并无期待的在PC端安装尝试了一下，只能说是效果一般。用也勉强能用，但最大的缺点是，无法看歌词。
于是，就找到了Lyricease，我愿称其为最好用，最好看的第三方网易云UMP
什么是UMP？
Universal Windows Platform(UWP)直译过来就是通用Windows平台，它是微软基于Win10开发的一个通用平台，微软希望自己所有的设备（包括Surface，Xbox，HoloLens等）都运行于同样的系统，然后让用户像使用苹果设备一样从自己官方的商城来下载安装应用，这样可以让应用市场中心化的同时保证系统的安全，稳定，高效。
clear
n = input('方程个数n=');
A = input('系数矩阵A=');
b = input('常数列向量b=');
if (size(A) ~= [n, n]) | size(B) ~=[n, 1]
disp('输入不正确，要求A是n阶方阵，b是n维列向量')
else if det(A) == 0
disp('系数行列式为0，不能用克拉默法则解此方程组')
else
for i = 1:n
B = A;
B(:,i) = b;
x(i) = det(B)/det(A);
end
end
end
clear
x1 = [1, 0, 2]';x2 = [2, 1, 1]';x3 = [2, 0,1]';x4 = [3, 1, 1]';
x5= [1, 1, 1]';
A = [x1, x2, x3, x4];
[R, jb] = rref(A); len = length(jb);
if len < 5
'The vector group is linearly dependent and serial numbers are'
jb
else
'The vector group is linearly independent'
end
运行结果；
ans =
The vector group is linearly dependent and serial numbers are
jb =
1 2 3
判定二次型的正定值。在Matlab中判定 $ f = 3x_1^2+2x_2^2+7x_3^2-4x_1x_2-4x_2x_3 $是否为正定二次型？
matlab代码：
clear;
A = [3, -2, 0;-2, 2, -2;0, -2, 7]; %输入二次型对应的对称阵
[q, d] = eig(A); %求矩阵A的特征值与特征向量
v = diag(d); %提取A的所有特征值
if all( v > 0) %判断矩阵A的特征值是否为正
disp('二次型为正定')
elseif all(v >= 0)
disp('二次型为半正定') %判断矩阵A的特征值是否均为非负
elseif all(v < 0)
disp('二次型为负定') %判断矩阵A的特征值是否均为负
elseif all(v <= 0)
disp('二次型为半负定') %判断矩阵A的特征值是否均为非正
else
disp('二次型为不定')
end
运行结果：
二次型为正定
6
计算二次型的正（负）惯性指数。在Matlab中计算二次型 $ f = x_1x_2+x_1x_3-3x_2x_3 $的正（负）惯性指数。
Matlab代码：
clear;
A = [0, 0.5, 0.5;0.5, 0, -1.5,;0.5, -1.5, 0];
%输入二次型对应的对称阵
[q, d] = eig(A); %计算矩阵A的特征值与特征向量
n = length(d); %计算矩阵A的阶数
zheng = 0; fu = 0; %将正负惯性指数初始值置为0
for i = 1:n
if d(i) > 0
zheng = zheng + 1;
elseif d(i) < 0
fu = fu + 1;
end
end
zheng, fu
