matlab 编程

第一题：

a=[2,5,4,6,3,78,1,6,5,21];

for j=1:9

for i=1:(10-j)

if a(i)>a(i+1)

b=a(i+1);

a(i+1)=a(i);

a(i)=b;

end

end

end

第二题

a=[5,2,3,6,4;7,8,9,6,12;11,12,13,14,22;44,41,53,26,31];

t=max(max(a))

for i=1:4

for j=1:5

if t==a(i,j)

i

j

end

end

end

第三题

s=1;

sum=0;

for i=1:20

s=si;

sum=sum+s;

end

sum

第四题

sum=0;

s=100;

%下降部分

for i=1:10

sum=sum+s

s=s/2;

end

%上升部分

s=50;

for i=1:9

sum=sum+s;

s=s/2;

end

sum

第五题

x=input('请输入数字：');

y=input('请输入数字：');

f=x^2+sin(xy)+2y

希望对你有些帮助。

1、用两种方法求解方程组的解

2x+y+z=10

3x+1y+z=14

2x+3y-1z=1

（注：a1代表学号倒数第一位数；a2代表学号倒数第二位数；a3代表学号倒数第三位数）

2、对函数y=[sin(112x+112)]/(113x)进行如下 *** 作：（1）求x趋近112时的极限值（2）求出x=112时的二阶导；（3）求从112到113的定积分

（注：a为自己的学号后三位数；b=a+1）

3、画图题:利用ezplot函数将以下y1=e112x;y2=x112;y3=112x函数的图像画在一个figure图中，并且给该figure图增加标题为“112同学的函数”，每个子图添加函数表达式的标注，横坐标为（-112,112）

（注：题中所有的a表示学号的后三位数

）

将前3题保存为以姓名拼音的M脚本文件

4、程序设计题

编写一个函数文件实现如下功能：输入一自然数x，若x为偶数，则分别计算出s1=1、s2=1+2、s3=1+2+3…sx=1+2+…+x，并把1到x按序存入向量r中，s1到sx存入相量s中，然后利用plot函数画出r与s的数学关系曲线；若x为奇数，则分别计算出s1=1/1、s2=(1+2)/2、s3=(1+2+3)/3…sx=(1+2+…+x)/x，并把1到x按序存入向量r中，s1到sx存入相量s中，然后利用plot函数画出r与s的数学关系曲线

[注：每个同学的曲线颜色、线型都不可相同，对于每个同学的要求在上机时候提出来]

5、电路仿真题

电阻的单位均为欧姆，测各支路电流波形，各元件电压波形，并用用显示器display显示其电流有效值

这里的关键在于mat2gray这个函数的用法。

mat2gray有两种调用方式：

I = mat2gray(A, [amin amax])

I = mat2gray(A)

第一种调用方式指定对应于光强0-1的A范围，而第二种方式则以A的最小值、最大值分别作为amin和amax。这样意味着，按照正常的理解，如果A的最小值不是0，本来没有完全黑色的点的，但在第二种调用方式下，最小值的点光强为0，即为全黑，类似的，最大值为完全白色。

如果用

H=mat2gray(a,[0 255]);

figure; imshow(H);

会得到和原始的a相同的图像。

%设计低通滤波器：

[N,Wc]=buttord()

%估算得到Butterworth低通滤波器的最小阶数N和3dB截止频率Wc

[a,b]=butter(N,Wc); %设计Butterworth低通滤波器

[h,f]=freqz(); %求数字低通滤波器的频率响应

figure(2); % 打开窗口2

subplot(221); %图形显示分割窗口

plot(f,abs(h)); %绘制Butterworth低通滤波器的幅频响应图

title(巴氏低通滤波器'');

grid; %绘制带网格的图像

sf=filter(a,b,s); %叠加函数S经过低通滤波器以后的新函数

subplot(222);

plot(t,sf); %绘制叠加函数S经过低通滤波器以后的时域图形

xlabel('时间 (seconds)');

ylabel('时间按幅度');

SF=fft(sf,256); %对叠加函数S经过低通滤波器以后的新函数进行256点的基—2快速傅立叶变换

w= %新信号角频率

subplot(223);

plot()); %绘制叠加函数S经过低通滤波器以后的频谱图

title('低通滤波后的频谱图');

matlab 程序设计

>

你可以参考下

>

我做第二个题目

function pdemodel

[pde_fig,ax]=pdeinit;

pdetool('appl_cb',1);

set(ax,'DataAspectRatio',[1 65 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',13);

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

pdetool('initmesh')

pdetool('refine')

pdetool('refine')

pdetool('refine')

pdetool('refine')

% PDE coefficients:

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

setuprop(pde_fig,'currparam',['10';'00';'0 ';'10'])

% Solve parameters:

setuprop(pde_fig,'solveparam',

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

'05','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')

以上就是关于matlab 编程全部的内容，包括:matlab 编程、matlab函数程序设计问题、matlab程序设计等相关内容解答，如果想了解更多相关内容，可以关注我们，你们的支持是我们更新的动力！