浙江大学远程教育学院。
2023年夏《可视化科学计算》上机操作。
1、 调试下列程序：
h0=figure('*******','none','position',[450 280 370 230],.
'name','图形对象属性的获取和设置实例');
t=0:pi/40:4*pi; hline1=plot(t,sin(t));
xlabel('t');ylabel('sin(t)')grid; ht=title('正弦曲线');
hline2=line(t+0.2,sin(t),'color','r');
haxes=gca;
h0_color=get(h0,'color')
haxes_color=get(haxes,'color')
hline1_color=get(hline1,'color')
hline2_color=get(hline2,'color')
set(haxes,'gridlinestyle
set(hline2,'linewidth',4,'color',[0.8 0.8 0.8]);
set(ht,'fontangle','italic','fontweight','bold');
set(gca,'children',[hline1 hline2
1. 写出程序的运行结果；
h0_color=0.8000 0.8000 0.8000
haxes_color=1 1 1
hline1_color=0 0 1
hline2_color=1 0 0
2. 修改程序，程序的标题为你的姓名，把正弦线由你学号顺序构成的各个位数的折线。（本题10分）
h0=figure('*******','none','position',‘7 1 3 0 0 9 0 1 2 0 0 1]’,张德军','图形对象属性的获取和设置实例');
t=0:pi:4*pi;
hline1=plot(t,sin(t));
xlabel('t');
ylabel('sin(t)')
grid;
ht=hline('713009012001');
haxes=gca;
2、 在同一画面中应用函数subplot(m, n, p)按上、下分别画出y=sin(2x),y=sin(x)*cos(2x)的图形，图形的标题是你单位的全称，横坐标中显示的是你的姓名，纵坐标是你的学号。（本题10分）。
zhdj_subplot(2,1,1)
ezplot('sin(2*x)')
title('兴达数码快印')
xlabel('张德军');
ylabel('713009012001')
subplot(2,1,2)
ezplot('sin(x)*cos(2*x)')
title('兴达数码快印')
xlabel('张德军');
ylabel('713009012001')
3、有一质点，初速度为v0,加速度为a，经过t时间后，质点所运动的位移。请写出matlab源**。 （本题10分）。
zhdj_clc;clear all;close all ;
v0=input('请输入初始速度：')
a=input('请输入加速度：')
t=0:0.1:5;
s=v0*t+0.5*a*t.^2;
plot(t,s)
xlabel('t/s')
ylabel('s/m')
4、有下列实验数据：
用三次多项式拟合，求出表达式的值，画出原数值图与拟合示意图，图的标题设置为你的学号，并测算x为时的y值（本题10分）。
zhdj_clc;clear all;close all;
x=1:9;
y=[10.2 14.1 13.9 12.3 10.5 11 12.4 13.5 16];
p=polyfit(x,y,3);
plot(x,y,'r')
y1=polyval(p,x);
hold on
plot(x,y1,'y')
y2=polyval(p,[7.5 12]) 当x为7.5,12时的y值。
、设计一程序，注意程序的流程是选择一张有你本人的**（你满意的），用matlab语句显示**后，将**转化为灰度图像，然后进行灰度变换，最后提取**的轮廓线。（本题15分）。
zhdj_clc;clear all;close all;
a1=imread('d:\'
imshow(a1)
a2=rgb2gray(a1);
a3=edge(a2);
imshow(a3)
、本题按以下步骤进行：
1）对信号y=sin(314t)+sin(628t)加上白噪声，画出图像。
2）进行离散傅立叶变换，并画出傅立叶变换后的频谱图、功率谱图。
要求图像横坐标为你的学号，纵坐标为你的地址，标题为你的姓名（本题15分）。
clearfs = 100采样率。
t = 0:1/fs:10-1/fs时间轴。
y = sin(314*t)+sin(628*t信号。
snr = 5
y = awgn(x,snr
figure,plot(t,y1绘制信号波形。
xlabel('713009012001');
ylabel('金华市婺城区北苑小区');
title('张德军');
f = 1/(0.005*pi定义采样频率。
n = 1200; %定义采样点数。
t = 0:n-1)./f定义时间序列。
y = sin(314*t) +sin(628*t信号表达式。
plot(t,y); 绘制信号时间图像。
grid on; %显示网格。
title('张德军信号时间图像。
xlabel('713009012001横坐标标注。
ylabel('金华市婺城区北苑小区纵坐标标注。
y1=fft(y进行快速傅里叶变换。
magy1=2*abs(y1)./n提取变换幅值并统一单位。
f=(0:n-1)*f./n定义频率序列。
plot(f,magy1绘制频幅图。
grid on; %网格。
title('张德军信号频谱图像。
xlabel('713009012001横坐标标注。
ylabel('金华市婺城区北苑小区纵坐标标注。
、结合工作单位的可视化计算的应用情况，自己命题应用matlab中所学的知识，解决一个实际问题。要求写出题目、题意分析、matlab指令或程序、最后答案及分析（本题30分）。
1）选题理由。
工资问题，对于每个人都是个不小的问题，2）题目。
工资计算。3）题意分析。
通过公司的营业额，来计算出每个人的最优工资数。
4）、matlab源**。
clear;clc
fora=30000:1000:150000;%a从30000每隔1000到150000循环。
x=0:0.1:1600;%x的范围。
y1=zeros(size(x));
y2=zeros(size(x));
y3=zeros(size(x));
y4=zeros(size(x));
n=length(x);
fork=1:n
ifx(k)>=a-151600)/11&&x(k)<(a-61600)/11
y1(k)=-2.2*x(k)+0.2*a-775;
elseifx(k)>=a-61600)/11&&x(k)<(a-25600)/11
y2(k)=-1.65*x(k)+0.15*a-425;
elseifx(k)>=a-25600)/11&&x(k)<(a-7600)/11
y3(k)=-1.1*x(k)+0.1*a-225;
elseifx(k)>=a-7600)/11&&x(k)<=a-1600)/11
y4(k)=-0.55*x(k)+0.05*a-100;
endend
y=y1+y2+y3+y4;%y表示f(x)
ymin,t]=min(y);%求f(x)的最小值ymin，t表示其位置。
ymin;x(t);%每月工资最优分配额。
endymin
x(t) 5）结果。
yminx(t)ymin=
ans=
6）讨论。这个程序，可以算出，发的最优工资数，当然，实际发工资的时候，并不一定能按这个，我相信对于很多公司可能适用。
可视化计算离线作业 1
浙江大学远程教育学院。可视化计算 课程作业2015年 冬 注意 所有图像的标题必须呈现足够你本人信息。1.10分 求解下列线性方程组的解 答案 a 1 1.5 2 9 7 0 3.6 0.5 4 4 7 10 3 22 33 3 7 8.5 21 6 3 8 0 90 20 b 3 4 20 5 1...
《科学计算可视化》作业 二 答案
科学计算可视化 作业 二 答案。1 科学计算可视化 visualizationinscientificcomputing,visc 是发达国家80年代后期提。出并发展起来的一个新的研究领域。1987年2月，美国国家科学 会在华盛顿召开了有关科学计算可视化的首次会议，与会者认为 将图形图像技术应用于科...
可视化编程作业
作业。课程名称软件工程 学生学院计算机学院 专业班级 10级4班 学号 3110006379 学生姓名陈泳鑫 指导教师崔洪刚 2013 年 1 月 6 日。题目 住院管理系统。1 角色分析 用户对象模型。用户对象有 操作员 系统管理员。医院住院管理系统的用例图。用例图说明 密码修改 用于修改当前用户...