程式碼
clear all
close all
clc
a=imread('nomo.jpg'); % 輸入影像
[c1,c2,c3]=size(a);
b=rgb2ycbcr(a); % 轉換成ycbcr
cb=b(:,:,2);
cr=b(:,:,3);
cb(cb<=125)=255; % cb的判斷式,是人臉為白色
cb(cb~=255)=0;
figure(1);imshow(cb)
cr(cr>=137)=255; % cr的判斷式,是人臉為白色
cr(cr~=255)=0;
figure(2);imshow(cr)
rr=and(cr,cb); % 找出二者都是白的
rt(:,:,1)=double(rr).*double(a(:,:,1));
rt(:,:,2)=double(rr).*double(a(:,:,2));
rt(:,:,3)=double(rr).*double(a(:,:,3));
figure(3);imshow(uint8(rt))
要是沒先把膚色給分類出來並加以定位人臉的的位置,
則會影響接下來的辨識結果.
切割膚色有很多的分法,這個切割器我是使用ycbcr的方法來做為偵察.
為什麼使用ycbcr呢:
因為RGB的色彩空間很容易受到光源的影響而有明顯的變化,
而使在設定顏色的分割值,會有很大的大的誤差.
又因為ycbcr有良好的亮度和色度分離的特特(彩度)所以我就使用它.
在此Y 代表亮度,Cb和CR都代表色度(彩度) 轉換公式如下
但在matlab你不用建立上式,只要打rgb2ycbcr('檔案')轉換就可以了
%----程式碼-----
clear all
close all
clc
a=imread('nomo.jpg'); % 輸入影像
b=rgb2ycbcr(a); % 轉換成ycbcr
cb=b(:,:,2);
cr=b(:,:,3);
cb(cb<=125)=255; % cb的判斷式,是膚色為白色
cb(cb~=255)=0;
figure(1);imshow(cb)
cr(cr>=137)=255; % cr的判斷式,是膚色為白色
cr(cr~=255)=0;
figure(2);imshow(cr)
rr=and(cr,cb); % 找出二者都是白的
rt(:,:,1)=double(rr).*double(a(:,:,1));
rt(:,:,2)=double(rr).*double(a(:,:,2));
rt(:,:,3)=double(rr).*double(a(:,:,3));
figure(3);imshow(uint8(rt))
%--------------結果----------
原影像圖,如下
分割後的cb影像,如下
分割後的cr影像,如下
只秀出膚色圖
以上就是建立膚色切割器的基礎公式
由上方的圖可知,這些影像還一些小雜點,在這我用形態學中的Opening operation來消除雜點
雜訊去除
形態學
Opening(斷開)運算：
opening運算是先被侵蝕後再膨脹，通常是用來平滑影像輪廓，截斷窄的細頸，消除細的突支，
其原理如下:
（一）膨脹：將物體的邊界往外膨脹若干像素的方法。 如下圖：
（二）侵蝕：與膨脹相反的動作，這是一種從物體的邊界上，將物體往內收縮若干像素。
如下圖:
斷開圖如下
標記法：
兩像素間關係有下述兩種性質：
(一)像素的近鄰(neighbors)：
在座標(x，y)處的像素P有水平和垂直方向的近鄰，其座標為
(x-1，y)，(x+1，y)，(x，y+1)，(x，y-1)，
這樣一組像素稱P點的4-近鄰(4-nieghbors)。
P的四個對角近鄰的座標為(x-1，y+1)，(x+1，y+1)，(x-1，y-1)，(x+1，y-1)，
這些點與4-近鄰在一起，稱為P的8-近鄰(8-neighbors)。
（二）連通性(Connectivity)：
像素間的連通性是用於確定影像中物體邊界和區域組成成分的重要概念。
要判斷兩個像素是否連通，必須要確定它們是否在某種意義中相鄰(如它們是否為4-近鄰)，
並且它們的灰度值是否滿足特定的相似準則，(如灰度值是否相同)。
在此我們用V來定義連通性的灰階集合，並介紹二種類型的連通：
(a)4-連通：
若某一像素Q 位於4-近鄰( P)中，且Q和P的像素值位於所定的灰度值(V)中，則稱4-連通。
(b)8-連通：
若某一像素Q位於8-近鄰( P)中，且Q和P的像素值位於所定的灰度值(V)中，則稱8-連通。
(三)標記法原理:
標記法是使用在二值化影像，白色為255，黑色為0，p表掃描過程中之像素值；
先由左而右，再由上而下掃描影像，在此只考慮4-連通分量。
當掃描到為0的時則移動到下一掃描位置，當掃描到p為255的，並根據其連通性檢查p點的左方和上方，
如果它們均為0的為話，給p一個新的標記值，如果這二個近鄰中只有一個為255，
則把它的標記給予p，如果它們都是255且標記相同，則把這個標記給予p。
如果它們都是255但標記不同，則把其中一個標記給予p，並標明這兩個標記是等價的。
掃描結束後所有值為255的點都被標記過了，但有些標記是可能是等價的，
所以再對影像作第二次掃描，將所有被記錄為等價的標記值改成相同的標記。
新增加的程式碼
se=strel('disk',2);%用open過濾膚色圖像
bw= imopen(rr,se);
figure(4);imshow(bw)
%用open過濾後的圖,再用標記法來找出最大區塊
l=bwlabel(bw,8); % 標記法
stats=regionprops(l,'Area'); % 找出相同的並統計
allArea=[stats.Area];
tt=max(allArea); % 找出最大的區塊
idx=find(allArea==tt);
bw2=ismember(l,idx);
figure(5);imshow(bw2)%過濾後的人臉圖像
%----------------------------------------------------
GUIDE
首先把一張rgb的照片轉成ycbcr(y代表亮度,cb代表色度(blueness),cr也代表色度(redness))
然後我們可以使用slider來調整y cb cr的範圍,使它只出現膚色的範圍
即可達到膚色分割了
clear all
close all
clc
a=menu('萬年曆','建月曆請按此','查某年某月某日為星期幾請按此');
if a==1
y=inputdlg({'請輸入你要建那一年的月曆'});
y=str2num(char(y));
m=menu('請選擇你要建那一月的月曆','1','2','3','4','5','6','7','8','9','10','11','12');
ans=[];
if rem(y,100)==0 & rem(y,400)==0 rem(y,4)==0
d1=[31 29 31 30 31 30 31 31 30 31 30 31];
else
d1=[31 28 31 30 31 30 31 31 30 31 30 31];
end
s=y-1;c=d1(1:(m-1));
for k=1:d1(m)
c1=sum(c)+k;
ss=s+fix(s/4)-fix(s/100)+fix(s/400)+c1;
ans(k)=rem(ss,7);
end
ans(ans==0)=7;
pp=ans(1);ans1=[];
k1=0;
for n=1:6;
for m1=pp:7
ans1(n,m1)=k1;
k1=k1+1;
end
pp=1;
end
ans1(ans1>d1(m))=0;
s=xlswrite('月曆.xls',{[y],' 年',[m],' 月'});
s=xlswrite('月曆.xls',{'星期日','星期一','星期二','星期三','星期四','星期五','星期六'},...
'a2:g2');
s=xlswrite('月曆.xls',ans1,'a3:g8');
as=sprintf(['你的月曆已經建好了\n\n','請到你的儲存路徑下打開[月曆.xls]這個檔\n\n']);
msgbox(as)
else
y=inputdlg({'請輸入你要查詢那一年'});
y=str2num(char(y));
m=inputdlg({'請輸入你要查詢那一月'});
m=str2num(char(m));
k=inputdlg({'請輸入你要查詢那一日'});
k=str2num(char(k));
ans=[];
if rem(y,100)==0 & rem(y,400)==0 rem(y,4)==0
d1=[31 29 31 30 31 30 31 31 30 31 30 31];
else
d1=[31 28 31 30 31 30 31 31 30 31 30 31];
end
s=y-1;
c=d1(1:(m-1));
c1=sum(c)+k;
ss=s+fix(s/4)-fix(s/100)+fix(s/400)+c1;
dd=rem(ss,7)+1;
di={'日','一','二','三','四','五','六'};
msgbox(['你查詢的',num2str(y),'年',num2str(m),'月',num2str(k),'日','是星期',di(dd)])
end
第一題
有一向量A=3i +5j +10k和向量B=5i -6j +2k 都通過原點,請
找出兩者相交的單位向量
% 第一題程式碼 %
close all
clear all
clc
a1=[3 5 10];
b1=[5 -6 2];
c1=cross(a1,b1) % 在這利用cross即可求出
第一題結果
c1 =
70 44 -43
二者相交的向量為70i+44j-43k
第二題
利用A=randn(10,10)求出mean(A),median(A,2),std(A)
std(A,1,2)
% 第二題程式碼 %
a=randn(10,10);
disp('平均值')
m=mean(a)% 求平均值(在此是指每一行的平均)
mean(a(:))% 對整個矩陣做平均
disp('各列之中間值')
me=median(a,2)% 求各列之中間值(在此2是指每一列中間值,要是把2給去掉則是對每行運算)
disp('分母為n-1的標準差')
s=std(a)% 對行向量做標準差
s1=std(a(:))% 對整個矩陣做標準差
disp('分母為n的標準差,並對列向量來運算')
s1=std(a,1,2) % 在這1指分母為n的標準差,2指對列向量做標準差
第二題結果
平均值
m =(每一行為平均值)
-0.1643 0.3951 -0.6530 0.0288 0.0838 -0.1272 -0.2105 -0.7598 0.6833 -0.2402
ans =(整個矩陣的平均值)
-0.0964
各列之中間值
me =
-0.4912
-0.7159
-0.2354
0.2010
0.4158
0.3165
0.2103
-0.2451
-0.8788
-0.2415
分母為n-1的標準差
s =(對行向量做標準差)
0.9348 0.9574 0.8332 0.7537 1.3951 0.8999 1.0765 0.5953 0.7629 1.1109
s1 =(對整個矩陣做標準差)
1.0006
分母為n的標準差,並對列向量來運算
s1 =
0.6429
0.8450
1.0733
0.9973
0.6032
0.9794
1.2223
0.9861
0.7078
0.8956
第三題
Prove that exp(i*theta)=cos(theta)+i*sin(theta) using matlab commands
此題我利用判斷式來判斷此公式是否相等,當結果為1為相等,0為不相等
% 第三題程式碼 %
theta=[10 20 30 40 50 60]/180*pi;%把角度變弧度,在此我用10~60度下去判斷
x=((exp(i.*theta))==(cos(theta)+i.*sin(theta)))
第三題結果
x =
1 1 1 1 1 1
所以可證明此式是相等的
第四題
建立一個R=eye(3)*5+4*ones(3)i求其abs(R), angle(R), real(R) and imag(R)
% 第四題程式碼 %
R=eye(3)*5+(4*ones(3))*i % 為一向量所組矩陣
ab=abs(R)% 求?對值
an=angle(R).*180./pi% 求其夾角
re=real(R)% 求向量之實數部份
im=imag(R)% 求向量之虛數部份
第四題結果
R =
5.0000 + 4.0000i 0 + 4.0000i 0 + 4.0000i
0 + 4.0000i 5.0000 + 4.0000i 0 + 4.0000i
0 + 4.0000i 0 + 4.0000i 5.0000 + 4.0000i
ab =(?對值)
6.4031 4.0000 4.0000
4.0000 6.4031 4.0000
4.0000 4.0000 6.4031
an =(其夾角)
38.6598 90.0000 90.0000
90.0000 38.6598 90.0000
90.0000 90.0000 38.6598
re =(其實數)
5 0 0
0 5 0
0 0 5
im =(其虛數)
4 4 4
4 4 4
4 4 4
% 第一題程式碼 %
b=(reshape((1:50),5,10)).^2 % 可以直接把1:50放入,再直接開平方
% 第一題結果 %
b =
1 36 121 256 441 676 961 1296 1681 2116
4 49 144 289 484 729 1024 1369 1764 2209
9 64 169 324 529 784 1089 1444 1849 2304
16 81 196 361 576 841 1156 1521 1936 2401
25 100 225 400 625 900 1225 1600 2025 2500
% 第二題程式碼 %
a1=magic(5); % 把大於10的變NaN
a1(a1>10)=NaN
% 第三題程式碼 %
a1=magic(5); % 把大於10的變0
a1(a1>10)=0
% 第二題和第三題結果 %
a1 =
17 24 1 8 15
23 5 7 14 16
4 6 13 20 22
10 12 19 21 3
11 18 25 2 9
a1 =
NaN NaN 1 8 NaN
NaN 5 7 NaN NaN
4 6 NaN NaN NaN
10 NaN NaN NaN 3
NaN NaN NaN 2 9
a1 =
0 0 1 8 0
0 5 7 0 0
4 6 0 0 0
10 0 0 0 3
0 0 0 2 9
%第四題程式碼%
close all
clc
clear
x=input('請輸入你要抽的次數');
draw_number(x);
% function draw_number %
function [C]=draw_number(no_of_draw) % 因為要比較所以我用FOR把資料帶入
no_of_draw
o=length(no_of_draw);
% draw ball numbers within ndraw times
for m=1:o;
C=zeros(1,5);
n=1;
while n<=no_of_draw(m)
ball=fix(rand*10);
if ball<2
C(1)=C(1)+1;
elseif ball<4,
C(2)=C(2)+1;
elseif ball<6,
C(3)=C(3)+1;
elseif ball<8,
C(4)=C(4)+1;
else
C(5)=C(5)+1;
end
n=n+1;
end
fprintf('抽%d次的結果為\n',n-1)
disp('第一類 第二類 第三類 第四類 第五類')
C
end
% 第四題結果 %
請輸入你要抽的次數
[100:100:1000]
no_of_draw =
100 200 300 400 500 600 700 800 900 1000
抽100次的結果為
第一類 第二類 第三類 第四類 第五類
C =
25 16 20 15 24
抽200次的結果為
第一類 第二類 第三類 第四類 第五類
C =
38 41 42 42 37
抽300次的結果為
第一類 第二類 第三類 第四類 第五類
C =
55 62 51 59 73
抽400次的結果為
第一類 第二類 第三類 第四類 第五類
C =
77 76 97 76 74
抽500次的結果為
第一類 第二類 第三類 第四類 第五類
C =
95 102 110 106 87
抽600次的結果為
第一類 第二類 第三類 第四類 第五類
C =
117 120 120 118 125
抽700次的結果為
第一類 第二類 第三類 第四類 第五類
C =
145 162 122 144 127
抽800次的結果為
第一類 第二類 第三類 第四類 第五類
C =
157 151 172 169 151
抽900次的結果為
第一類 第二類 第三類 第四類 第五類
C =
176 191 187 161 185
抽1000次的結果為
第一類 第二類 第三類 第四類 第五類
C =
192 210 189 231 178
%第五題程式碼%
clear all
close all
clc
s=input('請輸入妳要的二位數號碼\n');
n1=1;
while n1<=20
m=0;
while 1
a1=fix(rand*100);%抽號碼
m=m+1;
if a1==s
fprintf('\n電腦抽出的結果%d\n',a1)
fprintf('電腦抽了第%d次才把你輸入的%d給抽出來',m,s)
% 花了多少次才把妳輸入的號碼給抽出來%
break
end
end
n1=n1+1;
fprintf('\n總共有20個循環,這是你第%d個循環\n\n',n1-1)%抽第幾次
end
% 第五題結果 %
請輸入妳要的二位數號碼
14
電腦抽出的結果14
電腦抽了第133次才把你輸入的14給抽出來
總共有20個循環,這是你第1個循環
電腦抽出的結果14
電腦抽了第506次才把你輸入的14給抽出來
總共有20個循環,這是你第2個循環
電腦抽出的結果14
電腦抽了第90次才把你輸入的14給抽出來
總共有20個循環,這是你第3個循環
電腦抽出的結果14
電腦抽了第22次才把你輸入的14給抽出來
總共有20個循環,這是你第4個循環
電腦抽出的結果14
電腦抽了第175次才把你輸入的14給抽出來
總共有20個循環,這是你第5個循環
電腦抽出的結果14
電腦抽了第44次才把你輸入的14給抽出來
總共有20個循環,這是你第6個循環
電腦抽出的結果14
電腦抽了第61次才把你輸入的14給抽出來
總共有20個循環,這是你第7個循環
電腦抽出的結果14
電腦抽了第33次才把你輸入的14給抽出來
總共有20個循環,這是你第8個循環
電腦抽出的結果14
電腦抽了第61次才把你輸入的14給抽出來
總共有20個循環,這是你第9個循環
電腦抽出的結果14
電腦抽了第38次才把你輸入的14給抽出來
總共有20個循環,這是你第10個循環
電腦抽出的結果14
電腦抽了第77次才把你輸入的14給抽出來
總共有20個循環,這是你第11個循環
電腦抽出的結果14
電腦抽了第1次才把你輸入的14給抽出來
總共有20個循環,這是你第12個循環
電腦抽出的結果14
電腦抽了第13次才把你輸入的14給抽出來
總共有20個循環,這是你第13個循環
電腦抽出的結果14
電腦抽了第48次才把你輸入的14給抽出來
總共有20個循環,這是你第14個循環
電腦抽出的結果14
電腦抽了第74次才把你輸入的14給抽出來
總共有20個循環,這是你第15個循環
電腦抽出的結果14
電腦抽了第11次才把你輸入的14給抽出來
總共有20個循環,這是你第16個循環
電腦抽出的結果14
電腦抽了第22次才把你輸入的14給抽出來
總共有20個循環,這是你第17個循環
電腦抽出的結果14
電腦抽了第84次才把你輸入的14給抽出來
總共有20個循環,這是你第18個循環
電腦抽出的結果14
電腦抽了第263次才把你輸入的14給抽出來
總共有20個循環,這是你第19個循環
電腦抽出的結果14
電腦抽了第153次才把你輸入的14給抽出來
總共有20個循環,這是你第20個循環
若採用定存的方式，每個月存放一定的數額（例如每個月五百元），銀行年息5%
close all
clear all
clc
n1=0;
sum=0;
while sum<10000000 sum="sum+sum*(5/100/12);" sum="sum+500;" n1="n1+1;" style="color: rgb(51, 204, 0);">Ans:每月存五佰元存到10,000,000要花89年
上個星期去淡水,感覺還真是不一樣啊~^__^...以前都搭 metro...
%第一小題程式碼%
close all
clear
clc
n=0;
sum=0;
d1=1;
while sum<100000000 n="n+1;%每一天增加一元" sum="sum+sum*(5/100/365);%每天的利息" sum="sum+n;%每天的本金" style="color: rgb(0, 153, 0);">Ans:每天增加一元存到10,000,000要花10726天%第二小題程式碼%
close all
clear
clc
n=0;
sum=0;
d1=1;
while sum<100000000
n=n+5;%每一天增加五元
sum=sum+sum*(5/100/365);%每天的利息
sum=sum+n;%每天的本金
end
fprintf('每天增加五元存到10,000,000要花%d天\n',n/5)
Ans:每天增加五元存到10,000,000要花5529天
%第一小題程式碼%
clear;
clc;
disp('三角函數表')
for n=0:90
fprintf('角度為%2d 正弦值為%0.4f 正切值為%0.4f',n,sind(n),tand(n))
fprintf('\n')
end
%結果%
三角函數表
角度為 0 正弦值為0.0000 正切值為0.0000
角度為 1 正弦值為0.0175 正切值為0.0175
角度為 2 正弦值為0.0349 正切值為0.0349
角度為 3 正弦值為0.0523 正切值為0.0524
角度為 4 正弦值為0.0698 正切值為0.0699
角度為 5 正弦值為0.0872 正切值為0.0875
角度為 6 正弦值為0.1045 正切值為0.1051
角度為 7 正弦值為0.1219 正切值為0.1228
角度為 8 正弦值為0.1392 正切值為0.1405
角度為 9 正弦值為0.1564 正切值為0.1584
角度為10 正弦值為0.1736 正切值為0.1763
角度為11 正弦值為0.1908 正切值為0.1944
角度為12 正弦值為0.2079 正切值為0.2126
角度為13 正弦值為0.2250 正切值為0.2309
角度為14 正弦值為0.2419 正切值為0.2493
角度為15 正弦值為0.2588 正切值為0.2679
角度為16 正弦值為0.2756 正切值為0.2867
角度為17 正弦值為0.2924 正切值為0.3057
角度為18 正弦值為0.3090 正切值為0.3249
角度為19 正弦值為0.3256 正切值為0.3443
角度為20 正弦值為0.3420 正切值為0.3640
角度為21 正弦值為0.3584 正切值為0.3839
角度為22 正弦值為0.3746 正切值為0.4040
角度為23 正弦值為0.3907 正切值為0.4245
角度為24 正弦值為0.4067 正切值為0.4452
角度為25 正弦值為0.4226 正切值為0.4663
角度為26 正弦值為0.4384 正切值為0.4877
角度為27 正弦值為0.4540 正切值為0.5095
角度為28 正弦值為0.4695 正切值為0.5317
角度為29 正弦值為0.4848 正切值為0.5543
角度為30 正弦值為0.5000 正切值為0.5774
角度為31 正弦值為0.5150 正切值為0.6009
角度為32 正弦值為0.5299 正切值為0.6249
角度為33 正弦值為0.5446 正切值為0.6494
角度為34 正弦值為0.5592 正切值為0.6745
角度為35 正弦值為0.5736 正切值為0.7002
角度為36 正弦值為0.5878 正切值為0.7265
角度為37 正弦值為0.6018 正切值為0.7536
角度為38 正弦值為0.6157 正切值為0.7813
角度為39 正弦值為0.6293 正切值為0.8098
角度為40 正弦值為0.6428 正切值為0.8391
角度為41 正弦值為0.6561 正切值為0.8693
角度為42 正弦值為0.6691 正切值為0.9004
角度為43 正弦值為0.6820 正切值為0.9325
角度為44 正弦值為0.6947 正切值為0.9657
角度為45 正弦值為0.7071 正切值為1.0000
角度為46 正弦值為0.7193 正切值為1.0355
角度為47 正弦值為0.7314 正切值為1.0724
角度為48 正弦值為0.7431 正切值為1.1106
角度為49 正弦值為0.7547 正切值為1.1504
角度為50 正弦值為0.7660 正切值為1.1918
角度為51 正弦值為0.7771 正切值為1.2349
角度為52 正弦值為0.7880 正切值為1.2799
角度為53 正弦值為0.7986 正切值為1.3270
角度為54 正弦值為0.8090 正切值為1.3764
角度為55 正弦值為0.8192 正切值為1.4281
角度為56 正弦值為0.8290 正切值為1.4826
角度為57 正弦值為0.8387 正切值為1.5399
角度為58 正弦值為0.8480 正切值為1.6003
角度為59 正弦值為0.8572 正切值為1.6643
角度為60 正弦值為0.8660 正切值為1.7321
角度為61 正弦值為0.8746 正切值為1.8040
角度為62 正弦值為0.8829 正切值為1.8807
角度為63 正弦值為0.8910 正切值為1.9626
角度為64 正弦值為0.8988 正切值為2.0503
角度為65 正弦值為0.9063 正切值為2.1445
角度為66 正弦值為0.9135 正切值為2.2460
角度為67 正弦值為0.9205 正切值為2.3559
角度為68 正弦值為0.9272 正切值為2.4751
角度為69 正弦值為0.9336 正切值為2.6051
角度為70 正弦值為0.9397 正切值為2.7475
角度為71 正弦值為0.9455 正切值為2.9042
角度為72 正弦值為0.9511 正切值為3.0777
角度為73 正弦值為0.9563 正切值為3.2709
角度為74 正弦值為0.9613 正切值為3.4874
角度為75 正弦值為0.9659 正切值為3.7321
角度為76 正弦值為0.9703 正切值為4.0108
角度為77 正弦值為0.9744 正切值為4.3315
角度為78 正弦值為0.9781 正切值為4.7046
角度為79 正弦值為0.9816 正切值為5.1446
角度為80 正弦值為0.9848 正切值為5.6713
角度為81 正弦值為0.9877 正切值為6.3138
角度為82 正弦值為0.9903 正切值為7.1154
角度為83 正弦值為0.9925 正切值為8.1443
角度為84 正弦值為0.9945 正切值為9.5144
角度為85 正弦值為0.9962 正切值為11.4301
角度為86 正弦值為0.9976 正切值為14.3007
角度為87 正弦值為0.9986 正切值為19.0811
角度為88 正弦值為0.9994 正切值為28.6363
角度為89 正弦值為0.9998 正切值為57.2900
角度為90 正弦值為1.0000 正切值為Inf
%第二小題程式碼%
clear;
clc;
disp('這是個角度轉三角函數的程式(sin 和 tan)');
% 以下是可以改變度數之範圍及小數點之位數之程式 %
--xy=input('請選擇您要的小數點位數:取小數點第四位請輸入1,取到小數點第十四位請輸入2\n');
--y=input('請輸入您要轉換的角度或角度範圍\n');
dg(xy,y) % 呼叫function dg %
%function dg%
function [sine,tane]=dg(xy,a);
if xy==1
----format short;
----for n=a
--------fprintf('角度為%2d 正弦值為%0.4f 正切值為%0.4f',n,sind(n),tand(n))
--------fprintf('\n')
----end
else
----format long;
----for n=a
--------fprintf('角度為%2d 正弦值為%0.14f 正切值為%0.14f',n,sind(n),tand(n))
--------fprintf('\n')
----end
end
%執行結果%
這是個角度轉三角函數的程式(sin 和 tan)
請選擇您要的小數點位數:取小數點第四位請輸入1,取到小數點第十四位請輸入2
1
請輸入您要轉換的角度或角度範圍
[0:5]
角度為 0 正弦值為0.0000 正切值為0.0000
角度為 1 正弦值為0.0175 正切值為0.0175
角度為 2 正弦值為0.0349 正切值為0.0349
角度為 3 正弦值為0.0523 正切值為0.0524
角度為 4 正弦值為0.0698 正切值為0.0699
角度為 5 正弦值為0.0872 正切值為0.0875
>>
這是個角度轉三角函數的程式(sin 和 tan)
請選擇您要的小數點位數:取小數點第四位請輸入1,取到小數點第十四位請輸入2
2
請輸入您要轉換的角度或角度範圍
[0:5]
角度為 0 正弦值為0.00000000000000 正切值為0.00000000000000
角度為 1 正弦值為0.01745240643728 正切值為0.01745506492822
角度為 2 正弦值為0.03489949670250 正切值為0.03492076949175
角度為 3 正弦值為0.05233595624294 正切值為0.05240777928304
角度為 4 正弦值為0.06975647374413 正切值為0.06992681194351
角度為 5 正弦值為0.08715574274766 正切值為0.08748866352592
%第一題程式碼%
close all; clear; clc disp('這是個?體運動的程式') x=input('請輸入速度\n'); y=input('請輸入仰角\n'); newton(x,y); %呼叫function%
%function newton%
function [h,r,y1]=newton(v,th);
sa=size(v);
g=9.81;
h=(v.^2.*(sind(th)).^2)./(2*g);
r=(v.^2.*sind(th.*2))./g;
for n=1:sa(1,2)
s=linspace(0,r(n));
y1=tand(th(n)).*s-(g./(2.*v(n).^2.*(cosd(th(n))).^2).*(s.^2));
plot(s,y1);
text(r(n)/2,h(n),sprintf('最高為%2.5f',h(n)))
text(r(n),0,sprintf('最遠處%2.3f',r(n)))
hold on
end
%執行結果%
%第二題和第三題程式碼%
close all;
clear;
clc
disp('這是個?體運動的程式')
x=input('請輸入速度\n');
y=input('請輸入仰角\n');
z=input('請輸入高度\n')
newton1(x,y,z); %呼叫Function newton1%
%function newton1%
function [h,r,y1]=newton1(v,th,h1);
sa=size(v);
g=9.81;
t=(v.*sind(th)+sqrt((v.*sind(th)).^2+(4.*0.5.*g.*h1)))./g
h=h1+(v.^2.*(sind(th)).^2)./(2*g);
r=v.*cosd(th).*t
for n=1:4
s=linspace(0,r(n));
y1=(tand(th(n)).*s-(g./(2.*v(n).^2.*(cosd(th(n))).^2).*(s.^2)))+h1;
plot(s,y1);
text(((v(n).^2.*sind(th(n).*2))/g)/2,h(n),sprintf('高%2.5fm',h(n)))
text(r(n),0,sprintf('%2.3fm遠\t\t\n ',r(n)))
hold on
end
%執行結果%
t(著地時間) =
5.0536 5.6038 6.3395 7.2327
