用matlab写一个小程序

% t4.m参考程序(考虑了多人最高、最低等情况)：

% 请看懂！！！能编最好！

function t4()

clc

%---------start 初始化---------------

row=input('输入学生人数：')

col=input('输入课程数目：'樱做敏)

score=fix(rand(row,col)*49)+49

xh0=1:row

kch0=1:col

for ii=1:row

zcj_pre(ii)=sum_(score(ii,:))

end

disp('原始成绩表<第一列为学号、中间依次为各科成绩、最后一列为学生总成绩>：')

disp([xh0' score zcj_pre'])

%---------end 初始化--------------------

% 调用tongji1

[c_max,index_max,c_min,index_min,c_ave]=tongji1(score)

disp('平均分统计：')

fprintf(' 课程 平均\n')

% 按以上格式输出

disp([kch0',round(c_ave)'])

disp('课程最高分统计：')

fprintf(' 课程 最高 学号\n')

for jj=1:col

fprintf('%d %d',kch0(jj),c_max(jj))

for kk=1:length(index_max{jj})

fprintf('%d ',index_max{jj}(kk))

end

fprintf('\n')

end

disp('课程最低分统计：')

fprintf(' 课程 最低 学号\n')

% 按以上格式输出

for jj=1:col

fprintf('%d %d',kch0(jj),c_min(jj))

for kk=1:length(index_min{jj})

fprintf('%d ',index_min{jj}(kk))

end

fprintf('\n')

end

% 调用tongji23

[s_max,index_max,s_min,index_min,px]=tongji23(score)

disp('学生总成绩最高分统计：')

fprintf('最高总分 学脊枝号\n')

% 按以上格式输出

fprintf('%d ',s_max)

disp(index_max)

disp('学生总成绩最低分统计：')

fprintf('最低总分 学号\n')

% 按以上格式输出

fprintf('%d ',s_min)

disp(index_min)

disp('--------------------------------')

fprintf('总成绩排降序：\n')

fprintf('总成绩 学号\n')

% 按以上格式输出

disp(px)

% 存入zcj

zcj=px(:,1)%zcj在主函数空间

% 存入xh

xh=px(:,2)%xh也在主函数空间

function [c_max,index_max,c_min,index_min,c_ave]=tongji1(score)

% c_max课程最高分col*1向量

% index_max课程最高分对应学号col*1元胞

% c_min课程最低分col*1向量

% index_min课程最低分对应学号col*1元胞

% c_ave课程平均成胡兆绩col*1向量

[row,col]=size(score)

for jj=1:col

[c_max(jj),index_max{jj}]=max_(score(:,jj))

[c_min(jj),index_min{jj}]=min_(score(:,jj))

c_ave(jj)=ave_(score(:,jj))

end

function [s_max,index_max,s_min,index_min,px]=tongji23(score)

% s_max col门课程总成绩最高分

% index_max col门课程总成绩最高分对应学号

% s_min col门课程总成绩最低分

% index_min col门课程总成绩最低分对应学号

% px 排降序后的矩阵row*2 [col门课程总成绩学号]

[row,col]=size(score)

for ii=1:row

zcj(ii)=sum_(score(ii,:))

end

[s_max,index_max]=max_(zcj)

[s_min,index_min]=min_(zcj)

px=sort_down(zcj)

function y=sum_(x)%自编求和函数，x为向量

y=0

for ii=1:length(x)

y=y+x(ii)

end

function y=ave_(x)%自编求平均值函数，x为向量

y=sum_(x)/length(x)

function [x_max,index]=max_(x)%自定义函数

% x_max：向量x中的最大值

% index：向量x中的最大值对应的下标向量（有可能多个）

index=[]

x_max=x(1)

for ii=1:length(x)

if x(ii)>=x_max

x_max=x(ii)

end

end

for ii=1:length(x)

if x(ii)==x_max

index=[index,ii]

end

end

function [x_min,index]=min_(x)%自定义函数

% x_min：向量x中的最小值

% index：向量x中的最小值对应的下标向量（有可能多个）

index=[]

x_min=x(1)

for ii=1:length(x)

if x(ii)<=x_min

x_min=x(ii)

end

end

for ii=1:length(x)

if x(ii)==x_min

index=[index,ii]

end

end

function px=sort_down(x)

% px=[y_down index]

% y_down:x排降序

% index:x排序后的行标序列

[row,col]=size(x)

if row==1&&col~=1%保证x为列向量

x=x'

end

row=length(x)

index=1:row

px=[x,index']%构成排序前px

%采用选择排序2：

for ii=1:row-1

kk=ii

for jj=ii+1:row

if px(jj,1)>px(kk,1)

kk=jj

end

end

if kk~=ii

temp=px(kk,:)

px(kk,:)=px(ii,:)

px(ii,:)=temp

end

end

%%

clc

clear

m = linspace(1,100,100)

flag = 0

res = 0

for i = 3:length(m)

for j = 2:i-1

if mod(m(i),j) == 0

flag = 1

break

else

flag = 0

end

end

if flag == 0

r = [res,m(i)]

res = r

end

end

r(2:length(r))

我做第二个题目

function pdemodel

[pde_fig,ax]=pdeinit

pdetool('appl_cb',1)

set(ax,'DataAspectRatio',[1 6.5 1])

set(ax,'PlotBoxAspectRatio',[15 10 1])

set(ax,'XLim',[-15 15])

set(ax,'YLim',[-10 120])

set(ax,'XTickMode','auto')

set(ax,'YTickMode','auto')

% Geometry description:

pderect([0 1 100 0],'R1')

set(findobj(get(pde_fig,'Children'),'Tag','PDEEval'),'String','R1')

% Boundary conditions:

pdetool('changemode',0)

pdesetbd(4,'neu',1,'0','0')

pdesetbd(3,'dir',1,'1','0')

pdesetbd(2,'neu',1,'0','0')

pdesetbd(1,'dir',1,'1','100')

% Mesh generation:

setuprop(pde_fig,'Hgrad',1.3)

setuprop(pde_fig,'refinemethod','regular')

pdetool('initmesh')

pdetool('refine'宽扒)

pdetool('refine')

pdetool('refine')

pdetool('refine')

% PDE coefficients:

pdeseteq(1,'1.0','0.0','0','1.0','0:10','0.0','0.0','[0 100]')

setuprop(pde_fig,'currparam',['1.0''0.0''0 ''1.0'])

% Solve parameters:

setuprop(pde_fig,'solveparam',...

str2mat('0','18432','10','pdeadworst',...

'0.5'慎唯昌,'longest','0','1E-4','','fixed'山唯,'Inf'))

% Plotflags and user data strings:

setuprop(pde_fig,'plotflags',[1 1 1 1 1 1 1 1 0 0 0 1 1 1 0 0 0 1])

setuprop(pde_fig,'colstring','')

setuprop(pde_fig,'arrowstring','')

setuprop(pde_fig,'deformstring','')

setuprop(pde_fig,'heightstring','')

% Solve PDE:

pdetool('solve')

---