浙江大学继续教育学院
2012年《可视化计算》上机操作
姓名学号得分
１、在同一窗口中画出y=sin(2x),y=sin(x)*cos(2x)的图形（本题10分）。
>> t=linspace(0,2*pi,36);
>> x=sin(2*t);
>> plot(t,x);
>> y=sin(t).*cos(2*t);
>> hold on;
>>plot(t,y);
２、求矩阵A的逆矩阵，其中：A=[1 2 3;4 5 6;7 8 9] （本题10分）。
３、请写出Matlab指令或程序求出下列方程组的解（本题10分）。
6x1+x2+3x3=6
3x1+5x2+x3=1
x1+6x2+4x3=8
４、求函数x=*t+4在区间t[0,5]的数值积分值，写出matlab程序代码（本题10分）。
５、有下列实验数据：
用三次多项式拟合，求出表达式的值，画出拟合示意图，并测算x 为11、12时的y值（本题10分）。
x=[1 2 3 4 5 6 7 8 9];
y=[ 11 16];
y0=polyfit(y,x,3)
x1=11;
x2=12;
f = polyval(y0,y);
plot(y,x,'o',y,y0,'-r')
f1 = polyval(y0,x1))
f2= polyval(y0,x2))
６、求二元高次方程组x2+3y=4 的解析解（本题10分）
y2+4x=8
代码：
f1='x^2+3*y=4';
f2='y^2+4*x=8';
[x,y]=solve(f1,f2,'x','y');
x=vpa(x,4)
y=vpa(y,4)
图：
７、设计一程序，先将彩色图片转换为黑白图片，然后再转换为二进位影像（本题10分）。
im=imread('');
im=rgb2gray(im);
imshow(im);
im=double(im);
８、对信号y=sin(314t)+sin(628t)进行离散傅立叶变换，并画出它们的图像（本题10分）。
f = 1/*pi);
N = 1200; t = (0:N-1)./f ; y = sin(314*t) + sin(628*t) ;
plot(t,y); grid on;
９、结合工作单位的可视化计算的应用情况，自己命题应用matlab
中所学的知识，解决一个实际问题。要求写出题目、题意分析、Matlab
指令或程序、最后答案及分析（本题20分）。
题目：
随机产生两个3X4的矩阵A与B, 求所有同一位置上的较大元素构成
的新矩阵C,并得到矩阵C中的最大值的行与列
代码：
a=rand(3,4),b=rand(3,4)
for i=1:3
for j=1:4
c(i,j)=max(a(i,j),b(i,j)); endendc
[x,y]=find(c==max(c(:)))图：
附：上机考核范围（2012年-2015年1月间下列内容不变）：
1.应用matlab作图（二维、三维、标题、标注、坐标标记等）
2.矩阵运算（矩阵加、减、乘、逆）
3.线性方程组求解（线性方程solve与非线性方程fsolve）
4.数值积分（数值积分int与微分diff）
5.曲线插值与拟合(polyfit、polyval)
6.符号运算（符号矩阵、因式分解、极限）
7.图像载入与基础变换（旋转、影像显示在圆柱体、球体、黑白影
像）
8.傅里叶频谱变换（白噪声产生、信号迭加运算、付里叶变换及逆
变换、功率谱）
9.实际工程应用
